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INTRODUCTION 


In nova fert animus mutatas dicere formas 
corpora; di, coeptis (nam vos mutastis et illas) 
adspirate meis primaque ab origine mundi 
ad mea perpetuum deducite tempora carmen! 

—Publius Ovidius Naso [Metamorphoses] 

My soul is wrought to sing of forms transformed to bodies new and strange! 

Immortal Gods Inspire my heart, for ye have changed yourselves and all 
things you have changed! Oh lead my song in smooth and measured strains, 
from olden days when earth began to this completed time! 

—Ovid, Metamorphoses [Translation by Brookes More, Boston, 1922.] 

Metamorphoses are performed on the way from the visible world as it ap¬ 
pears to human eyes down to what we call the “fundamental” level. Fun¬ 
damental physics follows very simple and very clear rules. The rules of our 
fundamental laws of nature are symmetry and breaking of symmetries. 

Spontaneous symmetry breaking was shovra to occur also at the very 
fundamental level confirmed by the discovery of the Higgs boson at the 
Large Hadron Collider (LHC) close to Geneva. The Standard Model of el¬ 
ementary particle physics has finished its triumphal procession and shall 
be completed. However—the Standard Model is not the ultimate truth. 
Not only that there are observations in nature that cannot be explained 
within the Standard Model: about 95 % of our universe seems to be un¬ 
known and there is no sufficient explanation of why we live in a universe 
of matter. The missing pieces seem to be triggered by cosmology. There 
are still conceptual puzzles that lack an explanation: Why are there three 
families of matter fermions and why do they mix so strangely? And what 
causes electroweak symmetry to be broken? A further open issue can be 
solved by a minimal and somewhat symmetric extension of the Standard 
Models: neutrino masses. The question why they are so small still remains 
unsettled. 

We know since the invention of Quantum Mechanics that we live in 
a quantum world. Quantum corrections are of importance in any theo¬ 
retical description of our fundamental processes. Precision calculations 
are performed and applied to collider phenomena at very high accuracy. 
Quantum Electrodynamics was tested with an amazing precision. Quan¬ 
tum corrections also guide us through the following pages: we give an 
explanation of neutrino mixing by virtue of threshold corrections as they 
may arise in a popular extension of the Standard Model. For that pur¬ 
pose, we start with degenerate neutrino masses at leading order and get 
the observed deviations from that degenerate pattern by quantum cor¬ 
rections which also generate the mixing. The validity of this description 
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INTRODUCTION 


Thanks to the 
ClassicThesis style. 


can be excluded with an upcoming measurement of the neutrino mass: 
are rather heavy neutrinos excluded by experiment, it gets implausible 
for such quantum corrections being significant. In a concluding chapter, 
we give an explanation of large neutrino mixing for a hierarchical mass 
spectrum in contrast. 

Quantum corrections in those extensions are known to have an impact 
on electroweak symmetry breaking. We investigated their meaning for the 
parameter range which is favored by quantum models of fermion mixing 
and found a genuine effect of such corrections which has not yet been 
described in the literature. The results can then be used to constrain the 
parameter range from the requirement of the electroweak vacuum being 
stable and are complementary to existing constraints. 

This thesis is structured as follows: We start with an overview of the 
fundamentals of modern particle physics in Chapter 2 and set up the 
Standard Model and its popular extensions we deal with throughout this 
thesis. A special focus lies on supersymmetric extensions and their ap¬ 
plication on neutrino flavor physics that shall be discussed in Chapter 3. 
There we propose an explanation for the large observed neutrino mix¬ 
ing based on quantum corrections without referring to tree-level flavor 
models which is quite orthogonal to what is mostly performed in modern 
literature. Quantum corrections also possibly destabilize the electroweak 
vacuum: in Chapter 4 we first review the status of one-loop corrections 
to the effective Higgs potential and explicitly calculate such corrections 
for the dominating parts of the Minimal Supersymmetric Standard Model. 
Transferring this knowledge to the neutrino extension, we show in Chap¬ 
ter 5 that there is no fear of vacuum decay in presence of heavy Majorana 
neutrinos. Finally, we comment on a further possibility to explain large 
neutrino (and simultaneously small quark) mixing exploiting the hierar¬ 
chical mass patterns in Chapter 6. 

For the introductory part, we do not intent to give a complete and ex¬ 
haustive overview of the field. This is simply not possible in a limited 
amount of space. Nevertheless, we try to give as much information as pos¬ 
sible and elaborate some theoretical basement as necessity for the work 
performed within this thesis. The major achievements as outcome of the 
following pages is the attempt of a quantum corrected description of neu¬ 
trino mixing in view of any unknovra new physics (where as example 
of known new physics we take supersymmetry) and the influence of the 
quantum nature of our basic theory on the stability of the ground state of 
the theory. Notations are explained where they appear for the first time. 
Matrices in flavor space are denoted with bold symbols, m. Upshape labels 
are put wherever they fit, so e. g. m^. = nix ^ 

There are marginal notes appearing rather continuously. Their mean¬ 
ing is mostly to be seen as side remarks, where from time to time some 
commonly known and used notation is introduced for completeness. Gen¬ 
erally, this thesis shall be readable and understandable also while ignoring 
the page margins. 


PARTICLE PHYSICS IN A NUTSHELL 



2.1 THE STANDARD MODEL OF ELEMENTARY PARTICLES 

The Standard Model (SM) of elementary particles describes fundamental 
interactions of the smallest constituents of matter. We call elementary par¬ 
ticles elementary, because no substructure has been observed yet and they 
are sufficient to build up the matter in the universe. Actually, only about 
five percent of the matter content of our universe is made out of what is de¬ 
scribed by the SM, see e. g. [1]. The unknown matter species, “Dark Mat¬ 
ter”, only interacts gravitationally with our known matter, maybe there are 
weak gauge interactions [2, 3]. Gravity itself is not included in the SM of 
elementary particles, all what will be discussed in this thesis is physics 
without gravity on a flat space time. Cosmology, however, enters through 
the back door several times. 

The interactions of the SM are gauge interactions of the electroweak [4- 
6] and strong [7, 8] interactions (see Sec. 2.1.1), the self-interaction of 
the scalar sector leading to spontaneous symmetry breakdown [9, 10] in 
the electroweak interaction [11-14] (see Sec. 2.1.2) and the interactions 
of the SM scalar field with the matter fermions which give the flavor of 
the theory (Sec. 2.1.3). 


2.1.1 The Gauge Part 


Gallia est omnis divisa in partes tres. 
—Gaius lulius Caesar, De hello Gallico 

The gauge interactions of the SM are displayed by the group structure 
SU(3)c X SU(2 )l xU(l)y, 

where the strong and weak interactions are governed by the SU(3)(. and 
SU(2 )l factors, respectively. The remaining U(l) factor is of the weak 
hypercharge Y. Matter fields (fermions) are placed in either fundamental 
or trivial (i. e. singlet) representations of the gauge groups. In this way, we 
can set up the matter content of the SM: there are twelve quarks where 
six of them interact with all three gauge interactions and further six are 
accomplished to fill the Dirac spinor as they are the missing right-handed 
fields. Eurthermore two leptons interact under SU(2 )l x U(l)y and one 
lepton (right-handed electron) as SU(2 )l singlet carries only hyper charge. 
There is no right-handed neutrino needed in the SM because it would be 
a complete singlet under the full gauge group. We have 

Ug UA URr URg llRb 

^R,r ^R,g ^R,fe 




All Gaul is divided 
into three parts. 


The strong 
interaction couples 
only to colored 
particles, quarks 
and gluons, where 
the SU(2)l only 
interacts with 
left-handed 
fermions. 
Right-handed 
fermions are 
singlets under 
SU(2)l. 
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For Dirac spinors, 

\p — \p ^ _ 


The T° are the 
Gell-Mann matrices 
and T; the Pauli 
matrices. 


We thank Luca Di 
Luzio for sharing 
this observation 
with us. 


We have put the 
contributing 
authors in 
alphabetical order. 
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where r, g, b label the three color degrees of freedom (red, green, blue) 
and the left and right projections of Dirac fermions are given via 

= withP l = ^( 1 - 75 ) and PR = ^(l + r 5 ) 

for a generic Dirac spinor xp. 

For some reason, the fermion content of the SM is triplicated and each 
copy of the particle set (generation) differs by the particle masses. The 
interaction with the Higgs boson, which sets the masses, also generates 
transitions among the three generations. With only gauge interactions, 
we have for one SM generation of matter fields the kinetic Lagrangian 

£kin 

( S'^-i—r°G“’'"-i—f • W'^-i—TB'" 

2 2 2 

where the generators of SU(3)(, are denoted by r“/2 and those of SU(2 )l 

by Tj/2. The gauge couplings are labeled in an obvious manner and Y is 

the hypercharge operator. The spinor shall be used for a generic gauge 

multiplet; and are the corresponding SU(3)c, SU(2)l and 

U(l)y gauge vector bosons. 

Further generations can be just added in parallel because gauge inter¬ 
actions do not make inter-generation transitions, so 'Fpi'F —> 'Fjpi'Fj and 
i = l,...,nQ counts the number of of generations. An interesting obser¬ 
vation reveals itself in the “gaugeless limit” with gi_ 2,3 0 - neglecting 

fermion masses but obeying the gauge structure, the SM fermions have 
an enhanced [U(3)]^ symmetry (see also Chapter 6 ). The counting is sim¬ 
ple: one U(3) factor for each gauge representation because of three gen¬ 
erations. There are five gauge representations: left-handed quarks (color 
triplets and weak doublets), right-handed quarks (color triplets, twice for 
up and down sector), left-handed leptons (weak doublets) and the right- 
handed electron. In the gaugeless limit, the global symmetry gets even 
more enhanced: 

P gl,2,3“*0 

[U(3)]" -> U(45). 

The 45 is 3 X 15 which is just the triplication of one generation. In the 
presence of right-handed neutrinos, we have U(48) since the three right- 
handed neutrinos are complete SM singlets and simply added to the game. 

The gauge symmetries of the SM neither allow for gauge boson nor for 
fermion masses. Masses of gauge bosons violate gauge invariance and also 
Dirac mass terms of fermions are forbidden in the SM because left- and 
right-handed fields live in different representations of the gauge groups, 
so there is no gauge invariant way to construct them. We do not want 
to slaughter the sacred cow of gauge symmetry in order to introduce 
particle masses per brute force. A gauge invariant construction of par¬ 
ticle masses related to spontaneous symmetry breakdown of the gauge 
symmetry can be achieved via the Brout-Englert-Guralnik-Hagen-Higgs- 
Kibble [11-14] mechanism briefly described and introduced in the follow¬ 
ing section. 
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2.1.2 The Higgs Part 


I think, we have it. 

—Rolf Heuer 

The Higgs boson was the missing piece in the SM postulated 1964 and 
finally discovered in the first run at the LHC [15, 16]. It is a necessary 
ingredient to perform the spontaneous breakdown of electroweak gauge 
symmetry in the SM. Spontaneous symmetry breaking happens, once the 
ground state of the theory does not respect the initial symmetry anymore. 
In the SM we have an SU(2) x U(l) gauge symmetry, which is sponta¬ 
neously broken to the electromagnetic U(l). Unfortunately, the theory 
does not break itself—up to now, we have the gauge fields and fermions 
transforming under the fundamental representations. To break the gauge 
symmetry, we have to introduce additional fields. The most economic ex¬ 
tension of the field content is one further fundamental representation 
of a scalar H = (/t+,/t°), which is an SU(2) doublet. We allow for self¬ 
interaction of the scalar field and write down the following potential 

y (H) = + A (h'^h) , (2.3) 

where gauge invariance (no linear and cubic terms) and renormalizability 
(no monomials higher than four) dictate this structure. The sign of the 
/i^-term now decides whether the symmetry is broken or not, whereas the 
A-term has to be chosen positive for V to be bounded from below (for an 
extended discussion about the stability of scalar potentials see Chapter 4) . 
The minimization with respect to the neutral component (pi^ > 0) gives a 
vacuum expectation value (vev) v to the scalar doublet field 



We have a freedom of SU(2)-rotation (which is a gauge choice), therefore 
we can choose the vev to be in the neutral component since electric charge 
is not supposed to be broken. The unbroken generator is given by 

^(t3 + T)(H) = 0 

and the electric charge as combination of weak hypercharge Y and the 
third component of weak isospin T^w is given via the Gell-Mann-Nishijima 
relation 

Q = T3w + (2.4) 

Expanding around the minimum, we obtain the physical Higgs boson 
and the unphysical charged and pseudoscalar Goldstone bosons and 
X^, respectively: 

A(v + v« + i;r») 




Eureka! 

—Archimedes 


Gauge fields 
transform in the 
adjoint 

representation. 

The upper 
component carries 
electric charge +1 
where the lower 
one is electrically 
neutral, for this 
purpose the 
\J(1 )y- charge of 
that scalar has to be 
-|-1 according to 
Eq. (2.4). 


This way, we 
recover the 
unpleasant factor 
1 / from above 

which we keep to 
coincide with half 
of the literature. 
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Doing this 
procedure, the 
Goldstone bosons 
are “eaten” by the 
gauge bosons. 


The main point is, 
that the gauge 
boson masses can 
be calculated as 
combination of 
gauge couplings 
and vevs, which is 
true in any 
spontaneously 
broken gauge 
symmetry. 


Yukawa introduced 
a fermion-scalar 
interaction to 
describe 
pion-nucleon 
interaction in the 
same way. This 
coupling is still 
called Yukawa 
coupling though 
used in a different 
sense. 
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The Goldstone bosons remain massless and can be absorbed by gauge 
choice into the gauge bosons which acquire masses and a longitudinal 
degree of freedom. Through the kinetic couplings to the gauge bosons, 
ip^H), and the existence of v, they acquire partially a mass and 
mix to the physical W bosons (electrically charged) and the Z which is 
neutral and a mixture of and B^. The photon (the orthogonal state 
to the Z) remains massless. Details are omitted because this is common 
textbook knowledge, see e. g. [17, 18]. The masses are 

A / r 2 2^2 j aa -2 2 ^ 1~^^2 

= ^ ^ and — - —. 

The ratio of the two gauge couplings determines the weak mixing angle 
6^, which is the angle of the SO(2) rotation transforming into 

tanS^^, = gifg 2 - Inverting the mass relations, we obtain with 
the measured masses and gauge couplings v = 246 GeV which sets the 
electroweak scale. 

2.1.3 The Flavor Part 

The bodies of which the world is composed are solids, and therefore have 
three dimensions. Now, three is the most perfect number,—it is the first of 
numbers, for of one we do not speak as a number, of two we say both, but 
three is the first number of which we say all. Moreover, It has a beginning, a 

middle, and an end. 

—Aristotle 

The Higgs scalar is there, spontaneous symmetry breaking has happened, 
especially we have /t° = (v + (p°)/ V2. We exploit this fact to generate 
fermion masses via the same mechanism, coupling the Higgs doublet to 
the fermions via Yukawa interactions. To construct gauge invariant Dirac 
mass terms, we have to contract the SU(2) doublets with the Higgs dou¬ 
blet and append the singlet right-handed fermions: 

-4uk = • Hdnj - • Hunj + + h- c- , (2.5) 

where the dot product denotes SU(2) invariant multiplication (which gives 
the minus sign for up-type Yukawa) and the charge conjugated Higgs dou¬ 
blet 

= ('_?)■ 

The couplings of Eq. (2.5) with / = u, d, e are arbitrary matrices in 
flavor (i. e. generation) space where the indices i,j = 1, ...,ng count 
the number of generations. Up to now, we know Rq = 3 copies of SM 
fermions, where a fourth sequential fermion generation is excluded after 
the Higgs discovery [19]. 
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Eqs. (2.2) and (2.5) define two different bases which can be trans¬ 
formed into each other and produce fermion mixing phenomena. After 
the Higgs doublet acquires its vev, Eq. (2.5) gives masses to the fermions, 
where 

mf = ( 2 . 6 ) 


defines the mass matrices. The basis, in which is diagonal is therefore 
called mass basis. In contrast, the gauge interactions of Eq. (2.2) define the 
interaction basis. The interplay of gauge multiplets combining two a priori 
independent mass matrices (for the up and the down sector) results in the 
generation (flavor) mixing of the weak charged interaction. To elaborate 
on this feature, which prepares especially for Chapters 3 and 6, we rotate 
into the mass eigenbasis using bi-unitary transformations 

s{y^ ' =Y^ = diagonal. (2.7) 

In view of the gauge representation of the fermions as indicated in Eq. (2.5), 
we cannot rotate left-handed up and dovra fermions independently since 
both form a doublet. Obeying the gauge structure, we have three indepen¬ 
dent rotations in the quark sector: 


It can be easily 
shown that the left 
and right unitary 
matrices Sl r 
diagonalize the left 
and right Hermitian 
products, 

S^YY'^ {Sj[and 
SrY'Y (Sr)'^ with 

«r)'Sl,r = 1 . 


Qtd 


(2.8a) 


“r,; “ '^Rdj^Rj’ 

(2.8b) 

dK,i 

^R,i~ ■^r.o'^R,;- 

(2.8c) 


This set of transformations is not sufficient to simultaneously diagonalize 
the mass matrices m“ and m'^. We can determine and S“ from the 
up-type Yukawa coupling via Y“ = S^Y“ (S“)'; then we need a further 
unitary matrix V from the left to diagonalize 


yd ^ ^Qyd fd ^ yyd ^ ysQyd ^^dy 


The matrix V measures the misalignment of Yukawa couplings; if V = 1 
the alignment of up and dovra Yukawa would be exact and both up and 
down mass matrices are simultaneously diagonal. Now, we see the out¬ 
come of the transformation into the mass eigenbasis: the neutral current 
interactions (or SU(2 )l singlet-like) are unaffected and still flavor con¬ 
serving thanks to the unitarity of mixing matrices. In contrast, the charged 
current interaction Cqq = + h. c. reveals the mixing matrix V 

as leftover in the left-handed charged fermion current J^. Performing the 
transformations into the mass basis from above, we have 


Jr = “tr'^dL ^ (S?) ' (2.9) 

We identify in Eq. (2.9) the Cabibbo-Kobayashi-Maskawa [20, 21] (CKM) 
matrix as Vckm = 1^'- The CKM matrix describes the threefold mixing of 


8 PARTICLE PHYSICS IN A NUTSHELL 


CP is the combined 
charge-parity 
transformation. 

Parity 
transformations 
describe discrete 
transitions from 
left- to right-handed 
coordinate systems 
(and vice versa] via 
X —X. Charge 
transformations flip 
all charges similar 
to complex 
conjugation which 
flips the sign in 
front of the 
imaginary unit. 


The magnitudes of 
CKM elements can 
be displayed as 
follows \Vckm\ = 

f ' 


V ■ 



the SM generations and gives a possibility for CP violation, tvhich was the 
reason why Kobayashi and Maskawa extended the two-generation descrip¬ 
tion to a third generation. A mixing matrix of two flavors cannot violate 
CP because all complex phases can be absorbed in redefinitions of the 
fermion fields whereas a 3 x 3 unitary matrix has three angles and six 
phases from which only five phases can be removed because one global 
phase can stay arbitrary. The most convenient way of parametrizing the 
three rotations with one complex phase was introduced by [22] and is 
commonly used as “standard parametrization” [23]. This parametrization 
decomposes the CKM matrix into three successive rotation with one mix¬ 
ing angle for each rotation in the 2-3, 1-3 and 1-2 plane, respectively: 


^CKM — ^23(^23)^13(^13) 5 ckm)^ 12 (^ 12 ) 


r 




100 

0 ^23 S 23 

“•523 ^23 

r 

^12^13 


r 


^13 

0 


0 Si 3 e 
1 0 


i 5 c 


.A 




Cl3 

• 512^13 


r 


V 


^12 

“•512 

0 


“•512^23 “ C12'5235i36'^'^'™ ^ 12^23 “ •5l2'523'5l36‘^‘^’™ 

512523 “ Ci2C23Si3e'^CKM —C12523 “ 5i2C235i3e' 


( 2 . 10 ) 

5i2 0 

^12 0 
0 1 ) 

^^g-iScKM^ 

523^13 ’ 

^23^13 y 


with Cjj = COS 6ij, Sij = sin0;j and is the CKM CP-phase. We single 
out two important features of this parametrization which will be conve¬ 
nient in the further course of this thesis: (a) the separation into three 
rotations in three different flavor planes allows to keep track of the indi¬ 
vidual contributions in the final result (this can be seen from the upper left 
matrix elements, where ~ and (b) the CP phase sits in the 1-3 
rotation which for both quark and lepton mixing has the smallest angle. 

The SM as described so far has no room for lepton mixing. The Yukawa 
Lagrangian (2.5) can be exactly diagonalized for the charged leptons, be¬ 
cause we have two free rotations that can be absorbed into redefinitions 
of the lepton fields. Mass terms for neutrinos are not scheduled in the SM. 
As it is a minimal theory, there are no right-handed neutrinos since they 
are pure gauge singlets and do not interact. The only interaction they 
would have are Yukawa interactions with left-handed neutrinos. Flavor 
mixing, however, needs fermion masses. So the observation of neutrino 
oscillations (see as reviews [23, 24, and references therein]) already hints 
towards new physics beyond the minimal SM. More about neutrino flavor 
follows in Sec. 2.4. 

The CKM matrix has been measured with amazing precision [23] 


I^ckmI — 


^0.97427 ±0.00014 0.22536 ± 0.00061 
0.22522 ±0.00061 0.97343 ± 0.00015 


V 


0.00886 


+0.00033 

-0.00032 


0.0405 


+0.0011 

- 0.0012 


0.00355 ±0.00015^ 
0.0414±0.0012 . 

0.99914± 0.00005^ 

( 2 . 11 ) 
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De gustibus non est disputandum. 
—Jean Anthelme Brillat-Savarin 

Supersymmetry (SUSY) is the only symmetry extension of the S-matrix 
which is not a direct product of any internal symmetry group and the 
Poincare group of space-time as stated in the Coleman-Mandula (CM) 
theorem [26]. As loophole in the CM theorem, the Haag-Lopuszaiiski- 
Sohnius theorem [27] proposes supersymmetries as extension of the space- 
time symmetry (Poincare symmetry) in a way that fermionic generators 
(in the spinor representation of the Lorentz group) transform bosons into 
fermions and vice versa. 

The fermionic generators of SUSY obey a so-called pseudo Lie alge¬ 
bra [27], which is the anti-commutator relation 

{Q^Q“} = 2r^^P^5^“, (2.12) 

with the Dirac y-matrices as structure constants (together with a Kronecker- 
5). The indices N,M count the number of SUSY generators. iV = 1 cor¬ 
responds to one generator as in the MSSM. The operators are Ma- 
jorana spinors, where a is a spinor index; the Hermitian conjugate is 
j ', and the generator of space-time translations also known 

as 4-momentum vector. Conserved currents related to SUSY (“supercur- 
rents”) are spin 3/2 currents, where the conserved quantity of the energy- 
momentum is the energy-momentum tensor, a spin 2 quantity, see 
e. g. [28]. In this way, SUSY is a candidate to combine gravity and gauge 
theories—the field connected to the supercurrent is the spin 3/2 gravitino 
whereas the one related to the conserved energy-momentum is the gravi¬ 
ton. However, SUSY extends the Poincare algebra of space-time and there¬ 
fore also the structure of space-time itself has to be extended, which leads 
to the “superspace”. The Poincare group includes the four-dimensional ro¬ 
tations of the Lorentz group SO(l, 3) and translations in Minkowski space 

= + a\ 

with the Lorentz transformation and a constant vector a^. Generator 
of spatial translations is the 4-momentum P^, for which 

[P^„P^]=0, (2.13a) 

Pp] ~ ^ {SfjLv^P ~ Snp^v) (2.13b) 

hold with L^'' = i {x^d'^ —x^'d^) = x^P'^ — x'^P^ and 

= i (^g^Pincr _ gPPiva _ gva^pp gP^gvp'^^ ( 2 . 13 c) 

Eqs. (2.13) are called the Poincare algebra. 

The fundamental fields of the SM fit into irreducible representations of 


A direct product 
means that any 
internal symmetry 
generator shall 
commute with the 
generators of 
Poincare symmetry 
Fermionic 
generators, 
however, obey 
anti-commutation 
relations which lead 
to so-called graded 
Lie algebras (details 
in any good 
textbook about 
supersymmetry, 
e.g. [25]). 


The superspace is a 
conceptually 
different concept 
than the object 
introduced in [29]. 
We consider 
Minkowski space 
with the metric 
Sjrv 


Scalars (trivial 
representation), left 
and right chiral 
spinors, vectors, ... 
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can be defined 
via the y-matrices, 
— 4 [T/jjTv] ■ 

Left and right chiral 
superfields, Vector 
superfields, ... 

The spinorial 
coordinates are 
Grassmann 
numbers, 6^ = 0. 


The fermion field ^ 
is a Weyl spinor. 


The complex 
conjugated field 
is called right-chiral. 


The bar over is 
not to be confused 
with the Dirac-bar. 
It shall keep in 
mind that the 
component fields 
are conjugated. 


the Poincare group whose invariants are related to mass and spin. Extend¬ 
ing the Poincare algebra with the fermionic generators from Eq. (2.12), 
one gets in addition 

[P^Qa]=0, (2.14a) 

= (2.14b) 

Eqs. (2.12), (2.13) and (2.14) form the super-Poincare algebra. Particles 
of supersymmetric field theories fit into irreducible representations of the 
super-Poincare algebra. A supermultiplet contains bosonic and fermionic 
degrees of freedom; in a similar manner fermionic coordinates 6 are 
needed. Superspace coordinates are complex, 

yPi = -iOa^O, 

(2.15) 

yn = xf^ + ieaf^e. 

We have = (l 2 ,cr‘) (i = 1,2,3), the vector of Pauli matrices. All su¬ 
perfields fF are functions of the superspace coordinates, fF{x,9,9). 


CHIRAL SUPERFIELDS The lowest representation of the Super Poincare 
algebra are chiral superfields, that contain a scalar field cf) and a fermionic 
component Additionally, there is an auxiliary field F that can be elim¬ 
inated with the equations of motion (eom) because there are no kinetic 
terms for F in the SUSY Lagrangian. These eom result in scalar mass terms 
(F-terms). Eor completeness, we give the full superspace expansion of the 
left chiral superfield ^ = {(j),S,,F} and its complex conjugate 

4>{y, 9) = (p(y)-\- s/29^(y) + 99F(y), 

— I - — — (^ 2 . 16 ) 

^'fy, 9) = (j>*(y) + V^9f(y) + 99F*(y). 

It is convenient (and also convention) to work only with left-chiral su¬ 
perfields, so right-handed fermions of the SM are squeezed into the left- 
chiral representation via charge and complex conjugation. If we have a 
SM fermion /l and its scalar superpartner /l, they fit into the left-chiral 
Fl = {/l,/l}- Their right-handed colleagues are put into a left-chiral su¬ 
perfield as Fr = {f^,ff}. It is necessary to treat right-handed fermions 
separately, because they transform differently under the SM gauge group. 
The gauge representation and the Poincare representation, however, must 
not be mixed up. Poincare left-handed fields can be obtained via charge 
conjugation. Note that charge conjugation does not change the gauge rep¬ 
resentation from the singlet to a doublet representation. 


VECTOR SUPERFIELDS Chiral superfields are spin 0 and spin 1/2 fields. 
The spin 1 gauge bosons of the SM have to have a different super-Poincare 
representation. Vector superfields V are real fields, so = V, and can 
be constructed out of chiral superfields, see e. g. [30]. There is a gauge 
freedom (“supergauge”) in the space of vector superfields which allows to 
choose a particular gauge to reduce the most general representation of a 
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vector superfield to one vector field A^(x), one complex two-component 
spinor A(x) and one auxiliary field D{x) which again can be eliminated 
using the eom. This special supergauge choice is known as Wess-Zumino 
gauge [31] and we have 

^w-z(^> 0, 6) = 6a^6A^{x) + 966l{x) + 6X{x)66 + -0699D(x). 

(2.17) 

INTERACTING SUPERFIELDS The non-gauge interactions of chiral su¬ 
perfields can be written in the following superpotential: 


W(^) = hi<Pi + 


(2.18) 


The superpotential is a holomorphic function of chiral superfields, con¬ 
tains therefore only left-chiral (or only right-chiral) superfields and has 
mass dimension three. The supersymmetric Lagrangian can be obtained 
from the superpotential as the “highest component”, i. e. the coefficient in 
front of 00. This can be seen from the definition of the action [32] 


+ W (^) 5 (0) + W''' (^'^) 5 (0) 





s = 

d^x 

d^ 0 d^ 0 



_ 


The kinetic term can be put into (super) gauge invariant shape by 
inserting the gauge supermultiplet 

with some gauge coupling g, such that 




+ 


eeee 


W(^) 


-|-h. c. 


ee 


(2.19) 


The trilinear 
couplings of the 
superpotential are 
dimensionless and 
symmetric in 

the bilinear 
couplings mij have 
mass dimension one 
and the tadpole 
coupling hi 
dimension two. 

Integration over 
Grassmann 
numbers behaves 
like differentiation, 

“J d^e = d^/de^”. 


\x 


The notation | 
means that the 
coefficient in front 
ofX is taken. 


SUPERSYMMETRIC MASS TERMS AND SCALAR POTENTIALS There are 
Still auxiliary fields around. Eliminating the f-fields with dC/dF* = 0 
and dmdF; = 0 leads to 


dW'^ 


d^'l 
* dW 

77^ _ 


= -h* 




— -hi - niijclij - -fijk4>j4>k> 


which results in the scalar F -term potential 


Vp{ct>,cl>*)=F*Fi = 


dW'^ 


d4>l 


dW 


d^i 


( 2 . 20 ) 


( 2 . 21 ) 


The D-terms are eliminated analogously via dLIdD = 0 with 

= -gclijT^ctij, (2.22) 

for each gauge symmetry with coupling g and generators 7“.. Analogously, 
we get the D-term potential 

Voitlt, 4>*) = = g" . 


The notation | 
means that the 
derivative is 
evaluated at 
6=0 = 6 . 


(2.23) 
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Superpartners are 
abbreviated with a 
tilde over the 
symbol, f for a 
sfermion or squark 
q; similarly W, B; 

admixtures like 
neutralinos or 
charginos . 


The term soft 
breaking shall 
reflect the fact that 
all SUSY breaking 
couplings are 
related to the 
couplings of the 
superpotential and 
the full theory still 
is supersymmetric, 
only the ground 
state breaks the 
symmetry. Soft 
breaking terms are 
hence related to a 
vev and are 
dimensionful 
quantities that do 
not introduce 
quadratic 
divergences. 


In a SUSY theory, 
bosonic and 
fermionic 
zero-point energies 
exactly cancel to 
zero in contrast to a 
non- 

supersymmetric 

theory. 


SETTING THE LANGUAGE Fermions of the SM get scalar superpartners 
in a supersymmetric theory. Those are called scalar fermions or sfermions. 
vector fields of the SM get fermionic (Majorana) spinor partners that are 
denoted with the suffix -ino like gaugino, gluino, electroweakino. The su¬ 
perpartners of the Higgs scalars are fermions as well and therefore also 
called higgsinos (though Higgs fields are chiral superfields). 

2.2.1 How to Break Supersymmetry 

Unfortunately, SUSY has not yet been observed in fundamental interac¬ 
tions. If so, e. g. charged scalar particles with the mass of the electron 
must have been seen. Since no selectrons appear in atomic physics, no 
squarks have been detected in high energy collisions, gluinos and elec- 
troweakinos hide maybe somewhere, SUSY has to be badly broken. Actu¬ 
ally, SUSY breaking is constructed in a way to happen “softly”. However, 
mass terms for superpartners are needed to shift their masses into the TeV 
regime to cope with their hide-and-seek play. We do not go into the de¬ 
tails of collider phenomenology, maybe there are some stripes left in the 
SUSY landscape to find at least some superpartners at the electroweak 
scale (100GeV rather than 10TeV). The SUSY corrections we calculate 
and exploit in Chapter 3 anyway are non-decoupling contributions. So, if 
all SUSY parameters are shifted uniformly to higher scales, the results do 
not alter. In this way, we may use flavor physics as an indirect probe of 
SUSY breaking. Soft breaking is expected to be the result of a spontaneous 
symmetry breakdown and all mass terms and mass dimensional couplings 
are related to some vev. In this case, if SUSY is broken via a process as the 
Higgs mechanism, all masses generated by this breaking are of the same 
scale. However, breaking of SUSY is different from spontaneous breaking 
of any internal symmetry: broken SUSY leads to a non-zero vacuum en¬ 
ergy density since a supersymmetric ground state always has exactly zero 
energy. We are only interested in the phenomenological output of SUSY 
breaking that can be described very elegantly as shown below. There is 
a vast amount of concepts on the market which cannot be reviewed as 
it is to be seen complementary to most SUSY phenomenology. We also 
can only refer to a small subset of literature which comprises interesting 
ideas like dynamical SUSY breaking [33-37]. SUSY breaking occurs in a 
“hidden” sector where the highest component of a superfield acquires a 
vev—in case of chiral superfields one has f-type breaking [38]; in case of 
vector supermultiplets D-type breaking [39]. Neither in the description of 
O’Raifeartaigh nor Fayet-Iliopoulos a deeper reason for the vev is given 
as it is in the dynamical models. SUSY breaking has then to be trans¬ 
mitted from the hidden to the visible sector via some mediator fields— 
popular attempts are gauge mediation, see e. g. as review [40] (also in 
combination with dynamical breaking [41, 42], “supercolor” in contrast 
to technicolor [43-45]), and anomaly mediation [46, 47]. Minimal super¬ 
gravity [48] allows to combine and break local iV = 1 SUSY and grand 


2.2 SUPERSYMMETRIC EXTENSIONS 13 


unified theories. Finally, the determination of the Higgs boson mass allows 
to slightly discriminate between the different types of models [49] . 

For the phenomenological processing of soft SUSY breaking, we may be 
ignorant of the dynamics behind SUSY breaking and mediation of SUSY 
breaking. Instead, the soft breaking terms are added to the supersymmet¬ 
ric Lagrangian without deeper knowledge of their origin. 


^ — -^susY + '^soff (2.24) 

The soft breaking Lagrangian comprises mass terms for the scalar 
components of chiral superfields and gaugino Majorana mass terms 
for the fermionic parts of vector supermultiplets. Moreover, mimicking 
the polynomial structure of the superpotential there are trilinear, bilinear 
and linear terms in the scalar components of chiral superfields allowed 

-Csoft = - + h. c. ) 

(2.25) 

Si;^i<^j+Ci</)l+h.C. J . 

Certainly the terms of (2.25) are not allowed to break internal symmetries. 
The A- and S-terms are symmetric in their indices and obviously carry 
mass dimension one and two, respectively. C-terms are only present if 
there are tadpole terms in the superpotential which only may occur for 
gauge singlets. We see no connection to superpotential parameters in the 
SUSY breaking terms and stick to the notation of Eq. (2.25) instead of 
factorizing artificially the superpotential parameters as partially done in 
the literature [32] writing e. = Aij^/ 


, 1 1 


2.2.2 The Minimal Supersymmetric Standard Model 


We now have the ingredients to set up the Minimal Supersymmetric Stan¬ 
dard Model (MSSM)—^which indeed comprises broken SUSY. The mat¬ 
ter content of the MSSM is given by the chiral superfields of 3 x 15 SM 
fermions (2.1), the vector supermultiplets of the SU(3)c x SU(2 )l x U(l)y 
gauge interactions and the Higgs—which is part of a chiral superfield and 
has to be doubled [52]. The today’s language of the MSSM was basically 
set by [53]; a comprehensive overview of supersymmetry, supergravity 
and particle physics was given in [54]. 

The superpotential of the MSSM is given by 


W, 


MSSM 


- - Y,)Ha ■ • Ql,; 

(2.26) 


where the bilinear ju-term is the only dimensionful parameter of the super¬ 
potential and itself no SUSY parameter which causes conceptual problems 
and solutions to that problem [55, 56] . In order to have Yukawa couplings 
to both up and down type fermions, there are two Higgs doublets and 


Soft breaking does 
not induce 
quadratic 
divergences at 
one-loop [50], 


The signs of A-, B- 
and C-terms are of 
no meaning and 
depend on the 
convention. To 
interpret soft 
breaking masses as 
masses, their signs 
are fixed. 

We do not consider 
the 

“non-holomorphic” 
A-terms like 

[51]. 


Looking at the dates 
of the most 
important SUSY 
publications, we 
find the golden age 
of SUSY about more 
than thirty years 
ago. 

Fermion masses are 
with (h°) = Vu and 
{h°) = Vd given by 
m“ = VyJ'^lAi, 
m^ = vyY^lA2 
and 

m® = vyY^lAi. 


has hypercharge 
+ 1/2, Ha - 1 / 2 . 
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Hu couples to up 
type Helds, to 
down type fields in 
Eq. (2.26). 

The doublet 
superfields are 
correspondingly 
Ql = {qL.^L} and 
fh = 

Despite of the 
subscript j;, are 
left-handed Weyl 
fermions. 

Csoft together with 
F- and D-terms 
gives the mass 
squared matrices of 
sfermions and the 
ga ugino/higgsino 
mass matrices. 
Diagonaiization of 
these matrices 
result in flavor 
changing vertices. 
Mass and 
diagonaiization 
matrices are 
specified in App. A. 

Clarifications about 
the spinor notation 
in App. A. 

The F-terms give 
additional 
interactions of 
sfermions and 
Higgses and are of 
importance for 
sfermion mass 
terms and the 
analysis of the 
minimum structure 
of the scalar 
potential. The 
D-terms are 
determined by 
gauge couplings 
squared and give 
the quadrilinear 
terms in the Higgs 
(and sfermion) 
potential. 
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with different L7(l)y-charges (capitals denote chiral superfields), 



(2.27) 


The left-handed quarks and leptons form the SU(2)L-doublet chiral super- 
fields Ql = and Ll = with L7 l = {u^,uj, Dl = {4.41 

up and down (s)quarks and Ml = {vl, Vl}, £l = (s)neutrino and 

(s)electron, respectively. The SU(2)l singlets are in the left-chiral repre¬ 
sentations [7 jl = {u*,u^}, Djl = {d*,d^} and = {e*,e^}. Generation 
indices are suppressed, where in Eq. (2.26) i,j = 1,2,3. 


SOFT BREAKING IN THE MSSM SUSY has to be softly broken in the 
MSSM, so we set the soft breaking Lagrangian according to Eq. (2.25) 
with the fields of the MSSM and have 


^MSSM 

■"soft 


^ (*^ 0 ) ij + “R.i ('^'9 ij %d + ^R,i 4,; 

+ ^L,i + eR,i ('^e)i/^Rd 


+ 






iA+h.c. 


+ m^j/idl^ + m 2 jhu |2 + (s^hd-^u + h. c. ) 

+ 2 + h. c. j + — ^M 2 AA + h. c. ^ 

+ ^(M 3 l“A“+h.c.). 

(2.28) 


The scalar mass and trilinear terms are self-explanatory; is the Higgs 
B-term and gaugino masses are Ml 2,3 with labels according to the gauge 
couplings gi, 2 , 3 - The gauginos are written as Weyl spinors. 


THE 2hdm OF THE MSSM SUSY dictates the Lagrangian: the supersym¬ 
metric part is related to the superpotential which sets the interactions 
among chiral superfields; the SUSY breaking part is given by The 

scalar potential, however, does not only include Usor = —Usoft> but also 
f-terms and D-terms. We have two Higgs doublets, so the most general 
Higgs potential resembles the potential of a two-Higgs-doublet model 
(2HDM) [57, 58]: 

^ =^11 K^d + KK + (^12 K-hd + h. c. ) 

+ y ('■>d) U ^ U A, (/.;/.„) (/.>d 

r ^ r t ^ f ^5 r ^2 (2.29) 

+ A 4 (/t;/ta) (h ' HJ + ( y (Hu • Ha) 

- Ae (H+Hd) (h„ • Ha) - A 7 (h+H J (h„ • Ha) + h. c. ). 
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In the MSSM, the parameters of Eq. (2.29) are calculated by the methods 
described in this chapter above. Working this out, one finds that no con¬ 
tributions at the tree-level for the self-couplings A 5 5 7 exist. Moreover, the 
potential is constructed in such a way, that A 4 does not show up in neutral 
Higgs interactions. The mass terms are a combination of the p-parameter 
and soft breaking masses: 

= + + ^ 12 =^ 1 ^’ (2.30a) 

K 2 = -h = , ^4 = y • (2.30b) 

For a reasonable theory, the scalar potential has to be bounded from below, 
i. e. there are no directions with V —> —00. There are three simple condi¬ 
tions to be fulfilled in order to avoid unboundedness from below [58], 


> 0 , 


A 2 > 0 and 


A3 > —\/A3A2, 


(2.31) 


which are always fulfilled in the MSSM at tree-level with (2.30b). No more 
conditions are needed for the tree-level MSSM, because A 5 g 7 = 0. 

If the mass matrix formed out of mfj has one negative eigenvalue, the 
scalar components of the Higgs doublets acquire vevs 





assuming that electromagnetic U(l) stays intact. The individual vevs are 
fixed via the W mass, + v^ = v^ = 4M^ /g^, and the ratio is to be seen 
as free parameter 

tan^ = —. (2.32) 


The requirement of spontaneous symmetry breaking gives relations be¬ 
tween the tree-level parameters of the 2HDM potential (2.29) 


m?, = m ?2 tan /3-cos 2^ A^, 


mL = m?2Cot/3 H-cos2/3A2. 


(2.33a) 

(2.33b) 


We get the mass matrices for CP-even and CP-odd as well as charged com¬ 
ponents as second derivative of the potential with respect to the corre¬ 
sponding fields. Expanding the Higgs doublets around their vevs. 




h, = 


-Ad“ 


This can be seen 
writing = 

u—uO _ ^0*u— 

“u '^d '^u '^d ' 


We then have 
Vu = V sin p and 
Vd = V cos p. 


Conditions (2.33) 
ensure that the 
global minimum of 
(2.29) is 

determined by Vj, 
and Vd. 


we have eight dynamical fields (charged fields are complex) out of which 
three Goldstone bosons have to be eaten by the gauge fields; five physical 
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All we have at hand 
are g^, g| and vev 
relations, and 
tanp. 


More about 
effective potentials 
in Chapter 4. The 
idea is to calculate 
the Higgs potential 
at one- or two-loop 
order and obtain 
the masses as for 
the tree-level 
potential. 
We use FeynHiggs 
at some later point 
to determine the 
lightest MSSM 
Higgs mass. 


fields remain: two CP-even (/t° and H°), one CP-odd (A°) and the charged 
Higgses (H"*"). If CP is conserved and not spontaneously broken; otherwise 
/t°, and mix. The mass of the pseudoscalar A° can be related to the 
yet unconstrained tree-level mass parameter, 2 m ^2 = ttt^oSin2^, and is 
then also a free parameter of the theory. 

The MSSM predicts a rather light Higgs boson mj,o < if no radia¬ 
tive corrections are taken into account. Already one-loop corrections lift 
the lightest Higgs mass well above Mg [59-62] and are quite needed, if 
we want to explain the discovered Higgs boson mass at 125 GeV [15, 16]. 
The dominant radiative corrections at one-loop are related to the large 
top Yukawa coupling and originate in diagrams with stops or tops. If 
already one-loop corrections are large, two loops are of equal importance 
and have been calculated diagrammatically [63-66] as well as in the effec¬ 
tive potential approach including two-loop effects [67-72]. The indepen¬ 
dent approaches of diagrammatic and effective potential calculations were 
shown to coincide up to knovra differences [73]. Current up-to-date tools 
for numerical evaluation of MSSM calculations obtain those corrections 
(and some more) as FeynHiggs [66, 74-77] or popular MSSM spectrum 
generators as SoftSUSY [78], SuSpect [79] and SPheno [80]. The three- 
loop SUSY-QCD effects are also available [81, 82] and ready to use in the 
computer code H3m [82]. Three-loop corrections are not only important 
to reduce the theoretical uncertainty in the precise prediction of the light 
MSSM Higgs mass but also give important contributions for multi-TeV 
stops [83]. Spectrum generators may be combined and compared using 
the Mathematica package SLAM [84]. 

If superpartners are generically heavy and the Higgs scalars of the 2HDM, 
however, remain light, the MSSM can be matched at the full one-loop level 
to an effective 2HDM where the couplings are determined via SUSY pa¬ 
rameters at the decoupling scale [85]. 


2.2.3 Radiative Flavor Violation in the MSSM 


Evaluating running 
MS masses at the 
SUSY scale. 


The masses and mixings of the fermions in the SM (and MSSM) enter 
via the Yukawa couplings which are ad hoc parameters, though dimen¬ 
sionless. A problem or rather a puzzle in that respect is the question why 
the masses (Yukawa couplings) of the first two generations are so small 
compared to the third generation. Two-Higgs-Doublet models give a han¬ 
dle on the comparability of top and bottom mass via tan^, since in the 
MSSM mj/nib = tan^Yj/Y^ and mt(lTeV)/mb(lTeV) 60: assuming 
tan 13 = 60, the top and bottom Yukawa coupling are of equal size, Yj Y^,. 

It is intriguing to keep only the large third generation Yukawa couplings 
and postulate 


Y 


u,d,e 


^0 0 
0 0 
^0 0 


0 

J 
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Figure 1: Flavor changing self-energies in the MSSM: gluino-squark, neutralino/chargino-squark and 
neutralino-sneutrino/chargino-slepton loops (from left to right). Similar diagrams exist for 
the charged lepton propagator. 


via imposing flavor symmetries as Z 2 or U(2) for the first ttvo genera¬ 
tions, see e. g. [87-92]. The idea of vanishing zeroth order fermion masses 
with a mass generation at the loop-level was already pointed out by Wein¬ 
berg [93] and applied in the context of grand unified models [94, 95]; an 
exhaustive analysis of radiative fermion masses in grand unified theories 
can be found in [96]. Radiative SUSY mass models allow to produce cer¬ 
tain hierarchies in the mass matrices [97], induce chiral symmetry break¬ 
ing via soft SUSY breaking [98] and generate Yukawa couplings radia- 
tively [99]. Moreover, SUSY threshold corrections are important to obtain 
Yukawa unification [100]. Imposing non-minimal flavor violation (NMFV) 
in the MSSM, radiative flavor violation (RFV) can be used to suppress the 
SUSY flavor changing contributions and generate the quark mixing of the 
CKM matrix radiatively [87-90]. On the other hand, some contributions 
can be enhanced [102]. RFV in the MSSM has been extensively studied 
and constrains the parameter space giving additional relations between 
flavor observables [103-106]. 

We do not follow the very ambitious goal at the moment to simultane¬ 
ously generate the fermion masses radiatively and accomplish for the mix¬ 
ing. The motivation behind RFV is to see the flavor changing off-diagonals 
in the CKM matrix (2.11) as small perturbations arising from higher-order 
effects. Flavor mixing in this description enters via loops of supersymmet¬ 
ric particles. The general soft SUSY breaking Lagrangian of Eq. (2.28) 
has arbitrary flavor structure that can be confined using either minimal 
flavor violation (MFV) techniques [107-109] or RFV Key ingredient are 
flavor changing self-energies shown in Fig. 1 which can be decomposed in 
chirality-flipping and chirality conserving pieces 


Sfi{p) = Sfl{p^)P^ + S}f{p^)P, 


+p 


Ji \F 

'sjl{p^)p^ + sff{p^)p^ 


(2.34) 


with and 27 “(p2) = 

The application of RFV to the lepton mixing matrix is given in Chapter 3, 
where we calculate SUSY threshold corrections to degenerate neutrino 
masses and generate both mass splittings and mixing angles. In the further 
course of Chapter 3, we apply the mixing matrix renormalization of [110] 


An amusing 
application of the 
radiative mass 
mechanism was 
found in the 
observation that 

rrie/m^ 

0{a) [86]. 

A very brief 
overview about 
Grand Unification is 
given in Sec. 2.3. 
With the term 
NMFV we denote 
potentially arbitrary 
flavor structures in 
the soft breaking 
terms, especially 
the A-terms. The 
consequences of 
NMFV in Sugra 
theories were 
discussed in [101]. 


MFV means that all 
FV stems from the 
standard Yukawa 
couplings—in softly 
broken SUSY this 
means that 
e. g. trilinear 
A-terms are chosen 
aligned with 
Yukawa couplings. 

Seee.g. [110, 111]. 
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to the lepton mixing matrix via SUSY self-energies. This approach has 
already been used in the CKM renormalization of the MSSM [89, 90] 
and for the leptonic case [112], where the neutrino self-energies were 
omitted. They, however, may give sizable contributions if neutrinos are 
not hierarchical in their masses. 


2.3 EXTENDING THE GAUGE SECTOR 


Grand Uniffed 
theories also give 
natural 
explanations of 
puzzles like 
neutrino mass and 
charge quantization. 


Supersymmetry is a symmetry extension of the fundamental description of 
particle interactions. The symmetry group which is afflicted is the Poincare 
group of space-time. According to the Haag-Lopuszahski-Sohnius theo¬ 
rem, SUSY is the only non-trivial extension of the symmetries of the S- 
matrix. However, a trivial symmetry extension in terms of the Coleman- 
Mandula theorem is an enlargement of the internal symmetry group. This 
can be either accomplished by adding additional gauge group factors to 
the SM group or by embedding the SM gauge group into one symmetry 
group. The latter is known as Grand Unification (GU). Towards a higher 
symmetry description of nature, the most natural would be to combine 
SUSY and GU. We comment briefly on Grand and Partial Unification in 
the following to get a perspective on the unified picture. 


Grand Unification is 
basically driven by 
the fact that the 
three gauge 
couplings tend to 
unify—in the 
MSSM nearly 
perfectly at a scale 
QGUT = 2xl0i®GeV. 


We cannot even 
discuss the most 
generic aspects of 
Grand Unified 
Theories in any 
detail because we 
neither really need 
it in the course of 
this thesis nor do 
we perform any 
calculations in GUT. 
The D-factor is a 
kind of parity 
symmetry [116]. 

The labels s and a 
denote symmetric 
and antisymmetric 
represen ta tions. 


2.3.1 Grand Unification 


When shall we three meet again? 

—First Witch [William Shakespeare, Macbeth] 

The gauge groups and the gauge representations of the SM are chosen 
on purely phenomenological grounds. When the SM was proposed, there 
were also no hints for neutrino masses. In this respect, the particle content 
of the SM is minimal on the one hand and full of assumptions on the other 
hand. It is even more surprising that it perfectly fits into Grand Unification. 

For the purpose followed in this short section, we give the most intrigu¬ 
ing example of a Grand Unified Theory (GUT) which is astonishingly sim¬ 
ple, combines all different representations of fermions in the SM into one 
representation and gives an explanation for small neutrino masses. The 
smallest simple Lie group which includes the SM gauge group and pro¬ 
vides a single representation for the 15 + 1 SM fermions of one generation 
is SO(10), proposed by Georgi [113] and Fritzsch and Minkowski [114]. 
Moreover, SO(IO) can be decomposed either to SU(5) x U(l), which con¬ 
tains the Georgi-Glashow model [115], or 

SO(IO) ^ SO(6) X SO(4) SU(4) x SU(2 )l x SU(2)r x D, 

which hints towards partial unification and shall be explained in Sec. 2.3.2. 

With the SM fermions in the spinor representation of SO (10), there 
are three candidates for mass terms as 16 <E) 16 = 10^ © 120^ © 126^ to 
construct SO(10) invariant Yukawa couplings [117, 118]. Gonstraining 
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ourselves to symmetric representations (and therewith symmetric Yukawa 
couplings), the Yukawa Lagrangian is given by 

_^so(io) ^ ^ 126h, (2.35) 

where i,j = 1,2,3 count generations and H labels Higgs representations 
of scalars obtaining a vev. The coupling to the 126^ generates neutrino 
Majorana masses—for the right-handed neutrinos at a high scale and for 
the left-handed neutrinos via a seesaw type I-I-II combination at the elec- 
troweak scale [120]. 

The are a lot of issues to be addressed in SO(IO) GUT (also in SU(5) 
similar tasks arise). First of all, the large symmetry group has to be bro¬ 
ken to the SM group. To do so, there are in general several Higgs rep¬ 
resentations at work to perform the symmetry breaking steps. We want 
to live in a SUSY environment and therefore have to address one special 
issue of SUSY GUTs: compared to the SM group, SO(IO) has rank = 5 
(the rank gives the number of simultaneously diagonal generators; SU (iV) 
has rank = N — 1 and SO(2iV) has rank = iV). Breaking SO(IO) down 
reduces the group rank, which induces non-vanishing D-terms breaking 
SUSY at the same scale. To avoid broken SUSY at a high scale, one in¬ 
troduces another Higgs multiplet in the conjugated representation of the 
rank-reducing multiplet (126^) that cancels the other vev and keep SUSY 
intact [121]; additionally the extra field is needed to cancel chiral anoma¬ 
lies [122]. The minimal SUSY SO(10) Higgs content responsible for GUT 
breaking, neutrino masses and electroweak breaking would then consist 
of 210^ breaking SO(IO) to a partially unified group; 126^ and 126^ 
responsible for neutrino Majorana masses; and 10^ which contains the 
two MSSM Higgs doublets [121, 123, 124]. The minimal SUSY SO(IO) 
model still is in quite a good shape if RG corrections are included into 
a fit of flavor data (quark and lepton) [125] where it was disfavored by 
the fit without RGE [126]. Despite the large representations, the gauge 
coupling stays perturbative even beyond the Planck scale if threshold and 
gravitational corrections are taken into account [127]. 

2.3.2 Partial Unification 


When the hurlyburly's done. 

—Second Witch [William Shakespeare, Macbeth] 

SO(10) as a framework gives an excellent playground to study partial 
unification: an enlarged gauge group as it appears by deconstruction of 
the GUT group via symmetry breaking in the top-down approach. The 
partially unified picture allows to restore parity, e. g. gives an explana¬ 
tion why the weak interaction only couples to left-handed fermions, and 
simultaneously generates Majorana masses for neutrinos [129]. This in¬ 
termediate symmetry is known as left-right symmetry, SU(3)(, x SU(2 )l x 
SU(2)r X U(l)B_i and additionally gauges the B — L number. A slightly 
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The coupling to the 
120[f gives no 
relation between 
down type fermions 
[119]. 


An introduction to 
seesaw mechanisms 
in given in Sec. 3.1. 

We cannot discuss 
SO(IO) in detail 
but rather want to 
point out a viable 
framework, which 
gives seesaw 
neutrinos “for free”. 


In view of no SUSY 
particles at the LHC, 
one may wonder 
whether it is not too 
bad to break SUSY 
at a high scale. 


Interestingly, the 
fits of [125] prefer 
the non-minimal 
model with an 
additional 120^. 


The SUSY version 
and its connection 
to the more general 
picture of radiative 
corrections to 
neutrino mixing has 
been studied 
in [128]. 
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more unified group is SU(4)c x SU(2 )l x SU(2)r which also can be bro¬ 
ken out of SO (10) [130] and unifies quarks and leptons below the GUT 
scale where lepton number advances to the “fourth color” [131]. 


In case of an 
inverted mass 
spectrum (where 
neutrino number 3 
is the lightest) a 
slightly different 
Arjigj is found. 
Actually, the sign of 
Arjijj is still 
unknown which 
gives the ambiguity. 

The choice of three 
right-handed 
neutrinos is done 
on symmetric 
grounds, moreover 
motivated by 
SO(IO). 


Remarks about 
spinor notation in 
App. A. 


Approximate 
unitary means 


l7'T7 = l-FO(p2). 


2.4 NEUTRINOS IN THE STANDARD MODEL AND BEYOND 


„Liebe Radioaktive Damen und Herren" 
—Wolfgang Pauli, Dec 4 1930 


As pointed out in the introduction to the SM, Sec. 2.1, neutrinos are ex¬ 
actly massless in the model. However, extensions of the gauge sector as 
briefly mentioned in Sec. 2.3 naturally incorporate neutrino masses. On 
the one hand, it would be a puzzle if neutrinos were exactly massless. On 
the other hand, observations clearly contradict the SM in that point—we 
shortly refer to the Review of the Particle Data group and the appropriate 
references therein [23]. Experimentally, the physical mass squared differ¬ 
ences can be obtained 


Aml^ = 7.50+° ^ X 10“^ eV^ 

= 2.457± 0.047 x 10“^ eV^, 


where Am^. = and we restricted ourselves to the result of a 

j I 

normal hierarchy (ArUgg > 0) as follows from a global fit of neutrino 
oscillation data [132]. 

Actually, masses for neutrinos can be very simply added to the SM 
Yukawa Lagrangian that was given in Eq. (2.5), adding three right-handed 
neutrinos to the SM what we then call vSM,^ 


^tSM _ 




SM 

Yuk 


+ Y^jLui ■ Hvrj - + h. c. 


(2.37) 


The Majorana mass Mr for right-handed neutrinos can be added with¬ 
out harm, because they are gauge singlets anyway. Because of the same 
reason, its value is not restricted to the electroweak scale. In view of a 
partially unified scenario, we assume the scale Mr somewhat below the 
GUT scale, Mr 10^2 -+^ GeV. 

Deriving the neutrino mass matrix out of (2.37), one gets 


mass 





(2.38) 


where m° = is the Dirac mass matrix and the charge conjugated 
right-handed neutrino. We have switched to the less heavy Weyl spinor 
notation (and only deal with left-handed Weyl spinors); the Vj] are 3- 
vectors in flavor space. 

The mass matrix of Eq. (2.38) can be perturbatively diagonalized with 
an approximate unitary matrix 

1 The factor | in front of has to be included because + h. c. = . 
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with p = see e. g. [133], such that 








(2.39) 


where 

ml = Mr. 


(2.40) 


Mass matrices like the one of Eq. (2.39) follow directly from SO(IO) GU 
with an intermediate left-right symmetric breaking scale [134]. The left- 
handed Majorana mass can be achieved via couplings to an SU(2 )l triplet 
Higgs which acquires a small vev via a vev seesaw (see Sec. 3.1). 

The diagonalization of the light neutrino mass matrix ml determines 
a mixing matrix which shall play the same role as the CKM matrix in 
the quark sector. Neutrino mixing was first proposed by Pontecorvo [135] 
and further developed by Maki, Nakagawa and Sakata [136]; we refer to 
the leptonic mixing matrix thus as Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) matrix. However, differently to the quark case, the PMNS ma¬ 
trix is not the additional transformation needed to diagonalize the second 
Yukawa coupling. Without right-handed neutrinos, there is no lepton mix¬ 
ing and we can always diagonalize Y® (which defines the charged lepton 
basis). Doing so, we redefine the lepton fields 




(2.41a) 

(2.41b) 


analogously to Eq. (2.8). Without loss of generality, we can always choose 
Mr diagonal, (Sr^ Mr ', redefining 

^ = ^lij^RJ- C2-41C) 

All transformations are now fixed and stays an arbitrary matrix in 
flavor space—which can be smartly parametrized in terms of knowns and 
unknowns, see Eq. (3.7). The object we have to deal with is anyway not 
but ml, which is also a complex symmetric matrix. We then have 

ihi = I^pmns^l^pmns ~ diagonal. (2.42) 


The mixing matrix of the charged current is found to be indeed U 


PMNS 


withv[;^v" 


L,i 


L,i * 


4^ = -^W+eLr'^VL + h.c. 


s/2 

ig2 

'V2 


W^e'S^Y^ ( V [/' ■ 

V L) ‘■^PMNS'^LI 


(2.43) 


h. c. 


The diagonalization 
procedure is the 
same including a 
left-handed 
Majorana mass 
already at the 
Lagrangian level, so 
we put it there. 


MJ^ is a complex 
symmetric matrix 
that is diagonalized 
via Takagi diagonal¬ 
ization [137] with a 
unitary matrix S^. 


The primed and 
double-primed 
fields are in the 
mass basis. With 

sJ: [sty = b 

^PMNS in 

the W-vertex. Note 
that the PMNS 
matrix is defined 
“upside-down” 
compared to the 
CKM matrix of 
Eq. (2.9). 
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The magnitudes of 
PMNS elements can 
be displayed with 
the central values 


as lUpMJvsl — 

I m m • \ 


In SO(IO) the fields 
assigned to 
right-handed 
neutrinos are in the 
same representation 
as all the other 
matter fields. So it 
is necessary to have 
the same number. 


In general, ij —> liT 
with j > i. 


The PMNS matrix is parametrized conveniently in the same way as the 
CKM matrix in Eq. (2.10). Majorana neutrinos (since they are real) do not 
allow to absorb as many phases as Dirac fermions, so two more complex 
phases survive, I/pmns = ^ckm^ with a phase matrixP = diag(e^“i, 1) 
The absolute values are going to be determined with better and better pre¬ 
cision; within the 3cr intervals we have [132] 


ll^PMJVsl 


'^0.801... 0.845 
0.225... 0.517 
^0.246... 0.529 


0.514... 0.580 
0.441... 0.699 
0.464... 0.713 


0.137...0.158^ 
0.614... 0.793 
0.590... 0.776^ 


(2.44) 


NEUTRINOS MEET SUPERSYMMETRY: EXTENSION OF THE MSSM The 
most economic extension of the Minimal Supersymmetric Standard Model 
to incorporate massive neutrinos is the extension by right-handed neu¬ 
trino superfields. We supersymmetrize the vSM and call this vMSSM with 
the following superpotential 

W^mssm = Wmssm + (2.45) 

For aesthetic reasons, we introduce three right-handed neutrino super¬ 
fields Nr; = {^RiJ^Rj}) the same number as left-handed fields as moti¬ 
vated from SO(10) GU. 

Compared to the MSSM, the extension with right-handed neutrinos also 
comes along with additional soft SUSY breaking terms: one more soft 
mass matrix, a Higgs-sneutrino trilinear coupling and the sneutrino B- 
term. We add the following soft breaking Lagrangian: 






(4.1 • + (^ a ) i 


,/RA+h-c.). 


(2.46) 


We write the neutrino B-term in a way that suggests no connection to 
Mr although it can be seen as “Majorana-like” soft breaking mass (and 
therefore denoted here as to make clear that it carries mass dimension 
two). The usual way in the literature [138-141] is to write it down as 
B^ = byMj^ where by is a parameter of the SUSY scale. 

In general, flavor off-diagonal entries in the soft breaking contributions 
do influence observation of flavor changing neutral currents (FCNC) in 
charged lepton physics. Especially leptonic flavor violation in decays as 
ju —> ey has never been observed, so the tightest bounds on new SUSY 
contributions may kill most parameter configurations of the model. Never¬ 
theless, these constraints only affect the charged lepton sector: The flavor 
mixing parts of soft squared masses for the lepton doublet as well as the 
trilinear selectron coupling A® have to be negligible, at least for the first 
two generations. Bounds on 2-3 mixing are less stringent. If we impose 
minimal flavor violation in the charged sector and allow for large flavor¬ 
mixing contributions in the neutrino A-terms, lepton flavor violating FCNC 






2.4 NEUTRINOS IN THE STANDARD MODEL AND BEYOND 


23 


are safe. An aesthetic aspect of this description might be the underlying 
SUSY breaking mechanism, leading to MFV in the well-known part of the 
theory and somehow complete anarchy in the part, which is not accessible 
yet. At this point, we imply that such a mechanism is viable and leads to 
the observed amount of flavor mixing. 


We impose RFV in 
the lepton sector as 
described in 
Sec. 3.2. 


THE SNEUTRiNO SQUARED MASS MATRIX Due to the Majoraua struc¬ 
ture, the sneutrino mass matrix gets blown up—similar to the neutrino 
mass matrix in the non-supersymmetric vSM—and there are twelve phys¬ 
ical sneutrino mass eigenstates instead of only three in the MSSM, where¬ 
upon half of them are heavy—similar to the heavy (mostly right-handed) 
neutrinos: 




Ml*,* 

Ml*„* 

Mi*y 


Ml, 

Ml, 

Ml,* 

Ml,* 

Ml* 

Mrr* 

Ml 

Ml 


V-^R*L 

Mji*,* 

Mr*r* 

Mr*rJ 


^ Ml, 

Ml„' 



2 

^Ml^] 

' Ml,^ 

)' 



(2.47) 


in a basis v = '^r’'^r) > £mass _ where the 

individual 6x6 blocks have the following hierarchies in the orders of 
magnitude [140]: 


Note that in the 
MSSM there are no 
right-handed 
neutrinos. So the 
number of states is 
doubled twice. 


Each (X = L, R) 
is a 3-vector in 
flavor space. 



{Mi,y 


[ ^(^SUSy) ^{^SVSY^r)'\ 

^IrJ ^ V^(MsusyMr) 0{MI) J ’ 

(2.48) 


which has a similar hierarchy as the full neutrino mass matrix 


My ^ 


0 0(v) A 

0(v) O(M^) J' 


In an analogous treatment to the neutrino sector, we can approximately 
diagonalize Eq. (2.48) and get an effective light sneutrino squared mass 
matrix—the RR block does not change significantly. Especially, there is 
no distinct left-right mixing in the active sneutrino sector (because right- 
handed neutrinos and their scalar partners are heavy and integrated out 
well above the SUSY scale). However, a mixing of the left-handed partner 
fields with their complex conjugate is left, which contributes to Majorana 
mass corrections at one loop (see Sec. 3.2). 

The light sneutrino mass matrix has the following structure [140]: 

■^V( ~ ■^LL~ “1“ ^(^SUSY^R )’ (2.49) 

which provides a correction term to the MSSM sneutrino mass ^SUSy/-^R 
This term is absent, if there is no LR mixing in the sneutrino sector or the 
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The AL = 2 terms 
are the one 
~ which 

violate the global 
U(l)i charge. 


The full sneutrino 
mass matrix is 
derived in App. A. 


right-handed mass is sent to infinity. Especially it provides a seesatv-like 
connection between left-right mixing (i. e. trilinear couplings A.j,) and the 
heavy neutrino mass scale. Though this contribution ought to be small, it 
induces a mass splitting of order of the light neutrino masses. 

Performing the perturbative diagonalization, we find lepton number vi¬ 
olating terms in the 6x6 light sneutrino squared mass matrix: 



WL=2 


(Kl^oTJ ’ 


(2.50) 


where the AL = 0 block preserves total lepton number, while generation 
mixing is allowed, and the AL = 2 block violates lepton number by two 
units. 

Explicitly, the entries of the 3x3 sub-matrices are given by [140]: 


= mj + -M‘^ cos213+ (2.51a) 

- m>R {Ml + miy" , 

^lL=2 = ^T^R[^i+{^iy] <:2.51b) 

+ X>f [Ml + miy^ 

- 2mf Mr [mI + (m2)T] {bI) (m^ + m^) M^my 

+ O (a^susy^r ) ’ 
where = —jucot^m°* + 

Diagonalizing A4?^ of Eq. (2.50) yields the six physical light sneutrino 
mass eigenvalues, which are pairwise degenerate. In the literature for 
the one generation case [142] as well as for the general case [140] it 
is proposed to transform into the CP eigenbasis and to deal with real, self¬ 
conjugate mass eigenstates. To perform this transformation, we use 





such that A4?^ = V^M^V now is in the CP basis. If diagonalizes the 
matrix A4?^, so does with respect to A4?^. 

In this basis, the Eeynman rules take a particularly convenient form, 
where one only has to evaluate “half” of the mixing matrix, since Zj^^^ ^ = 
for i = 1,..., 3 and s = 1,..., 6. The Eeynman rules and details con¬ 
cerning the mixing matrices are given in App. A. 

Anyhow, for our analysis, we only perform a numerical diagonalization 
for which the perturbative approach is an overkill and may be used to un¬ 
derstand the structures behind. In principle, we can directly diagonalize 
the full 12 X 12 sneutrino mass matrix. Numerical cancellations and in¬ 
stabilities can be avoided using higher working precision. Still, it is more 
convenient to work in the effective theory with only light sneutrinos and 
Eq. (2.50). Discussing the anatomy of flavor changing contributions, we 
shall later switch to the full theory. 


VARIANTS OF NEUTRINO FLAVOR PHYSICS 



Physics of neutrino masses is physics beyond the SM. The SM per se has no 
room for massive neutrinos. Weinberg’s “Model of Leptons” [5] is a min¬ 
imal model describing lepton physics back in the sixties. However, since 
the observation of neutrino oscillations, we know that this simplest model 
cannot be true. In the SM there are no right-handed neutrinos, Majorana 
mass terms among left-handed fields are forbidden by gauge symmetry. 
We discuss several possibilities how to generate effectively a left-handed 
Majorana mass term at tree-level in Sec. 3.1 respecting the gauge struc¬ 
ture of the SM. In any case, the following effective operator is the only 
possible dimension-five operator [149] 

^dims = ^(irH)c(H-L^.), (3.1) 

where the dot product denotes SU(2)-invariant multiplication and C the 
charge conjugation matrix. The interaction is suppressed by a heavy scale 
A that is not restricted by the physics of the SM. In this construction, 
the neutrino mass matrix is given by the combination v^A/ (2A), where 
the Ajj are dimensionless couplings that in general mix flavor. Eq. (3.1) 
violates explicitly lepton number (AL = 2) and can be tested by the ob¬ 
servation of neutrinoless double beta decay [150-152]. The presence of 
the non-renormalizable term (3.1) is sufficient to explain neutrino masses 
within the field content of the SM [153], however, to build a UV complete 
theory so-called seesaw mechanisms were elaborated [129, 154-157]. 

The existence of such an operator with an a priori arbitrary flavor struc¬ 
ture sets the stage for lepton flavor physics. Without loss of generality, we 
work in the basis where the charged lepton Yukawa couplings are diag¬ 
onal and get the PMNS matrix from the diagonalization of the neutrino 
mass matrix only: 

^ = Kmns ^^pmns = diagonal. (3.2) 

Note that A in general is a complex symmetric matrix, so is a uni¬ 

tary matrix. The eigenvalues (proportional to the masses) can be defined 
complex, depending on the proper definition of UpM^g (see Secs. 2.1.3 
and 2.4). There are additional phases that cannot be absorbed into redef¬ 
initions of the fields, which would be possible for Dirac neutrinos [158- 
160]. In the presence of the Weinberg Operator (3.1), Neutrinos are Majo¬ 
rana fermions . It is a matter of taste whether to choose complex masses 
(as done e. g. in Sec. 3.2.1) or to assign Majorana phases to the mixing 
matrix (as discussed in Sec. 2.4). The latter choice makes it obvious that 
Majorana phases can never be observed in neutrino oscillations [160]. 


At that time, there 
was no solar 
neutrino 

problem [143-148] 
but the MNS matrix 
was already 
proposed [136], 


Eq. (3.1) is usually 
called “Weinberg 
Operator” since 
Weinberg 
introduced it. 

In Dirac space 
C = ir2ro andC 
has the properties 
C'i' = C'^=C"1 = 
-C. 


We refer to this 
choice as charged 
lepton basis or 
interaction basis, 
because for charged 
leptons mass and 
interaction states 
then are the same. 

This 

diagonaliza tion 
procedure is known 
as Takagi diagonal¬ 
ization [137]. 
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There are also 
possibilities to 
generate neutrino 
Majorana masses at 
the loop level, see 
e.g. [161, and 
references therein]. 



Because of the 
presence of the 
Dirac Yukawa 
coupling to 
left-handed leptons, 
we call vr 
right-handed 
neutrinos. 


0-2 = 


0 1 
1 0 

0 -i 
i 0 


A rough estimate gives a glance at the order of magnitude of the high 
scale A: let us assume the couplings A; ~ 0{1) and the resulting neu¬ 
trino masses ~ 0(0.1 eV), then with v ~ O(lOOGeV) we get A ~ 
0(10^3 GeV): 


0.1eV = 0.1 


(100GeV)2 
10^3 GeV 


= 0.1 


lO'^-lO^eV 

10^3 


3.1 NEUTRINOS AT TREE-LEVEL 

There are three possible ways to construct a UV complete theory lead¬ 
ing to the operator (3.1), which can be seen by rewriting it in SU(2)- 
invariant ways and inserting the missing multiplets into the effective op¬ 
erator [162]. 

TYPE I We introduce the SU(2) metric e to express the invariant mul¬ 
tiplication by matrix products (Lorentz-invariant multiplication in Dirac- 
space is understood without further notation) 

and find as candidate for the full Lagrangian 

A = A Cvr + (v^)^Cvr. (3.3) 

The additional fields Vr have to be singlets under the SM gauge group in 
order to give (3.1) after having integrated them out and couple to left- 
handed leptons via a Dirac-like Yukawa coupling Yy [163]. Singlets are 
not protected by any gauge symmetry and can acquire a Majorana mass 
term as shown in Eq. (3.3). The mechanism behind this Majorana mass is 
unknown or not specified, especially it is not related to electroweak sym¬ 
metry breaking, therefore unrestricted by scale considerations. In the case 
of only one right-handed neutrino, Mr is identical to the scale A. If there is 
a stronger hierarchy in right-handed masses, the Vr; have to be integrated 
out successively [164, 165]. For symmetry reasons, we shall use the type I 
seesaw mechanism with three right-handed neutrinos, motivated e. g. by 
SO(IO) GUTS. 

TYPE II A similar deconstruction can be done using a scalar triplet mo¬ 
tivated by an embedding into left-right symmetric models [129, 154, 157, 
166]: 

where the cti are the generators of SU(2). The UV complete theory can be 
built with a scalar triplet 

2 5° 


0-3 = 


1 0 
0 -1 
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coupling both to the lepton and Higgs doublets such that 

Ai = (iAcTiCL) + [H^eUiH) + M^5*5i. (3.4) 


The inclusion of scalar triplets acquiring a vev has to be treated with 
care since triplet vevs spoil one famous and important relation of the SM, 
the p-parameter [167] 

^ M| COS^0VE 

At tree-level in the SM, p = 1, but radiative corrections show deviations 
from that prediction even at one-loop [168]. Since the p-parameter is 
precisely measured (p = 1.00040 ± 0.00024) [23] and calculated with 
high precision in the SM [169], any deviation induced by triplet vevs has 
to be small in a sense that the triplet vev itself has to be small. On the 
other hand, if a Higgs triplet of SU(2 )l generates a Majorana mass for 
left-handed neutrinos at tree-level, its vev is restricted to be rather small 
compared to the electroweak scale anyway. The Higgs triplet ((5°) = v^) 
e. g. alters the relation to 


M, 


P = 


w 


M| cos^ 0. 


= 1 


4vi 


w 


for general scalar representations acquiring vevs V; the p-parameter is 
given by 


, E.(4rKA + i)-3y,2)|v,|2q 


where Y; and T; are hypercharge and weak isospin of the i-th Higgs multi- 
plet and C; = 1 for complex and Cj = | for real representations, see [170] 
and [171]. 


TYPE III The third variant needs to introduce exotic fermions (calling 
singlets not exotic), namely triplets under SU(2) [172]. 

^dims ~ -I C (LAcJiH) . 

The triplet fermions f = TiCJ^ again are Majorana fermions coupling to 
the standard leptons via an appropriate Yukawa coupling: 

£jjj = Yj. (L^acTjH) CTi+MjTjCTi. (3.5) 

While seesaws of type I and II can be naturally incorporated in popular 
extensions of the SM and follow immediately from common grand unified 
scenarios as SO(IO) [114, 119, 120, 173, 174], the triplet fermion exten¬ 
sion follows a different philosophy and shall not be considered any further 
in this thesis. 


v's /v 
Y 


A 



Which is the case 
for triplets making 
Majorana neutrino 
masses [158]. 
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3.2 RADIATIVE NEUTRINO MIXING 

The tree-level formulation of neutrino masses and mixing is in general suf¬ 
ficient to explain a non-trivial mixing pattern. Take e. g. the type I seesaw 
mechanism: working with a diagonal right-handed mass matrix Mr, the 
neutrino Yukawa coupling cannot be constrained and stays a rather arbi¬ 
trary matrix. Inverting the decomposition of the effective mass operator 


= -v^Y^M-^Yl (3.6) 

yields the famous Casas-Ibarra relation [175] 

(3.7) 


Generalizations 
with an arbitrary- 
number of singlet 
neutrinos are 
straigh tforward, 
see [158], though 
they do not permit 
any relation like 
Eq. (3.7). 


with an (arbitrary) complex orthogonal matrix TZ. The diagonal Matrix k 
contains the light neutrino masses, my = v^k, and TZ in general has three 
complex mixing angles. Together with three right-handed masses, there 
are nine free parameters that do not change the low-energy phenomenol¬ 
ogy. Moreover, Eq. (3.7) is only a reparameterization of the unknown— 
there is no hint to the origin of the mixing angles in I/pmns- 

Radiative corrections to neutrino masses and mixing are generically less 
intensively studied than tree-level realizations. On the one hand, quantum 
effects from the renormalization group (RG) evolution play an important 
role comparing low-energy observables with flavor models at the high 
scale. On the other hand, threshold corrections at an intermediate scale 
bring completely different aspects of neutrino flavor into the game. 

Threshold corrections and RG corrections interplay in the determina¬ 
tion of the corrected neutrino mass matrix 


^cb + lAcm^cB > <^ 3 . 8 ) 

where denote the corrections and is the tree-level mass ma¬ 

trix [176-179] . Gapital indices are in the interaction or flavor basis (A, B = 
e, p, t) and I = is the sum of RG and threshold corrections. The 

threshold corrections can be calculated diagrammatically via self-energy 
diagrams as discussed in Sec. 3.2.1, the RG corrections are to be obtained 
from the integration of the renormalization group equation for the neu¬ 
trino mass operator [179-183]. 


BRIEF DISCUSSION OF RG EFFECTS The Contributions from the renor¬ 
malization group are known to give a sizable effect for quasi-degenerate 
neutrino masses [180, 184-189]. Especially the choice of the same GP 
parity for two mass eigenstates may lead to large mixing at low scale 
irrespective of the original mixing at the high scale [190, 191], known 
as infrared fixed points [189]. The effect from the renormalization group 
severely depends on the Majorana phases: for a vanishing Majorana phase, 
maximal mixing patterns get diluted on the way to the high scale [192] 
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for quasi-degenerate (mo ~ 0(1 eV)) neutrino masses. If, on the contrary, 
the phase is large or the overall mass scale is much smaller than 1 eV, 
maximal mixing is preserved. Likewise, zero mixing (as follows from the 
assignment = m 2 = m 3 ) is conserved [193] irrespective of the Ma- 
jorana phase difference — a 2 |. We therefore neglect contributions that 
preserve specific mixing patterns for the low-energy threshold corrections. 

In the following sections, we shall discuss different aspects of threshold 
corrections to neutrino masses and mixing. First, we start with an exactly 
degenerate mass pattern at the tree-level and figure out, whether thresh¬ 
old corrections have the power to generate the observed differences in 
mass squares and the non-trivial mixing. Exact degeneracy comes along 
with trivial mixing. Threshold corrections both lift the degenerate masses 
and therewith mix different interaction states. Second, we discuss a sce¬ 
nario implementing a non-trivial mixing pattern already at tree-level with 
threshold corrections modifying the seesaw mass of Eq. (3.6) at the loop 
level. Third, we renormalize the mixing matrix directly and resum the en¬ 
hanced contributions ~ niy.S'^j/ Am^j in case of quasi-degenerate spectra. 

3.2.1 The Case of Degenerate Neutrino Masses 

Degenerate neutrino masses are probably an amusing gimmick of nature. 
Where in the early times of the SM no masses for the electrically neutral 
fermions were foreseen, the observed flavor oscillations provide sizable 
but small mass differences of the individual neutrino species. Any direct 
measurement of neutrino mass still lacks the discovery [23, 195-198]. If 
the overall mass scale is much larger than the mass differences, the 
neutrino spectrum is quasi-degenerate. The masses are easily calculated 
and expanded in Amfj/ 

(3.9a) 

\my^ \ = ^++ -Am^^, (3.9b) 

\niy^ \ = \j{niy^)^ + TUy^ + -Amg^. (3.9c) 

The measurement of the “effective electron neutrino mass” which is done 
via the endpoint of the f spectrum gives a robust determination of 

|We; 1^ Am^^ (normal hierarchy), 

i 

where the indirect measurement using neutrinoless double-/)-decay may 
suffer from destructive interference due to Majorana phases (or other 
new physics). Complementary to measurements in the laboratory, esti¬ 
mates or bounds on the neutrino mass can be achieved by cosmologi¬ 
cal observations. Cosmology constrains the sum of relativistic neutrino 
masses (and by the same time counts the number of active neu¬ 

trino states) but the bounds are not quite robust. Let us take a look into 


29 


In the last two 
scenarios, the 
masses may not be 
degenerate. 


The content and the 
results of this 
section was 
published in [194]. 


We define the mass 
square differences 
Arn^- = rr? — rrP" 

IJ Vi Vj 

and only discuss the 
normal hierarchy. 
For inverted 
hierarchy, the 
lightest neutrino 

( 0 ) 

mass IS = my . 

^|Wei|^mf, with 
i 

U the neutrino 
mixing matrix. 
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Majorana neutrinos 
require SO (3). 


Note that exactly 
degenerate 
neutrinos have no 
mixing at tree-level. 


the Planck report [199]: depending on the fit model, several bounds are 
proposed; the strongest is < 0.23 eV, relaxing spatial flatness one 

finds < 0.32 eV. In any case, all analyses are based on the standard 

cosmological ACDM model. 

While a positive direct neutrino mass measurement in the near future 
will immediately put us into the quasi-degenerate regime, cosmology dis¬ 
favors this possibility with increasing significance. In any case, quasi-de¬ 
generate neutrino masses need a cautious treatment from the flavor sym¬ 
metry point of view, see e. g. [200-203]. Exact degeneracy, however, is 
a direct consequence of SO(3) or SU(3) invariance. We now want to 
keep the fundamental flavor symmetry for neutrinos intact and figure out 
whether radiative breaking has the power to produce the observed devia¬ 
tions from degeneracy and simultaneously the mixing matrix. 

Corrections to degenerate masses can be treated very easily performing 
a rediagonalization of the corrected mass matrix (3.8). In general, the tree- 
level mass matrix is not diagonal in the flavor basis. Transforming 
into the mass eigenbasis (at tree-level) using the tree-level mixing matrix 
results in a non-diagonal corrected mass matrix whose off-diagonal 
elements stem from the off-diagonal threshold corrections: 


^ab = 


( 0 ) 


5ab + {ni, 


(0) _L^(0) 


m 


n T 

J ^ab’ 


(3.10) 


where small indices a, b now are meant to be in the mass basis and 


T -Vt 

^ab - 2^^AB^Aa Bb ' 
AB 


Eq. (3.10) reveals two interesting observations: first, if any = 
= m, simplifies tremendously to (e. g. = —m^ 2 ^ = nig^^) 


r 


m = m 


1 + 2UaiUbiIab 
0 

‘^UaiUbbIab 


0 2UaiUb3Iab ^ 

~ 2Ua2Ub2^AB 0 

0 1 +2UA3UB3lABy 


(3.II) 


and there is only one off-diagonal entry left (remember that Majorana 
mass matrices are symmetric)—^which can be eliminated using one free ro¬ 
tation [176,177, 179]. Second, in the case the observed fla¬ 

vor mixing is directly a result of non-universal (and flavor non-diagonal) 
threshold corrections. Eor exact degeneracy, there are three free rotations 
which can be used to diagonalize I [194]. 


CP PHASES AND MAJORANA NEUTRINOS The case of Majorana neu¬ 
trinos does not allow us to rotate away as many CP phases as for Dirac 
fermions. In general, there are two more phases left, such that the com¬ 
plete diagonalization matrix can be written as a product of a unitary ma¬ 
trix with three angles and one phase and a phase matrix: Uy = with 
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P = diag i). xhe Majorana phases ai 2 can then be absorbed in 

a redefinition of the masses instead of a redefinition of the fields: 

m(o) = mo diag 1 ). 


We have chosen the phases in a particular way to have a real and positive 
m 3 . 

Under the assumption of CP conservation in the Majorana phases, we 
take ai 2 ^ {0,±^} and assign the relative CP parity of the respective 

u -1 ( 0 ) ( 0 ) ( 0 ) 

mass eigenstate to the mass eigenvalue, so e. g. mj = —m^ = as 
discussed in the following. 

THE CASE = —m^ 2 ^ = mg^^ The two-fold degeneracy leaves a 

freedom of rotation which is in the 1-3 plane (U —> UR 13 with a real 
rotation matrix and two fixed angles in [/) and the full diagonalization 
of in Eq. (3.11) can be done solving the equation 

YjUaiUb3Iab = 0, (3.12) 

A,B 

which, for fixed mixing angles in U can be done by choosing appropriate 
1^. An interesting exercise now is to find out, how many nonzero contri¬ 
butions in I are needed to fully reproduce the masses. The squared mass 
differences can be easily calculated from Eq. (3.11) 


Examination of 
other configurations 

(0) (0) 

as m\ =m\ = 

(0) (0) 

— m' or —m' = 

(0) (0) 

rrig — rrig are 
qualitatively the 
same and can be 
treated analogously, 
which is not done 
here because there 
are basically no new 
insights from other 
sign assignments. 


Am^, 

ab 


(1 +2[4,UBa4B) - {l+2UMblAB) 


^2 


(3.13) 


Shifting all corrections by an overall flavor universal constant, I ab = 

Iab ~ Iq^ab’ does neither change the mixing angles nor does it affect the 
ratio Am^i/ Arn^i [179]. The effect on the mass parameter m factors out 
and manifests itself also as a shift and can be absorbed by a redefinition: 
m = (1 + 21o)tti. 

With only flavor diagonal threshold corrections 1aa = ^a — ^e) we 
can try to fix the third, free mixing angle 613 in terms of the other two 
and the 1 ^: 


•5i3 — ^ 23^23 


^12 

C 12 le “ '^23^11 ~ ^23^T 


(3.14) 


which is a generalization of the result from [179] which was given for 
= 0. Let us now take the mixing angles as input values and determine 
one of the flavor diagonal corrections in terms of the others 

T 2 r I 2 r I ^12^23^23 f r t \ /'■-i i 

“ •^23^11 ^23^T d- ■ (3.15) 

C12'5i3 

It is, however, not possible to fit the large ratio Amgg/Amgg 33 to¬ 
gether with the rather large value 5^3 0.15 instead of @33 0° without 

introducing non-perturbative values of even for a highly degenerate 
neutrino mass spectrum. The situation changes in the presence of the 


As defined in 
Sec. 2.1.3 we again 
use the shortcuts 
sij = sin Oij and 
cij = cos 6ij. 


Solutions exist for 
some U > 1. 
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Figure 2; The allowed ranges for and for and within their 1 a ranges (the mixing 
angles are taken at the central values). The left plot shows the dependence on the lightest 
neutrino mass tuq in the vertical direction. There is one class of solutions where all three 
non-vanishing elements of / are close to zero. (Taken from [194].) 


We obtain the 
corresponding 
result 
from [176, 179], 

^13 = 

— tan 6-12 cot 2023 m 
the limit —>■ 0. 


Dirac CP phase ^ 0, of course. Then the situation also gets much 
more complicated. An easier and not uninteresting option is to study the 
influence of flavor changing threshold corrections. With a dominant 
we can indeed find viable solutions with a reasonable 6 ^ 3 , 

(V 23523 --f pT 2 023) tan 012 

5i3 — -^-;-—-, (3.16) 

^^^23 "^M-T ^tn 2023 -^e 

and the correct masses even for not too heavy (i. e. not too degenerate) 
neutrinos (m 0.1 eV). The interesting results are shown in Fig. 2. There 
is one class of solutions where all non-zero /^ are close to zero and give 
the correct and 033 . 


We are working 
with Majorana 
neutrinos, therefore 
the mass matrix 

as well as the 
threshold correction 
matrix I is 
symmetric. 


( 0 ) ( 0 ) ( 0 ) , , 

THE EXACT DEGENERATE SITUATION! m] = Let US as¬ 

sign the same CP parities to all three neutrino mass eigenstates—Eq. (3.10) 
does not show any zero entry, but obviously can be diagonalized by diag¬ 
onalizing the perturbation / only: 



^hi 

h2 

h3 

= ml+ m 

h2 

^22 

^23 


i^h3 

^23 

h3 j 


(3.17) 


and m is the common neutrino mass. 

We can now play an amusing game: Experimentally, 023, the atmo¬ 
spheric mixing angle, is measured to be roughly maximal I®231 ^ f with a 
small deviation of a few degrees. A maximal mixing in the 2-3 plane can 
be achieved via the rotation matrix 


U23 = 


fl 

0 


0 0 ^ 

1 1 

V2 42 

1 1 

' 72 727 
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And up to recently, all measurements of the third mixing angle @13 were 
consistent with zero. Taking these two phenomenological observations as 
starting point, we arrive at a determination of 632 in terms of I 22 and 
/12 only 

1 f 2 ^ 2/12 A 

012 - arctan —-— , (3.18) 

^ “■‘11 J 


where we have exploited @13 0 in order to approximate /13 ^ 1 12 - After 

performing the 2-3 rotation with U 23 , we are left with 


I' = uyu^23 


r 


hi 

tl2+h3 


42 

2-^22 


42 


(l2_(l3 ^ 

42 

0 


0 


y 


(3.19) 


The eigenvalues of m'^ can also be calculated in terms of the same three 
Ijj with m 3 = m. Altogether, there are four free parameters left (m. In, I 22 
and I 12 ) required for fitting three masses and one mixing angle ( 0 i 2 )- The 
other two mixing angles were set to phenomenologically motivated dis¬ 
tinct values (013 = 0 and 623 = 7 i/ 4 ) and shall receive small corrections 
in the following. 

We relax the restrictions required for 623 = n/4 (I 33 = I 22 ) nnd 0 i 3 = 0 
(li 3 = I 12 ) in a way that we first parametrize deviations: 


■^33 — h2 + 
■^13 — -^12 + 


(3.20) 


The matrix of threshold corrections is then written as 



^ hi 

■^12 

I 12 + 5 

I = 

h2 

h2 

hs 


yh2 + 5 

hs 

h2~^ ^ y 


(3.21) 


where we have also lifted the artificial requirement I 23 = I 22 which had 
no influence on 623 before anyway. We also need the full freedom of all fla¬ 
vor non-diagonal corrections to fit three mixing angles and three masses 
(namely Iii,l 22 >^i 2 >^ 23 > ^ ^^d e) and assign the “unperturbed” mass pa¬ 
rameter m to be the lightest neutrino mass m = . Any flavor-universal 

contribution in the threshold corrections can be again simply added as a 
shift in the diagonals: Iy^ = — Iq for A = e, p, t. 

As a proof of principle, we perform a brief numerical analysis for two 
benchmark scenarios, where one corresponds to a possible discovery of 
neutrino mass at the Katrin experiment (m^^^ = 0.35 eV) and the other 
one lies at the lower edge of quasi-degenerate neutrino masses (but still 
allowed by the tightest ACDM cosmology bounds, = 0.1 eV). The 
results are shown in Tab. 1. It is gratifying to see that the entries of I are in 
both cases of the size of a typical radiative correction < 0 ( 1 / 100 ) (note 


33 


where recently 
means < 3 years: 
the first (high 
sigma) non-zero 
measurements date 
from late 2012 
[204]. 


We implicitly set 
the whole 2-3 block 
to the same values: 
hs = ^33 = hi- 


The Katrin 
collaboration states 
the possibility of a 
discovery with 

> 0.35 eV and 
a 95% exclusion 
with < 0.2 eV 
[197]. 











34 VARIANTS OF NEUTRINO FLAVOR PHYSICS 


Table 1: Values of the threshold corrections needed to obtain the observed mixing 
angles and mass splittings for a common neutrino mass of 0.1 eV and 
0.35 eV. 


mo = 0.1eV mo = 0.35eV 


^11 


3.54 X 

CO 

1 

o 

1 


3.00 X 

1 

o 

t—1 

'^12 


1.19 X 

10 “^ 


1.02 X 

10 “^ 

22 


4.67 X 

10“2 


4.01 X 

10 “^ 

23 


5.43 X 

10“2 


4.67 X 

10 “^ 

e 


2.28 X 

10“2 


1.96 X 

10“2 

5 

^0.354 

6.73 X 

10 “^ 

^0.300 

1.56 X 

10 “^ 


1.19 

1.20 


1.02 

1.03 

I 

1.19 

4.67 

5.43 

X 10“^ 

1.02 

4.01 

4.67 


^ 1.20 

5.43 

6.96^ 


^ 1.03 

4.67 

5.97 


For radiative 
neutrino mixing, on 
one hand one has to 
ensure that there is 
no mixing at the 
tree-level. On the 
other hand, 
working with 
degenerate 
neutrinos imposing 
SO(3) orSU(3) 
flavor symmetries 
fulfills this 
requirement 
automatically. 


that we only tvant to generate tiny deviations from the degenerate pattern 
in a regime where the physical masses are only slightly non-degenerate) 
and show a hierarchy as 1 < 2 < 3 for labeling the generations. This 
observation can be used in any new physics model with flavor changing 
low-energy threshold corrections. 

A crucial point in the discussion is the behavior of the generic thresh¬ 
old corrections with the lightest neutrino mass. In case the overall mass 
scale drops below 0.1 eV, the spectrum loses the degeneracy property 
which is reflected in values ~ 0.1 as can be seen in Fig. 3. Corrections 
are needed that are not of the size of typical perturbative corrections. The 
hierarchical regime (m^y^ <C 0.1 eV) needs a special kind of flavor symme¬ 
try breaking where the degenerate patterns only needs a symmetry that 
guarantees equal masses. For a given symmetry breaking chain, the hierar¬ 
chy can be exploited to construct the mixing matrix out of the mass ratios 
(see Chapter 6 and Ref. [205]). 

THRESHOLD CORRECTIONS IN THE VMSSM: NMFV Non-minimal Fla¬ 
vor Violation (NMFV) in the soft SUSY breaking terms allow to play an 
amusing game: is it possible to generate the observed neutrino flavor mix¬ 
ing via soft SUSY breaking? The philosophy behind this reasoning was 
explained in [206]. In this way, flavor violation enters by means of su¬ 
persymmetric loop corrections. Since the origin of SUSY breaking stays 
unknown, we connect the flavor puzzle in the SM to the “hidden sector”. 

The superpotential of the vMSSM was introduced in Sec 2.4, Eq. (2.45). 
With the term “vMSSM” we denote the Minimal Supersymmetric Stan¬ 
dard Model (MSSM) extended with right-handed neutrinos and no further 
specification of a UV theory, especially the right-handed Majorana masses 
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Figure 3: Graphical representation of the individual threshold corrections to the neutrino mass. The 
first row shows the dependence on the absolute neutrino mass scale: the larger the 
smaller the corrections can be. The upper right plot shows the deviations from equal values: 
5 = /i3 — /i2 and e = 133 —122 es defined in Eq. (3.20). The lower line shows the interplay 
of the individual Ijy. compared to 1 23 , similar plots can be done for other combinations. The 
lower right plot shows the normalized to the largest contribution 733. (Taken from [194].) 


are just present irrespective of any symmetry breaking mechanism which 
generates them. We assume degenerate soft breaking masses, 


th^ — — Mgygyl, 


in order to avoid large leptonic FCNC observables as —» fjp or Ij —> 

(with (i,j = 1,2,3 and i < j). The flavor violating contribution 
then lies in the trilinear sneutrino-Higgs couplings only, if we ignore 
the flavor structure of the neutrino B-term. 

The SUSY threshold corrections to the neutrino mass matrix can be 
calculated in terms of neutrino self-energies [140] 


The assignment 
ml = mj is only 
for convenience and 
not a necessary 
choice. 
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(3.22) 
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with the decomposition of the neutrino self-energy 


Ey\p) =E) 
1 ] / 1} 




sy 

i] 




(p")Pr+ 


sy 

o 


ip^)p^ 


(3.23) 


Eq. (3.23) is the 
application of 
decomposition 
(2.34) for Majorana 
fermions. 


Bq and Bi are the 
standard 
Passarino-Veltman 
two-point loop 
functions, see 
App. A. 


For Majorana neutrinos, the self-energy is flavor symmetric (27jj = 27^; and 
the coefficients in front of the left and right projectors (Pl and P^ respec¬ 
tively) are related via complex conjugation. The neutrino self-energies are 
evaluated at p^ = 0 because we can neglect the neutrino masses com¬ 
pared to the superheavy particles in the loop (the same is true for all 
supersymmetric corrections to SM fermion self-energies). 

The generic flavor changing self-energies have already been calculated 
in [128] and are in agreement with [140]. We are interested in the influ¬ 
ence of the soft breaking sneutrino parameters, so we give the sneutrino- 
chargino and -neutralino self-energy (mixing matrices defined in App. A) : 




v2 


(3.24a) 


4(27r)2""^'' 

{S 2 ^ 2 k ~ (^PMNs)i'i(t^PMNs)T;> 

1 




4(271)2 Zfc 


—1 


/2 


(3.24b) 


\s2zl - 


where summation over repeated indices is understood. 

We now look for suitable values of in order to produce the struc¬ 
ture of the generic threshold correction as elaborated above, i. e. solve the 
equation mljj = — mdjj in terms of AT. Because the depen¬ 

dence on the neutrino mass parameter is quite interesting, we vary 
the neutrino mass in a wider range, shown in Fig. 4. 

The results are qualitatively very stable under variation of the free SUSY 
The same is true for parameters. In any case, we need large neutrino A-terms to get the struc- 
the application in of the threshold corrections as for the generic discussion. Effectively, 

Sec. 3.2.3. combination A’^ / Mgysy drives the corrections. For the analysis pre¬ 

sented here we vary the values of the following variables randomly in the 
given intervals: 


Msusy ^ [500,5000] GeV, 

Ml £ [0.3,3] Msusy) 

Mse [1,5 ]Msusy> (3-25) 

pe [-15,15]TeV, 
tan/3 e [10,60]. 


As expected, for low values of the absolute neutrino mass ttiq where 
the deviation from the degenerate pattern is large, the SUSY threshold 
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2000 3000 

^SUSy/ 


5000 



0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

mo/eV 


100 


Figure 4: The left plot shows values of for the variation of parameters specified in (3.25). The 
lightest neutrino mass was chosen in the regime plotted on the right side, where we 
rescaled all trilinear soft breaking couplings with the SUSY scale, = Ajj /Msygy The values 
of ai 2 and 013 are roughly the same since they differ only by a small parameter as described 
in the generic discussion. (Taken from [194].) 



Figure 5: Correlations between the different elements of found to reproduce neutrino masses and 
mixings. The left plot shows the off-diagonals with respect to A '^2 where the right shows the 
correlation with the diagonal entries. (Taken from [194].) 


corrections measured in the values of have to be large as shovra in 
Fig. 4 where we plotted the ratio a^j = A'^j / MgusY- The left-hand side of 
Fig. 4 compared to the right-hand side shows that basically this ratio is the 
parameter which drives the corrections and has the same shape as the 
dependent on tuq where the size of the AT depending on the SUSY scale 
also is sensitive to the parameters of the theory. 

3.2.2 Renormalizing the Seesaw 

In the previous Section and in the following we discuss supersymmetric 
threshold corrections to neutrino mixing with NMFV soft breaking trilin¬ 
ear couplings. There is, however, a part of the corrections which is in¬ 
dependent of soft breaking terms (a genuine F-term contribution to the 
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The expression 
diag(j«:i) means a 
diagonal matrix 
with entry Xi in the 
i-th diagonal 
element. 

As will turn out, g 
is something like 

giy2/(i67r2). 


sneutrino mass matrix) and nevertheless alters the flavor structure of the 
renormalized light neutrino mass, if the heavy neutrinos are not degener¬ 
ate. This type of corrections gives a logarithmically enhanced contribution 
to the one-loop seesatv mass formula with a logarithm of the heavy neu¬ 
trino mass: 


Am^ ~ YyM-^ 


1 + g diag 


log 


M, 


SUSY 


1 


/ 

R,i y J 


(3.26) 


where g is some loop suppressed coupling factor. This case becomes impor¬ 
tant if the right-handed neutrinos are hierarchical. The influence on the 
flavor pattern of Eq. (3.26) disappears completely for degenerate right- 
handed neutrinos; in that case log(M|ygY/^ r)/^r is proportional to the 
identity matrix, and the full flavor structure of the seesaw mass term sits 
in the product of neutrino Yukawa couplings 

For the description of the effect, it is sufficient to discuss a reduced 
parameter set of the vMSSM with a left-right mixing in the sneutrino 
sector only from the F -term contribution to the mass matrix ~ : 


C ^IvuI"v*y,y;vl + IvuI"v*y;yt^r 

- VdP*VL^vVR + Vr + h. c. 

+ + Aoft> 

Aoft = ^Lm^VL + VR^v^R. 


(3.27a) 

(3.27b) 


(3.27c) 


and without soft breaking trilinear and bilinear couplings h^v*AyV^ and 
respectively. In the following, we also neglect the p-term contri¬ 
bution in the sneutrino mass matrix (which is multiplied with the smaller 
V(j anyway). 

The sneutrino squared mass matrix can then be expressed in terms of 
the right-handed Majorana mass Mr, the Dirac mass m° = ^Yy and 
the generic soft SUSY breaking mass mg = M^usy us simplify the 
discussion with m? = = m?l): 

t V S 





(3.28) 


Mil- 


^tR = 


■^RR ~ 






mt 
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'Mr 


dT D* 

^y 


-m°M* 

2 y Ry 


mt 


iMRM*+m°mf+ m2 




D* 
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(3.29) 
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note that = Mr. We have chosen a basis N = 


with 


R. we nave eiiuseii a uasis iv — ^ vr, 

3-vectors in flavor space suppressing generation indices. 

Due to the large hierarchy between mg and Mr, the “right-handed” su¬ 
perpartners of the neutrinos have basically the mass Mr, where the “left- 
handed” ones live at mg. 


CORRECTIONS TO TYPE I SEESAW The tree-level seesaw mass formula 
carries a structure where the inverse right-handed Majorana mass is sand¬ 
wiched between Dirac Yukawa couplings, YyM~'^Yy. The SUSY one- 
loop corrections to the seesaw mass operator carry the same sandwich-like 
structure, ~ v^Y.^, (something) y!J, but between the two Yukawa couplings 
something more happens. Two Dirac Yukawas are essential to include two 
chirality flips in order to intermediarily have right-handed (s)neutrinos 
and therefore the suppression with 1 / Mr in the loop. Moreover, we need 
one lepton flow flip that has to be induced by one Majorana mass insertion 
as shown in Fig. 6, such that both diagrams show a AL = 2 transition. We 
estimate the diagram to the right of Fig. 6 to be 

Si It 

~ Y Md- 

1671^ ^ ^ ml 

We calculate the loop diagram in the mass insertion approximation, 
where the fields running in the loop are interaction eigenstates and do 
not refer to mass eigenstates. The mixing occurs via the couplings. In a 
world without /i-term, there is no Higgsino mixing (and = 0) and the 
diagram simplifies accordingly: 


-iSjf = 

U 4 

flR 


f 


(27r; 




IK JK 

[(p-q)2-m2^J 
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^—*0 U^l 

y_ (q^ + ^mg)MR^,Y.^Y;, _ 

fc=l (q2-m?)q2 - q)2 _ J 

„2„2 .. n, ,, 
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(3.30) 

For vanishing external momenta, the loop integral reduces to Cq and 
the one-loop self-energy yields 


„2„2 Hr yxyv 

'■ 2-i m2 _m2 


6471^ ^ m? 

fc = l VRi 
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log 
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T<,i 


mt f m? 
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(3.31) 
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Figure 6 ; Tree-level plus one-loop contribution to the type I seesatv mass formula for the light neutrinos. 


The mass splitting 
of the light 
sneutrinos is of the 
order of the mass 
splitting of light 
neutrinos—and 
therefore 
completely 
negligible 
compared to the 
SUSY scale mg. 


The tree-level 
mixing matrix may 
be governed by 
flavor symmetries. 


where we have used standard conventions for the loop functions (our defi¬ 
nition of Cq is given in App. A), approximated Mr expanded 

in small mass ratios like mg/m?^, where mg is either mg or 

The light sneutrinos can easily be taken degenerate, since their masses 
are dominated by the soft SUSY breaking mass parameters and deviations 
are of the order of the light neutrino mass squared. The heavier sneutrinos 
get their mass basically from the Majorana mass term. 

However, by assuming a similar hierarchy in the neutrino Yukawa cou¬ 
pling as for charged particles’ Yukawa couplings (quarks and leptons), to 
end up with a quasi-degenerate light neutrino mass spectrum, the heavy 
eigenvalues can differ over several orders of magnitude. 

In the case of degenerate SUSY masses, i. e. the limit rriy^ —> mg —> mg, 
Eq. (3.31) reduces to a logarithmically enhanced contribution to the tree- 
level neutrino mass operator k = YyM~^Y^: 


Ak 

^ 6471^ 


1 -V-VI.R ± y, 


(3.32) 


where is a diagonal matrix with the entries log j ^Vr*.- 

Up to small corrections 0{m^), niy^ = Mr. Therefore, one has to rediago¬ 
nalize the neutrino mass matrix 


m. 


= v2(k:^ + Ak^ 

f 


= vlYy diag 


log(ms/m^,J 


V 


m 




gi 

6471^ m 


A 




(3.33) 


which for exactly degenerate right-handed masses m^^^ ^ has 

the same flavor structure as the tree-level version: ~ In general, the 

tree-level PMNS matrix diagonalizes the combination YyM^ Yy 

which gets altered by the logarithm in the diagonal matrix in between 
the Yukawas. The physical PMNS matrix therefore can be significantly 
changed by the logarithmic structure and deviate from any pre-specified 
pattern. 
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Figure 7: Scatter plots for the corrected neutrino mixing angles 6jj = 6i2> ^13> ®23 for quasi-degenerate 
right-handed masses (left panel) and a strong hierarchy in the heavy spectrum (right panel). 
An interesting observation is that 6^2 nearly scatters over the complete range, especially for 
Yukawa couplings and irrespectively of the heavy spectrum. Also from an initially zero 0^3, 
this mixing angle could be as large as 30°. We have constrained the larger value of the Dirac 
neutrino Yukawa coupling to be perturbative (generally defined as < dir). 


Eq. (3.33) now defines the new mixing matrix as 


^pmns'”t^pmns 


u 


PMNS 


K, 


Ak, 


I U'' 

o'pmNS- 


Moreover, since the combination YyM~^Y^ is destroyed, the tree-level in¬ 
version given by the Casas-Ibarra formula makes the one-loop corrections 
sensitive to the arbitrary flavor structure of Y,^„ parametrized in three com¬ 
plex angles of TZ in Eq. (3.7). 

A similar relation was found for the non-supersymmetric contribution 
from the Higgs and Z boson [207, 208]. 

We look for deviations from a tribimaximal mixing pattern at the tree- 
level, especially 633 ^ 0, and figure out how much deviation from this 
very specific pattern is possible. The annoying part is the treatment of the 
a priori arbitrary structure of Yy. To be as generic as possible, we ran¬ 
domly scatter values for the three complex mixing angles in the complex 
orthogonal matrix TZ and consider for comparison two scenarios: one with 
roughly degenerate right-handed masses rtiy^^ within 

two orders of magnitudes (exact degenerate masses do not alter the mix¬ 
ing as discussed above) and second a large splitting with the three right- 
handed masses in the range 10^ ... lO^^GeV. Results are shown in Eig. 7. 
Since the Casas-Ibarra parametrization also needs the light masses as in¬ 
put (and always gives correct results), we fix = 0.3 eV. 

The presence of hierarchical right-handed masses therefore alters any 
preset mixing pattern. Without the need of flavor changing soft breaking 
terms, we, however, can still only achieve rather mild deviations from the 
initial 633 = 0. The other angles are allowed to scatter over the full range. 
The biggest deviation from 633 = 0 is also only possible for drastically 
large values of \TZij\. Since the angles of TZ are complex, its magnitudes 
of matrix elements are allowed to be larger than 1. This behavior gets 


The Casas-Ibarra 
parametrization 
relies on the 
tree-level 
formulation. 


Actually, we can 
assign any 
tree-level mixing in 
principle. However, 
(tri)bimaximal 
mixing is very close 
to the observed 
structures and can 
be motivated by a 
vast set of flavor 
symmetry models. 


Nevertheless, 9ij 
can be as large as 
30°. 








42 VARIANTS OF NEUTRINO FLAVOR PHYSICS 




Figure 8 ; Renormalization of the lepton mixing matrix by flavor changing self-energies at external legs 
according to the method by Denner and Sack [110]. 


also reflected into the definition of which anyhow challenges flavor 
model building (if one seeks an explanation of the strange values). To be 
conservative, we also constrain the resulting Yy to be perturbative and 
only allow solutions with max | Yj| < 4n. 


3.2.3 Renormalizing the Mixing Matrix 


The flavor diagonal 
terms in this 
procedure are 
infinite anyway as 
can be seen from 
below. 


In contrast to the 
quark mixing 
matrix, up and 
down sector are 
interchanged by 
convention. 


The mixing matrix of the charged current has a contribution from both 
up-type and down-type fermions as elaborated in Sec. 2.1.3 and 2.4. In a 
general basis (charged leptons not necessarily diagonal), the PMNS matrix 
is given by the combination 

l^PMNS = (^l) • 

Corrections to the mass matrices alter both and S^. The off-diagonal 
contributions can be absorbed by the following procedure into the mixing 
matrix, the diagonal contributions on the other hand can be absorbed into 
the wave function renormalization. 

We proceed with the diagrams of Fig. 8 and multiply the renormaliza¬ 
tion factors from the left (neutrino leg) and right (electron leg) being 
aware of the Hermitian conjugate mixing matrix in the vertex. 

= (l + Alfi) V'») (l + AUJ)’ 

[/(O) + (auIY , 

where [7^°^ denotes the unrenormalized, “bare” mixing matrix that is de¬ 
termined from tree-level flavor physics as in Sec. 3.2.1. Altogether, the 
renormalized interaction vertex with the W-Boson can then be written as 

^ ^T^Pl {u^°'>’^+Dl+dY , (3.35) 

where D^, are the correction matrices concerned with the left and right 
leg respectively. We have simply transferred the description from [89]. 
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The full contributions can be easily calculated attaching the generic self¬ 
energies Eq. (3.23) to the vertex and exploiting the equations of motion 
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The charged lepton contribution can be simplified and expanded in small 
mass ratios like m^/m^: 


AUl 
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1 Ue),LR 

1 y.(e) 
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-1 Ue)M 
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1 y.(e) 

171^^23 
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(3.37) 

^ y.(e),Ri 
31 

-1 y.{e),RL 

0 

J 



where we have neglected the chirality conserving self-energy contribu¬ 
tions which are suppressed by 1 / Mg^sY compared to the chirality flipping 
ones, see Refs. [89, 112]. 

In fact, this expansion does not work for neutrinos, especially not for 
the case in which all neutrinos are nearly of equal mass. There, we have 
to deal with the complete expression of AIJ^^\ Once again, we can safely 
neglect the part in Eq. (3.36a). Eor trivial tree-level mixing, = 5^j, 
we are left with 


^L,fi — 


AU 


(v) 






-l/i 


m: 


• m: 


(3.38) 


Eor the numerical analysis later on we shall keep all the contributions for 
convenience. Nevertheless, the simplifications of Eqs. (3.37) and (3.38) 
are enlightening to see where the dominant contributions come from. 

In principle, Eq. (3.34) shows two independent contributions, that are 
numerically quite different as well. The correction from the charged lep¬ 
ton leg suffers from a similar hierarchy in masses (with a result shown 
in Eq. (3.37)) as the renormalization of the quark mixing matrix does in 
Ref. [89] and lies numerically in the same ballpark. So AU^ may only 
account for minor corrections in the sub-percent regime. Nevertheless, 
the corrections from the neutrino leg are significantly enhanced due to 
the quasi-degenerate property of the neutrino mass spectrum. Only for a 
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pu{pi) = miu{pi) 

for Dirac spinors u. 


This means that 
corrections from the 
charged lepton leg 
only give CKM-like 
corrections and 
cannot account for 
the large leptonic 
mixing as also 
discussed in [112], 
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mOI/eV 

Figure 9; The neutrino mass spectrum gets more and more degenerate for increasing ^ as can be 
seen from the left side. The right side shows the increase of the enhancement factors fij, for 
better visibility f 2 s was multiplied by a factor of ten to separate it from /13. 


is a 

dimensionless 
quantity- 
parametrizing the 
loop correction and 
therefore 
suppressed at least 
by one loop-factor 
l/16n^ and the 
couplings inside the 
loop. 


Strongly hierarchical spectrum, the enhancement is gone similar to the hi¬ 
erarchical charged leptons. In this case, the large leptonic mixing follows 
from the method given in Chap. 6 and [205]. We see from eq. (3.38) an 
enhancement factor [209] 


ffi - 


my^rriy. 


(3.39) 


writing the self-energy = m^.crj.. The increase of fji with the neu¬ 

trino mass (and therefore the degree of degeneracy in the spectrum) is 
shown in Fig. 9. In this way, we can express the corrections in a model- 
independent way to figure out how large they have to be in order to give 
the right mixing and have 

O.40) 

iH 


with an arbitrary tree-level mixing matrix Estimates for aj. for no 
mixing in [7^°^ and tribimaximal mixing are shown in Fig. 10. 

Eq. (3.39) can be large in the case of a “large” neutrino mass as it 
might be measured by the KATRIN experiment in the near future. For an 
absolute neutrino mass around the discovery limit 0.35 eV and the mea¬ 
sured mass squared differences 7.50 x 10“^GeV^ and Amg^ 

2.46 X 10“^ GeV^ we have e. g. /21 1634 and /31 50. This property 

may account for a larger 1-2 mixing compared to the 1-3 mixing, when 
starting from the no-mixing hypothesis at tree-level. 

The determination of the three mixing angles obviously only needs to 
take into account the upper triangle: namely the absolute value of the 
t/i3 entry (sometimes called fixes sin 633 which can than be used to 
get sin 0^2 front U 12 and sin @23 from t/23- This procedure is well-known 
and was basically used during the years to precisely measure the mixing 
angles in absence of 1-3 mixing. 
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m^VeV m™/eV 

Figure 10; The model-independent analysis of the generic corrections as expected shows a decreas¬ 
ing behavior with increasing neutrino mass: the more degenerate the spectrum, the smaller 
corrections are needed. We show the result for = 1 left and for tribimaximal mixing at 
the tree-level right. 


To figure out what is possible by means of those corrections, we take 
a look at the case without tree-level flavor mixing and get the following 
(and pretty obvious) dependence on the corrections from the electron and 
neutrino leg: 


Ul2=\ 


i;.+ 
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c: 

<; 

1 _ 1 

Ul 3 =\ 



> 

c: 

<; 

1 _ 1 

U23 = 

[AUl] 


1 1 
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c: 
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(3.41a) 

(3.41b) 

(3.41c) 


If, by contrast, we consider a model with only no tree-level Oi^ (as in a 
tribimaximal mixing scenario), radiative corrections will lead to a non¬ 
vanishing 


Ci3 = [A(;t];2cS’+ [*!'!] 


'• u‘“’ 

I 13 ^33 


U 


( 0 ) 

11 




(3.42) 


even if the tree-level structure of no 1-3 mixing is preserved in the loop 
and there is no loop-induced mixing from the charged leptons. Take for 
this all = 0 in Eq. (3.42) and also = Oj then Ui^ is 

uniquely determined by the tree-level and the loop correction to the 
2-3 mixing 


The tribimaximal 


mixing is given by 


[210] 

^TBM 


Nl 

1 

0 ^ 
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Ve 
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1 
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V ^ 


V2 J 


U,, = u[l^[AUl]l^. (3.43) 

Therefore the size of ^3 is fixed to /Uj™ 0.15 x -/3 0.26 

in order to give the right L/33, which already is quite large for a typical loop 
correction due to the enhancement factor /32. The correction [^lfi]23 
is calculated in the underlying full theory which gives the self-energies. 
However, if the flavor structure is preserved in the loop, the size may 
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"We assume the 
“right-handed” 
sneutrino mass also 
a parameter of the 
SUSY scale, because 
SUSY breaking a 
priori has nothing 
to do with heavy- 
neutrinos. In any 
case, fh^ is 
negligible in the 
mass matrix if the 
right-handed 
Majorana mass is 
much heavier than 
a typical SUSY 
mass, 


We artificially set 

A® = 0. 

We take the 
following 
parameter ranges: 
Msusy e [0.5,5]TeV, 
Mj e [0.3,3] Mgygy, 
M2 e [0.5, 5] Mgygy, 
p e [“5,5] TeV and 
tan p e [5,60]. 


be estimated the following: ^3 ~ ^2™/(Idtr^) x /32 

0.22, where the factor estimates the loop suppression and/32 ^ /31 

50 was chosen for = 0.35 eV as described above. This estimate even 
indicates that the new physics in the loop needed to generate t/13 also 
may follow tribimaximal mixing. 

Corrections of this type have the power to completely mix up all imagin¬ 
able tree-level mixing patterns resulting in the observed PMNS structure. 
Especially, it is feasible to get the PMNS matrix out of a theory with no 
flavor mixing in the standard model fermions, where the flavor enters in 
the loops as one has the power of sufficiently large loop corrections. 


NEUTRINO MIXING MATRIX RENORMALIZATION IN THE MSSM Let US 
again discuss the effect of radiative neutrino mixing in extensions of the 
MSSM. The flavor changing self-energies are the same as previously dis¬ 
cussed and given in Eq. (3.24). Differently from the situation in Sec. 3.2.1, 
we now do not insist on degenerate neutrinos but in principle allow for 
an arbitrary mass spectrum before imposing the SUSY loop corrections. 
However, the spectrum is assumed to be sufficiently degenerate and the 
mass squared differences are the measured ones given in Eq. (2.36). 

Again, we want to figure out how much the flavor off-diagonal sneu¬ 
trino soft breaking terms can account for lepton mixing. We do not want 
to induce large flavor violation in the charged lepton sector, therefore we 
take the soft masses flavor-universal 


m2 = ml = fhl= 


Besides soft masses there are two more flavor matrices in the soft break¬ 
ing Lagrangian (talking about the vMSSM): the trilinear sneutrino-Higgs 
coupling A'' and the bilinear for the superpartners of right-handed neu¬ 
trinos. Literally, is the soft breaking counterpart of the right-handed 
Majorana mass in the superpotential—and therefore suspected to carry 
information from the high scale. As an order-of-magnitude estimate, we 
therefore follow the suggestions of [139, 140] and assign B^ = byM^, 
where by is supposed to be a parameter of the SUSY scale. 

We now allow for non-minimal flavor violating A- and B-terms and find 
values for the off-diagonal contributions in such a way that the flavor 
changing self-energies reproduce the weak current mixing matrix accord¬ 
ing to Eqs. (3.35) and (3.36), see also [206]. The charged lepton contribu¬ 
tion Eq. (3.37) is not only negligible compared to the enhanced neutrino 
one, we also suppress flavor changing soft breaking terms in order to avoid 
large charged leptonic ECNC. 

There are a few more parameters which we cannot and do not want to 
constrain: besides the gaugino masses and M2, the ;r-parameter of the 
superpotential as well as tan/3 and M^usy ^re arbitrarily and randomly 
chosen. We additionally scan over the lightest neutrino mass, and go 
as low as lOmeV. The renormalized masses, however, are not constrained 
to fit the corresponding physical masses for the purpose of Eig. 11. We 



3.2 RADIATIVE NEUTRINO MIXING 47 


lO"" 

10 ^ 

10^ 

10^ 

> 

cu 

52 10 -* 

3 

10^ 

10^ 

10^ 


10 ® 

10^ 

10 ® 

10 ® 



10 ® 

10 ® 

10 ® 

1 




500 1000 1500 2000 2500 3000 3500 4000 4500 5000 

^susy/^^^ 



500 1000 1500 2000 2500 3000 3500 4000 4500 5000 

^susy/^^^ 


0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

m™/eV 


Figure 11: The determination of soft SUSY breaking sneutrino parameters, Ajj and shows an 

interesting behavior: the two a priori independent parameter sets are very much the same! 
For the production of those plots, we varied several input parameters as the SUSY scale 
Msusy ^nd the lightest neutrino mass (which is also shown in the plots on the left 
and right respectively). The additional parameters which have been varied basically play no 
role in the ratio Ujj = Ay/MgysY ^^d bij = b'^j/Msusy ^hat reproduce the neutrino mixing 
matrix. 


want to study the flavor mixing contribution of the self-energies. Interest¬ 
ingly, the results for and are identical! Why that? 

The reason behind this unsettling observation lies on one hand in the 
choice of the neutrino parameters, on the other hand it is a physical re¬ 
sult of the flavor structure of the self-energies. To understand Fig. (11), 
we confess our input choice: the neutrino Yukawa coupling was chosen to 
unity (y’^ = 1) and the scale of right-handed masses appropriately to have 
the right order of magnitude for the effective light neutrino mass. This 
choice has degenerate neutrino masses with mass m^\ Because individ¬ 
ual contributions for different flavors do not mix (i. e. ~ Al,bp,the 
flavor-diagonal self-energies to not alter the result for the flavor-changing 
ones and can therefore be adjusted in such a way that the renormalized 
masses equal the physical ones. 

We can see from the structure of the dominant diagrams shown in The heavy- 

propagator o/Vr 

counts as -rjy • 

Ml 
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Figure 12: Self-energies renormalizing the neutrino Majorana mass with contributing (27^, left) and 
(27^, middle). The diagram on the outer right (S^) shows a suppressed contribution of 
The cross shows the AL = 2 insertion. 


Fig. 12 that the scaling of both A- and B-term contribution is the same: 
The A-term diagram scales as 

27^~A"^MRy"=yX/MR, 

Mr 

and the B-term contribution 


2j ~ 


(y^)TMR^B2 —MrY^ = 


M" 






/M, 


R> 


We also made use 
of the replacement 


with = Yyl and Mr = MrI. As can be seen from the third diagram of 
Fig. 12, there are also contributions that are suppressed with more inverse 
powers of Mr, e. g. 


\iC 

Zj ~ 


Ml Ml ^mI 


J V V 

M^ 



Ml • 


The dominant behavior already showed up in the effective mass matrices 
of light sneutrino states, Eq. (2.51a), the AL = 2 contribution exactly has 
one (symmetrized) term from the LR transition. 


m° Mf 


Mi 




^ DlVt 


M, 


R 


and one with the “soft SUSY breaking Majorana mass” B^, 

Ml + (m^)^] Bl {Ml + Ml) MRm°^ 


ml Mp 


v^y^B^ 

u-ty y 

Mi 


In a way, we are free to choose either A^ or B^ to generate neutrino mixing 
radiatively. In scenarios where there is a connection between A"^ and A®, 
such as GUT-inspired models and LR-symmetry, large sneutrino A-terms 
lead to large FCNC contributions for charged leptons. In those cases, the 
sneutrino B-term may be the appropriate choice for radiative lepton flavor 
violation. 
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Figure 13: Dyson resummation of the neutrino propagator at the W vertex. 


3.2.4 Remarks on a Resummation of Large Contributions 


The realization of quasi-degenerate neutrino mass eigenstates—as would 
be the case for the mass squared differences Am^ being much smaller 
than the light neutrino mass: Am^ <C m ^^—may lead to arbitrarily large 
corrections due to Eq. (3.39) shown in Fig. 9. Those large contributions 
in general have to be resummed. To clarify the meaning of those radia¬ 
tive corrections we perform a Dyson resummation of self-energies on the 
neutrino propagator, shown in Fig. 13. However, unlike the usual case, 
we have to deal with flavor non-diagonal self-energy matrices that carry 
chirality structures as well. 

The general formalism of propagator dressing for inter-generation mix¬ 
ing, Majorana and unstable fermions was already explored by Refs. [211- 
215], in the following we review those results of Dyson resummed prop¬ 
agators and apply them to our formulation of mixing matrix renormaliza¬ 
tion by the means of SUSY corrections as described in Sec. 3.2.3. 

As in Eq. (3.23), we decompose the self-energy matrix into its Lorentz- 
co variants: 


“Usual” means a 
situation like in 
QED: the resummed 
propagator is just 
given by the 

geometrical series 
i 

p—m—S ' 


s = p {slP, + + SlP, + S^Pi,, 

where we have omitted flavor indices. To keep track of the proper count¬ 
ing, we remind us that the “scalar” components carry mass dimension one 
whereas the “vectorial” part has dimension zero. 

Performing a Dyson resummation for the most generic propagators, we 
have to pay attention for both Dirac and flavor matrices. As usual, we 
write the resummed propagator as 


(3.44) 


sIPl - S^Pr) - (m + S^P, + S^Pr) ’ 
which can be decomposed in the following way: 

iS(P) = ^{ [P (1 - + ("> + )] ^’^ 


iS(p)= 


— m — E 


The notation is 
intuitive: —-—= 

p—m—I] 

means 

[pi - m - r] . We 
have to keep in 
mind, that E is not 
only a matrix in 
Dirac space, but 
also in flavor space. 


(3.45) 
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w^here the denominator D(p^} can be obtained as 

Dip^)=[p-m-p - {slP, + ] X 

[p+m- (siP, + slPil) + (i:f Pr + S^^Pj,) ] (3.46) 

=p"{^- (i - ) - (^++^3 • 

In the case of Majorana fermions, we have (using 27^ = S^* = S^ and 

27^ = ^ 27^) 

D{p) =p^\l-S^f'-\m + S^\^ . (3.47) 


The notation might 
be a bit confusing: 
while E are 
self-energy 
components, E are 
the components of 
the dressed 
propagator. 


However, for general self-energy matrices the situation gets worse, since 
we have to deal with matrix algebra. The result for the dressed propagator 
is given by [211, 214]: 

So- ip)=P Pr + P ) 0 Pl + Pk + (^f) .. Pl, (3.48) 

and the coefficients are obtained by the expressions 


sl = {A,)-\p^-Cj,C,)-\ (3.49a) 

si = {A^rHp^ - (3.49b) 

si = {Ai,)-^B,{A,)-\p^-Cj,C,)-\ (3.49c) 

Sl = {A,)-^BM~Hp^ - (3.49d) 


where 


CR — Bf,{Ajl) Ci—Bi{Ai) 


(3.50) 


and 





f 


(Al.r)o = 5ij - ( 27 ,%) .., (Bi,K)o = m°5,j + (27f J .. (3.51) 

The inverse of a matrix can be computed using the adjoint adj(M): 

1 


M“^ = 


det(M) 


adj(M), 


such that 

(p2 -M)“^ 


adj(p^ -M) 
det(p^ — M) ’ 


with IVl — CrCi^, Cj^Cr. 

Eqs. (3.48), (3.49), (3.50) and (3.51) define the full propagator, in¬ 
cluding all flavor mixing effects from the self-energies. In order to redo 
the calculation which resulted in Eq. (3.36a), we have to recall and clarify 
the meaning of the mixing matrix renormalization with external legs. The 
external-leg contribution to the full matrix element of the flavor transition 
is given by 


(y) 


-^kg=“(P/)^/; 


p-mj 


(3.52) 
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w^here w^e only include flavor changing contributions with j ^ / to the 
mixing matrix counterterms AUf of Eq. (3.34). The flavor diagonal parts 
vanish in the field renormalization of the external particle (and also renor¬ 
malize the mass). In this description, the external leg is only partially am¬ 
putated. Let us remind the procedure to arrive from the n-point Green’s 
function at the S-matrix elements via truncation of external field lines. The 
truncated Green’s function are given through multiplication with inverse 
propagators, see e. g. [216, and original references]: 

G„(pi,...,p„) = G2^(pi,-pi)---G2^(Pn,-p„)G„(pi,...,p„). (3.53) 


We then have the truncated 2-point function including the self-energy cor¬ 
rection 


--27 


(x) 


— m 


/I 


f} rh — 


m,- 


t J trunc 


= -nif) --27 

^ 10 — nif \ - 


W_ [ 

"fj n - 


m, 


-mj), 


which is G 2 = ■ In Eq. (3.52), we can only truncate “half” of the ex¬ 

ternal lines, since one is internally connected to the vertex. Multiplication 
from the right with the inverse propagator, — nij:'^, and appending the 
wave function u(p^) yields the expression (3.52), modulo the LSZ factor. 

The same can be done using instead the dressed propagator of Eq. (3.44) 
decomposed in (3.48). We have infinitely many self-energy insertions as 
in Eig. 13 with n —> 00 and Eq. (3.52) turns to 

= ^(Pf) (1^ - m - I7W) (3.54) 

with m and matrices in flavor space. The flavor diagonal part goes 
into the field renormalization and redefines the masses as 

Pf {Pf 27 }}^) - - -^ = Uf {pf ), 

p-ruf - 27 ^^ 

which is a generalization of the wave function renormalization conditions 
of [212] or rather [217]. This is only correct as approximation, if the 
off-diagonal self-energy terms are small; in general, the inverse matrix 
depends on those. We also keep the diagonal flavor dependence in the 
full propagator, so this condition is not to be seen strict. Eor the flavor 
changing part, a 27^^^ factor remains (with u{pf)p = u{jpf)mf') and we 
have 

which, inserted into the full matrix element for the flavor changing transi¬ 
tion in the W-vertex, gives 


We do not explicitly 
write the field 
renormalization 
constants because 
the flavor changing 
corrections are 
finite. 


We identify the 
mass eigenstates as 
those corresponding 
to the renormalized 
pole mass. 


tiyp,) [Ai)<’'>]^, uM = 

MPf) [r^)..u^{pi), 
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Figure 14: We have inverted the formulae for the mixing matrix renormalization to find the proper 
values of off-diagonal -terms. The same values of Ajj leading to the right amount of mixing 

in the formulation of Eq. (3.38) starting from trivial mixing give negligible contributions 
once the flavor-changing self-energies are resummed (labeled as 1/res)- The enhancement 
factor is not present anymore, once diagrams are resummed. On the other hand, if one 
takes the resummed expression to determine the off-diagonal -terms, the same values in 
Eq. (3.38) give—as expected—diverging results (1/ren)- 



tviththe (unrenormalized) vertex , and the “dressed” 

^ (v) 

counterterm matrix AUl turns out to be 


AU. 


(v)' 


J fi fj 


+ m^ 




+ m^ 




(3.56) 




tvhere we made use of the Majorana-specific relations = 27^* and 
27^ = Sf* and readopted the obvious notation of Eqs. (3.36a) and (3.38). 
The dressed propagator components are given in Eqs. (3.49) and 
summation over j is implied of course in Eq. (3.56). A similar relation 
holds for the charged lepton contribution—^which on may resum as well, 
though it is not parametrically enhanced as the neutrino leg. 

The interesting results including the Dyson resummation of the neu¬ 
trino propagator compared to the unresummed formulation are displayed 
in Eig. 14. We adjust the diagonal A-terms according to the mass renor¬ 
malization of Eq. (3.22) in order to fix the diagonal neutrino masses. The 
off-diagonals are then chosen to generate the PMNS off-diagonals and 
therewith the leptonic mixing angles. As can be seen from Eig. 14, the 
resummed counterterms j^AL/^'] .. immediately drop down where the un¬ 
resummed corrections get enhanced by Amj-. Eor the purpose of 

Eig. 14, we fed the A-terms responsible for the appropriate PMNS mixing 
into the either resummed or unresummed formulae. The plot on the left 
side shows proper mixing if the unresummed counterterms are used, the 
plot on the right side has proper mixing for the resummed version. On 
the other hand, using the resummation, the unresummed counterterms 
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explode when the enhancement gets sizable (note that they are paramet¬ 
rically enhanced by large A-terms anyway). 

Note, that the depend on the momentum squared, which is for the 
neutrino contribution the mass squared of the outgoing neutrino (in this 
case, we cannot setp^ = 0, see Eqs.(3.49) and (3.50), (3.51)). Moreover, 
it is worth noting that the dressed propagator components carry negative 
mass dimension, = — 2 and = — 1, which is crucial for the 

counterterm (3.56) being dimensionless and the full dressed propagator 
of dimension minus one. 

Let us perform a simple and rough order of magnitude estimate: The 
only mass scale that appears outside the self-energy loops is the scale 
of light neutrinos. Especially the dimensionful self-energies should not 
exceed the tree-level masses. Moreover, all quantities related to loops are 
suppressed by a factor of ~ 10“^. Estimating the values of Eq. (3.56), we 
find that the mass scale drops out and an overall loop-factor is left, so 
the suppression might be ~ 10“^! The leading terms of the E are not 
loop-suppressed. However, the situation for the dressed self-energies is 
not so clear as for the undressed ones. Looking at Eqs. (3.49) and the 
definitions of the matrices A, B, and C shows that the loop-suppression 
factor is subdominant, since the self-energies only enter as corrections to 
order-one parameters in A and B and the overall suppression factor is 
again the general loop-factor. 

Although the possible enhancement of Eq. (3.38) is absent in Eq. (3.56), 
we add up more contributions including the diagonal ones by summing 
over j = 1,..., Hg numbers of neutrinos. Therefore, we also collect flavor 
non-diagonal contributions from sources, that were not available with the 
method described in the previous section. Especially, there can be result¬ 
ing effects to mixing angles that are not there at tree-level—take for exam¬ 
ple a theory with extended flavor structures as the vMSSM. In such theo¬ 
ries, one can in principle introduce additional flavor violating terms like 
the trilinear A-terms—however, by dressing the propagators and there¬ 
fore the mixing matrix, it is not anymore a one-to-one correspondence 
that each A^j generates a corresponding mixing matrix element Uij (with 
t and j the same) and we get contributions to one element from the others. 
Since the counterterm elements in Eq. (3.56) are of the form EE, where 
E denotes the dressed self-energies, one easily sees that for example the 
1-3 element gets a contribution 


AU 


(y) 


J 13 


~ EijEj3 — EiiEi3 


E^ 9 E‘ 


12^23 


E^ o E' 


13^33- 


Even if 27^3 as well as Ei^ vanishes — there is the term from (1,2) x (2,3) 
left: Ei 2^23 ^nd the third mixing angle gets a finite renormalization by 
the other two mixing contributions. 

The interesting behavior of the Dyson resummed corrections to the neu¬ 
trino mixing matrix is displayed in Eig. 15. Here, we notice two things: 
first, in the case without resummation the strength of the corrections 
(~ AT / MgusY ^nd ~ i)T / Mgysy) decreases with a more and more de- 
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Figure 15: We compare the Dyson resummed (left panel) and not resummed (right panel) solutions 
for ^4^^-terms (red crosses) and -terms (blue stars). In the first row, we also varied the 
scale of the neutrino Yukawa coupling from ~ 10“®... 1. Where only blue stars are seen are 
the points similar to Fig. 11 where and -terms are the same. This only holds for the 
y.v = 1 case because of the different scaling with the Yukawa coupling of the contributions 
as discussed in 3.2.2. 


Still, we assume for 
simplicity 


generate neutrino mass spectrum. This observation is not striking but 
rather what we expected, because of the increasing enhancement factor 
fij = For the resummed contributions there is no depen¬ 

dence on the neutrino mass—also as expected, because the mass depen¬ 
dent factor was absorbed in the resummation. Second, if we do not restrict 
the overall neutrino Yukawa coupling to jy = 1, but let it flow in some 
regime (for this display, we have chosen jy e [10“^, 1]), the needed cor¬ 
rection strength can be significantly reduced into the percent regime or 
even lower. 

One can argue whether it is reasonable to follow the previously de¬ 
scribed procedure in order to cope with the enhancement which results 
from the quasi-degenerate property of the neutrino mass spectrum. On 
one hand, we are dealing with the flavor off-diagonal effects of the self¬ 
energies as small effects, renormalizing the masses only via the diagonal 
contributions. On the other hand, the self-energies themselves may still be 
perfectly perturbative and the large contributions are rather due to the en- 
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hancement factor of Eq. (3.39) leading to large entries of in Eq. (3.40). 
Ignoring the charged lepton contribution, we have 

= (l+A[/f). (3.57) 

Summing up similar contributions, it is probably better to include the full 
A [7^ and not only the propagator part, so 

[/r^^ = [/(0) (^1 +A[/f+ (Al/f)%...) (l-A[7f)”\ (3.58) 

The matrix AUj has only off-diagonal elements in the description accord¬ 
ing to Sec. 3.2.3 and the inverse is easily calculated with |^A[7]^'j.^. = 
fijCJij = u^j (we take the self-energies as real and remind the reader of 
the symmetric property of Majorana self-energies) 


1 

^12 

^13 

1 

1 “h 

1 -r ^23 

^12 ~ ^13^23 ^13 “T ^12^23 

~^12 

1 

^23 

detu 


1+^13 U23-Ul2til3 


“^23 

1 j 


1 

1+^12 J 


(3.59) 


with det u = 1 + + ^13 + ^ 23 - interesting to note, that with = 

u^j/ (1 +ufj) at most mixing matrix elements of 0.5 can be generated 
from trivial tree-level mixing only! Additional contributions that interplay 
in Eq. (3.59) tend to diminish the results further. We conclude, that for 
large ~ AT the proper mixing matrix elements get severely reduced 
and due to the fact that there now exists an upper bound on t/[“, we 
can never generate the observed neutrino mixing in this procedure. In the 
resummed propagator approach, we still have been able to find values of 
AT which reproduce the PMNS mixing. 

To conclude this digression on our attempts to resum the previously 
discussed large contributions to neutrino mixing, we state that neither of 
the described ways shall be applied at all. Eurthermore, the perturbative 
approach to the mixing matrix renormalization which is valid for quark 
mixing (and has been successfully applied in the MSSM [89]) is actually 
wrong if large mixings or quasi-degenerate masses are considered. The 
approach relies on the neglect of (Au[)^ terms which formally can be 
included as in Eq. (3.58). However, the mass and mixing renormalization 
cannot be disentangled anymore once there are substantial mixing effects. 
Instead, the inverse propagator of Eq. (3.44) defines actually a new mass 
matrix. Diagonalization of this matrix leads to a proper renormalization 
of mass eigenvalues and mixing matrix elements. This procedure has al¬ 
ready been addressed in Sec. 3.2.1 and now has been shovra to also be 
the proper procedure for non-degenerate masses. However, a not exactly 
degenerate mass spectrum at tree-level already intrinsically obeys a mix¬ 
ing pattern. We have argued in Sec. 3.2.1 that degenerate masses are 
somehow more likely to produce such a quasi-degenerate physical mass 
spectrum radiatively. In case the masses are rather hierarchical, radiative 
corrections are small anyway. 
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SUMMARY OF CHAPTER 3 In this chapter, we investigated the influence 
of radiative corrections to the origin of neutrino mixing. The large mixing 
angles in the lepton sector compared to the quark CKM matrix suggest 
a tree-level mechanism at work that generates this large mixing. How¬ 
ever, we showed that especially in the case of degenerate neutrinos at 
tree-level where there is no mixing at all, one-loop threshold corrections 
have the power not only to lift the degeneracy but also to reproduce the 
large mixing angles. A supersymmetric theory has generic flavor violation 
in soft SUSY breaking. If we impose an MSSM extended by right-handed 
Majorana neutrinos, we can generate neutrino mixing stemming from soft 
breaking trilinear couplings. Without degenerate masses, there has to be 
a tree-level mixing at work, which again can be trumped by radiative cor¬ 
rections. The SUSY threshold corrections to the leptonic mixing matrix are 
enhanced by the quasi-degenerate pattern of neutrino mass eigenstates— 
as long as their masses are not strongly hierarchical. Additionally, it was 
shown that hierarchical right-handed Majorana masses also enter the ra¬ 
diative corrections and alter the tree-level mixing pattern even without 
soft SUSY breaking. In any case, taking (SUSY) one-loop corrections into 
account, the proposed tree-level flavor structure whatever it may be gets 
significantly changed. The most powerful corrections to do so are, of 
course, flavor changing threshold corrections. Those allow to connect the 
observed flavor with SUSY breaking in a supersymmetric context. This 
way, we may solve the flavor puzzle by solving the SUSY breaking puzzle. 


THE FATE OF THE SCALAR POTENTIAL 



The following introduction follows basically standard textbook descriptions — 
see e.g. [17, 18, 218]. To get a comprehensive overview, it is nevertheless 
necessary to carefully take a look into those (or other) standard references. 
It is not the purpose of this thesis to review textbook knowledge — however, it 
is somehow compulsory and convenient to have a short survey over the field 
before getting started into the phenomenology of effective potentials. 

We observe in nature the spontaneous breakdotvn of a fundamental 
symmetry. The breaking of electroweak symmetry directly reflects the fact, 
that W and Z bosons are massive. In the SM, this property is added by 
hand by virtue of a scalar field having a potential whose minimum does 
not respect the original symmetry anymore. That way, we have sponta¬ 
neous symmetry breaking (SSB) at the cost of introducing a new field and 
adjusting the model parameters (i. e. the /r-term of the Higgs potential) 
manually. The theory does not give SSB itself. Does it? 

In a seminal paper, Coleman and Weinberg (CW) showed that SSB can 
(and does) occur, once quantum corrections to the classical potential are 
switched on [219]. They started with a symmetric theory, calculated the 
effective potential explicitly and pointed out, that the one-loop effective 
potential itself leads to spontaneous symmetry breaking. 

At roughly the same time, Jackiw derived the effective potential directly 
from the path integral and generalized the formulation to multiple scalar 
fields [220]. Lee and Sciaccaluga came from a different direction and 
showed, that the one-loop effective potential can be obtained from inte¬ 
grating one-loop tadpole diagrams [221]. This method does not only pro¬ 
vide an efficient and straightforward way to calculate effective potentials 
diagrammatically doing the resummation of all diagrams via the integral, 
but also glimpses towards something called the renormalization group im¬ 
proved (RGI) effective potential. Improving the effective potential means 
stabilizing it with respect to the choice of the renormalization scale. It 
turns out, that not only the RGI potential is scale invariant, but also any 
of its derivatives—and thereby the derived n-point functions. 

We briefly recapitulate the ground-breaking achievements of Goleman 
and Weinberg, Jackiw, as well as Lee and Sciaccaluga, and discuss the 
meaning of the RGI effective potential. Later on, we apply the machinery 
to the MSSM, reproduce the GW result and discuss its implications on 
the stability of the electroweak vacuum in the presence of light stops and 
sbottoms (where light means TeV-ish). 


The main purpose 
of the following 
introduction is to 
set the language 
and notation. 
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Note that in held 
theory, the 
Lagrangian C is the 
Lagrangian 
density—to end up 
with “L = T-V” 
we have to 
integrate over 
space: L = ^ d? xC. 


0 is commonly 
called “classical 
held” because it 
behaves like a 
classical held value 
rather than an 
operator. 


In general, not 
every possible value 
of ^ can be 
achieved with 
translationally 
invariant J. There 
are some 
“forbidden” regions 
following from 
n on -perturba tive 
effects [223], 


For simplicity, we 
suppress here 
consistently the 
label 2 PI in the 
notation of the IPI 
potential. 


4.1 THE EFFECTIVE POTENTIAL AND ITS MEANING FOR THE GROUND 
STATE 


In classical field theory, the meaning of the field theoretical potential V {(j)) 
of a field cj) is analogous to the potential energy of a particle in classi¬ 
cal mechanics: V{(j)) denotes the potential energy density of that field <p. 
Things have to be defined properly in a quantum field theory where the 
classical field (j) turns into a quantum field. 

In a quantum world, the “potential” turns to be the generating functional 
for one-particle irreducible (IPI) Green’s functions and can be derived in 
the path integral formalism as was shown by Jackiw [220]. For clarity, we 
distinguish between the IPI-potential Vipj(0) and the effective potential 
I4ff(^) of a field (j) according to Dannenberg [222]. In that view, 
has the same meaning as the classical field theoretical potential—where 
the field (j) = 4’{x) has to be replaced by a constant value (p{x) = 
which can be seen as the vacuum expectation value 


^ = 


{0\(t>\0) 


( 0 | 0 ) 


(4.1) 


in presence of external sources J. In this way, the ground state of the 
theory is determined by the value of 0 in the limit of vanishing sources, 
V = 01 value, the potential V of a real scalar field in 

= ^ (5^0)^-V(0) 


is minimized: 


dV 
d 0 


- 0 . 


(f) = (p=V 


(4.2) 


The classical potential I4i is identical to the potential written in the La¬ 
grangian density—in quantum field theory language it is the “tree-lever 
potential. For a real scalar field this could be (including only renormaliz- 
able terms obeying a global Z 2 reflection symmetry 0 —> — 0) 


^Cl = ^(<^) — m^0^ + A0"^. 


In the case < 0, the minimum of the potential does not respect the 
original symmetry anymore, so it is spontaneously broken. The IPI poten¬ 
tial now is the sum of all IPI diagrams that contribute to V as a function 
of 0. On the classical level, “tree-level”, “effective” and “IPI” potential are 
equivalent and given by the same expression. Sticking to renormalizable 
operators, this potential can at most be a quartic polynomial. The loop ex¬ 
pansion also adds IPI diagrams with more externals legs to the potential. 
In general, those diagrams include an arbitrary number of external legs: 
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w^here is the IPI n-point Green’s function (evaluated at vanishing 
external momenta, p = 0). Looking at Eq. (4.3), it is obvious that this 
potential is the generating function for all IPI n-point diagrams which can 
be obtained by differentiating with respect to </>. 

The IPI (effective)^ potential as given in Eq. (4.3) interpreted as the 
non-derivative part of the Lagrangian density going beyond tree-level is 
a powerful tool to survey the structure of SSB. Treating the one-loop po¬ 
tential as a function of its minimal value determined via Eq. (4.2) de¬ 
scribes the ground state of the quantum theory—and is therefore a better 
approximation to the proper ground state. 

THE LOOP EXPANSION Summing up all one-loop diagrams with an ar¬ 
bitrary number of external legs also goes beyond fixed-order perturbation 
theory: each vertex which adds one (in a 0^) or two (cj)^ theory) external 
legs comes along with one power of the coupling. The n-point contribu¬ 
tion is actually 0(A"). Nevertheless, it is a feasible way to use an on first 
sight artificial expansion in an O (1) parameter a which is defined via the 
condition 

£((^,5^(^,a) = a“^r((/),5^(^). 

The important point here is that a multiplies the complete Lagrangian as it 
is, irrespective of shifts and scaling in the fields. A power-series expansion 
in a is equivalent to an expansion in the number of loops: each vertex 
obviously brings one factor whereas every internal line or propagator 
of a diagram contributes a —the propagator is the inverse of the quadratic 
terms, that also come along with a“^. The power P of a which multiplies 
each diagram can be easily calculated as the difference of internal lines I 
and vertices V {P = I — V) and the number of loops is given by 

L = P + l= 7-y + l, 

because tree-level diagrams shall take as the total Lagrangian. 

The loop expansion still lacks knowledge of higher order terms (namely 
two-loop contributions) if we stop the expansion at first power. Moreover, 
non-perturbative effects like bound state formation that also lead to spon¬ 
taneous symmetry breaking [225, 226] cannot be included in that formal¬ 
ism. 

DISTINCTION BETWEEN EFFECTIVE AND IPI POTENTIAL Eor Complete¬ 
ness, we also briefly review the discussion about the so-called “convexity 
problem” which was clarified by Dannenberg [222] and comprehensively 
reviewed in the appendix of Sher’s Physics Report [227]. The problem, in 
fact, does not exist. However, connected to the question of convexity is the 
existence and interpretation of an imaginary part of the effective potential 
addressed by Weinberg and Wu [228]. 

1 The language used for this potential is not unique throughout the literature. Since more 
commonly effective instead of IPI is used, we will adiabatically change notation. 
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The subtle 
difference between 
generating function 
and functional will 
be pointed out in 
Sec. 4.1.1 about the 
path integral. 

The easy thing with 
the effective 
potential is, that it 
is an ordinary 
function of the 
classical value f 
and minimization 
can be done 
without pain. 


The argumentation 
basically follows 
[219] whereas the 
original idea behind 
the loop expansion 
is hidden in [224]. 

Actually, a can be 
identified with the 
Planckian constant 
h which also 
multiplies the full 
Lagrangian and 
exactly equals one, 
ft = 1, in natural 
units. 



Tree diagrams have 
a~^ because they 
only include one 
vertex or two 
vertices and one 
internal line. 
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Figure 16: The convex envelope for the non-convex potential with non-degenerate minima is given 
by the dashed interpolation left. In case both minima are degenerate, the convex effective 
potential is flat in between the extrema (right). 


The state |n) is the 
vacuum state in 
presence of external 
sources J: 

\n)'-^\o). 


The decay rate 
which is supposed 
to be calculable 
from the imaginary 
part has to be seen 
for the decay from a 
spatially averaged 
(i. e. constant) Held 
value which is 
localized in Held 
space (like the 
minimum of some 
non-convex 
potential) into a 
spatially localized 
conHguration like a 
bubble. This is only 
valid during the 
phase transition. 
This decay rate is 
not to be confused 
with the decay of 
the false vacuum 
itself [231, 232]. 


The rigorous distinction between IPI and effective potential follows 
from the path integral and shall be elucidated in Sec. 4.1.1, whereas 
the pictorial interpretation shown in Fig. 16 can be understood quantum- 
mechanically: we define the the effective potential V'eff(0) as the expecta¬ 
tion value of the energy density (n|.^|n) where (p = (0|^(x)|Q). This 
is the real and convex potential for the spatially constant field value <^, 
see [228] and [229, 230]. The potential which is field theoretically calcu¬ 
lable, however, is a different though related object. For the sake of clarity, 
we call this explicitly If we now consider a potential as drawn 

in Fig. 16 on the right-hand side, where there are two minima at two 
distinct values in field space (say ^ = ±cr), the energy density for field 
values |(^| < cr in between the two minima is minimized by a superpo¬ 
sition of the two vacuum configurations ^ = a and 0 = — cr and the 
average energy density is given by y(cr) = y(—cr), which is the dashed 
line in Fig. 16, right. For most purposes in quantum field theory, however, 
we are interested in configurations where the state is localized in field 
space (a homogeneous state) but may be averaged (constant) in space- 
time, = (^. Such localized configurations are unstable because the 

operation of field localization does not commute with the Hamiltonian 
and therefore they decay, where the decay rate can be estimated from the 
imaginary part of the (IPI) potential [228]. Configurations like this are 
the case if we are interested in perturbative excitations around a local 
minimum—as done in electroweak theory where we expand around the 
minimum. A similar configuration may happen where we have to expand 
around the local minimum of Fig. 16 on the left: it may have happened 
that the universe while cooling dovra ended up in the higher minimum. 
The transition to the deeper minimum is then a phase transition and the 
relevant potential is given by V^pj rather than Vgff, which corresponds to 
the dashed curve and only has the global minimum. Using Vgff, no phase 
transition could happen. Instead, V^pi can be taken as the analytic contin¬ 
uation of 14ff) see [227] and [233, 234]. 

The occurrence of an imaginary part in the effective potential is related 
to the non-convexity of the classical potential, which can be easily seen 
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rewriting the CW potential that will be introduced in Eq. (4.29) as done 
in [228] 

Vi-ioop(0) = + ^ [V"((^)2lny"(0)] +P((/)), (4.4) 

where V"{<p) is the second derivative (i. e. the mass term) of the classical 
(tree-level) potential. Now, develops an imaginary part when the 

argument of the logarithm drops negative, thus V"{(j)) < 0 which means 
V is non-convex where V^.ioop is complex. 

An illustrative example of a decaying configuration was given in [236] : 
let us distribute electric charges in a way that in a given (and arbitrarily 
large) space-time region there is a constant electric field. In QED, the con¬ 
stant electric field will pair-produce electrons and positrons from the vac¬ 
uum with a non-vanishing probability which shield the external charges 
and destroy the constant field. Such a configuration obviously is unstable 
which follows from the imaginary part of the analytic continuation. 

VACUUM DECAY The semi classical description of a false vacuum decay¬ 
ing into the true vacuum configuration was first given by Coleman him¬ 
self [237] from which this brief review heavily draws. Quantum correc¬ 
tions to the semiclassical theory were then considered in a follow-up pa¬ 
per [238] which is essential to calculate transition rates. 

In many field theories as extensions of the SM with more scalars. Grand 
Unified Theories, supersymmetric theories, but even the SM itself false 
vacua as shown in Eig. 16 may appear. The question whether we live in 
a false or the true vacuum is of cosmological importance: an infinitely 
old universe inevitably ends up in the true vacuum (global minimum). 
However, we know that our universe is not infinitely old and we know the 
history of hot big bang: as the universe cooled down, we may have ended 
up in the false vacuum. 

A situation like this is drawn in Eig. 17, where we assume to have ended 
up for some reason in the classically stable equilibrium state 0 = 

The state of lowest energy, however, is given hy cj) = (j)_. This state can 
classically only be reached by thermodynamic fluctuations—or via barrier 
penetration in a quantum world. Eor simplicity, we only consider a single 
real scalar field with a given potential energy density U 

(4.5) 

The decay probability F of the false vacuum per unit volume V follows 
the exponentially suppressed solution for barrier penetration 

r/V =Ae-’^^^[l + 0{h)], (4.6) 

where the quantum theory will give corrections 0{h) and the factor V fol¬ 
lows from the fact that the emerging solutions (named “bounce solution”) 
are not translation invariant and all spatial translations have to be inte¬ 
grated over. As pictorial view, Coleman states that for a decaying vacuum. 


The problem only 
occurs, if a wrong 
(unstable) vacuum 
state was used to 
build up the 
effective potential 
as was pointed out 
in [235], 
Nevertheless, the 
resolution of the 
problem still lies in 
the use of 
vocabulary. 

The term “false 
vacuum” describes 
exactly what it is: 
the true vacuum is 
expected to be the 
global minimum of 
the potential energy 
(density) whereas a 
local minimum 
mimics the vacuum. 
By knowledge of 
quantum mechanics, 
classically stable 
expansions around 
that minimum will 
decay with certain 
probability to the 
deeper ground 
state. 
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Figure 17: Instable potentials: on the left side, we have a false vacuum configuration at <p^ where 
the true vacuum lies at (j)_. The potential on the right side shows a quantum mechanical 
analogue, where the local minimum at qo is classically stable. 


the product of F/V with the four-volume of the past light cone has to 
be of order one to observe this decay. The coefficients A and B have to be 
calculated in the underlying theory which, in general, is not even possible. 

The presence of a bounce solution can be seen from the quantum me¬ 
chanical tunneling process as an allegory for the barrier penetration of 
the vacuum state. Consider the quantum mechanics of a single particle in 
one space dimension with a potential drawn in Fig. 17 on the right side. 
Classically, the equilibrium at qg is stable forever, whereas the quantum 
particle can penetrate the wall until it escapes at q = cr and propagates 
freely for q> a. The analogue of Eq. (4.6) is given by 

r=Ae“®/'*[l + 0(ft)], (4.7) 


where the division by unit volume is a division by one and B can be calcu¬ 
lated using the WKB method 


B = 2 


ra 

dqVW. 


■''Jo 


(4.8) 


d dL dj^ 

dt dq dq 


The potential 
energy in Fig. 17 
was adjusted in 
such a way that the 
total energy E = 0. 

Traditionally, 
5 jdq^/2(E-V)=0 
is considered, 
which is Eq. (4.9) 
with fixed 
integration limits, 
E = 0andV ^ -V. 
This actually causes 
the transition to 
imaginary time. 


A generalization in more than one dimensions was given by [239] but 
makes no difference for the qualitative discussion (details in [237]). The 
interesting observation Coleman made, is that, using the Euler-Lagrange 
equation for the one-particle Lagrangian L = |q^ — V (q). 


dV 



and energy conservation |q^ + V = E, q = const. In order to minimize B, 
we have to consider the variation 


dqVW = 0, 


•'40 


(4.9) 


whose solutions are determined by the differential equation 

d^q dV 
dT dq ’ 


(4.10) 
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w^here we have substituted the time by imaginary time t = i t. The energy 
conservation now turns to 

l(^ \ -v=o 


2 ydT 
and Hamilton’s principle 

drlg = 0, 


with the “Euclidean” version of the Lagrangian, = +^- Choos¬ 

ing the time where the particle “reaches” cr to be t = 0, the classical equi¬ 
librium q = qo is only reached in the limit t —> —oo and the “velocity” 
vanishes at t = 0: dq/dT = 0. The solution for B from Eq. (4.8) turns 
then to be 

~q rO 

dqV2V = dTLg. 

J Qq fc' — 00 

The equation of motion is invariant under t — > — t, so the particle bounces 
at T = 0 and goes back to qg as t —> +00. The coefficient B of the “bounce 
solution” is just given by the Euclidean action 


These requirements 
can be identified 
with boundary- 
conditions for the 
equation of motion 


(4.10), ^ 
and 
liniT 

go- 


= 0 


»±00 


T =0 


q(T) = 


B = 



= Se. 


(4.11) 


The difficulty in the generalization to a quantum field theory lies in the 
calculation of the Euclidean action for the bounce solution. To find the 
bounce, take the equations of motion for the scalar Lagrangian (4.5) in 
imaginary time description 




ci> = u'{4>) 


(4.12) 


and impose the bounce boundary conditions 


lim 0 (t,x) = </)+ 

T—>±00 


(4.13a) 


and 

d d) 

-^{0,x) = 0. (4.13b) 

OT 

Because of SO (4) invariance under Euclidean rotations, condition (4.13a) 
can be identified with 

lim (j){T,x) = </)_!_, 

|j?|—>00 

which gives the pictorial description of a spatially growing bubble of the 
decaying vacuum, where the false vacuum stays unchanged in regions far 
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Practically, one 
would obtain 4>{p) 
as numerical 
solution for this 
differential 
equation and 
integrates 
Eq. (4.14) 
numerically to get 
B. 


A has to have mass 
dimension four and 
is a dimensionful 
parameter related 
to the scale of the 
problem, so the 
expectation is a 
rough estimate of 
feelings. 


We may at least 
refer to the 
comfortable 
situation of 
observers observing 
the universe today 
as necessary 
condition for a 
stable vacuum. 


away from a bubble induced by quantum fluctuations. The final result can 
be simplified using SO(4) invariance and writing + |xp = p^, then 
Eq. (4.11) turns to 


/* 00 


B = 2n^ 


dp 


Idcj) 
2 dp 


UW 


(4.14) 


with boundary conditions 


d (f) 
dp 


= 0 , 

p=0 


and 


^(oo) = (^+. 


The equation of motion (4.12) determines <p{p)'. 

d^ (j) 3 dcj) dU 
dp^ p dp dtp 

What is left is the determination of the coefficient A in front of the expo¬ 
nent. However, the precise value of A does not play a significant role for 
the estimate of the decay rate. Therefore, only rough numerical estimates 
exist, where the exact analytical calculation [238] yields 


I det'[—d^+ 

^^ 4^V det[-d^ + U"{(j>^)]' 

The bounce solution is given by cp, and is the position of the false 
vacuum. The determinant is genetically difficult to evaluate (where det^ 
means the determinant without zero eigenvalues). Therefore, the value 
of A is usually estimated by the height and width of the barrier, because 
one expects A to be of the fourth power of the scale related to the barrier. 
The uncertainty of A can then be estimated to be ~ which can be 

related to an effective uncertainty of B [227]. Sher takes as estimate on 
A the barrier height^. For all practical purposes, a standardized height is 
used related to the electroweak scale andAf« (100 GeV)^ [240-242]. Sta¬ 
ble configurations are then determined hy S^[p]/h > 400 for a universe 
lifetime of 10^° years. 


OTHER IMPLICATIONS OF FALSE VACUA The question about stability, 
instability or metastability of any ground state, especially the electroweak 
vacuum we are believing to live in, is of major importance for the existence 
of observers in a current situation. Call this anthropic principle, in any case 
the fact that observers observe a certain broken phase which seems to be 
somehow stable makes us believe that it has to be like this. Previously, we 
discussed estimates on the vacuum life time and completely masked out its 
meaning. If we find configurations for the effective potential of a theory 
which predicts extremely short lifetimes of false vacua, we may argue 
that in view of a 13.8 billion year old universe any state with a lifetime 
much less than this number for sure has been decayed. Configurations 
of charge and color breaking false or true vacua as may occur in various 
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SUSY models, where the drastically enhanced scalar sector also contains 
colored scalars, will be postponed to the MSSM section 4.2. 

One interesting question to deal with is the influence of gravitational ef¬ 
fects: imagine the transition of a false into a true vacuum releasing a vast 
amount of energy. This energy portion is expected to grow with the vol¬ 
ume of the newly formed bubble. The bubble grows as the gain in volume 
energy supersedes the surface energy. Without gravity, a growing bubble 
can always be generated as long as it appears to be large enough. Includ¬ 
ing gravity into the game, we have to deal with the Schwarzschild radius 
of the energy density inside the bubble which also grows with the energy. 
Coleman and De Luccia [243] have shown (and proven) with very clear 
arguments that gravitation plays no significant role—as far as it was un¬ 
derstood at that time (with no new insights in the modern literature). The 
obvious thing is that conservation of energy is not violated by the tunnel¬ 
ing process and in the end the total energy of the bubble (inner negative 
volume energy and positive surface energy) has to be exactly zero. There¬ 
fore the gravitational held cannot do anything. In general, gravity has the 
tendency to prevent false vacua from decay and reduces the decay proba¬ 
bility. Finally, Coleman and De Luccia discuss the possibility of living in the 
leftovers from an early decay of a false vacuum where the world outside 
the bubble corresponds to the true vacuum. In this special scenario (decay 
into zero cosmological constant) the probability for decay gets enhanced 
by gravitational effects. 

The existence of several vacua (even in the pure SM there is a second 
minimum in the effective potential popping up around the Planck scale) 
allows also for degenerate vacua and thus coexisting phases. From the 
demand of such coexisting phases, Froggatt and Nielsen [244] predicted 
about 16 years prior to the discovery the SM Higgs mass to be 135 ± 9 GeV 
which just at the edge of the errors covers the recently discovered Higgs 
boson [15, 16] and they significantly reduced the uncertainty of the (at 
that time recently measured) top mass and shifted the central value down 
to 173 GeV (from 180 ± 12 GeV) with 5 GeV uncertainty. The foundation 
of this prediction is the “multiple point principle” (MPP), first applied to a 
coexistence of gauge couplings at the Planck scale [247, 248].^ Similar to 
coexisting phases in condensed matter physics (the presence of ice, water 
and vapor at the triple point of water), specific quantities take only very 
specific values after the surroundings are fixed. The vacuum stability issue 
was used to give lower bounds on the SM Higgs mass [249-263]. 

4.1.1 Path Integral 


Correlated to the 
total energy density 
is the question of 
the cosmological 
constant. In a pure 
quantum Held 
theory on fixed 
Minkowski space 
time, arbitrary 
constants can be 
added and 
subtracted to any 
potential energy 
density such that 
there is no absolute 
zero energy. Two 
non-degenerate 
vacua make a 
cosmological 
constant either 
appear or disappear 
after vacuum decay, 
adjusting the 
parameter in the 
either true or false 
vacuum we choose 
to live in. 


There are many 
papers predicting 
many values of the 
Higgs mass, so 
Ref [244] is for 
sure not the only 
one being right 
within the errors. 
The MPP also 
allows to form a 
dark matter 
candidate from top 
bound states with a 
reduced top 
mass [245, 246]. 


The most stringent derivation of the effective potential without explicitly 
calculating Feynman diagrams follow from the path integral definition. 


2 H. B. Nielsen answered the question about the reason behind the MPP that there are no 
deeper reasons behind most principles therefore they are principles [private communica¬ 
tion, PASCOS 2014]. 


Jackiw calls the 
Feynman 
diagrammatic 
calculation “an 
onerous task” which 
it indeed is, as 
described in the 
section about sea 
urchin diagrams. 
Sec. 4.1.3. 
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Already the two 
loop case gets much 
more involved. 


Z[J] gives all 
Green’s functions. 
The relevant ones 
for scattering 
amplitudes, 
however, are only 
the connected ones. 
Also for the 
definition of the 
potential, which 
gives the 
(self-) interaction, 
only connected 
Green’s functions 
are important. 
(Disconnected 
Green’s functions 
do not describe 
interactions.) 


In general, 

J d4x[-Veff((^) + 

... ], which for 
(p = const, reduces 
toEq. (4.19). 

Vgff((p) is now an 
ordinary function of 

4>- 


Jackiw calculated the effective potential for n self-interacting scalar fields 
up to two-loop precision [220]. The expansion of the path integral allows 
to keep track of the loop-expansion, as each term comes along with a 
factor h counting the number of loops. 

The effective potential is defined via the effective action which follows 
from the generating functional W[J(x)] of connected Green’s functions 


OO / ■\n r 


i W[J] = ^ d'^Xi • • • (xi,...,x„)J(xi) • • • J(xj. 

n=0 d 

(4.16) 

The “partition functional” Z[j] gives the unconnected Green’s functions 
Z[j] — I (4.17) 


and the vacuum transition amplitude from t = —oo to t —> +00 in pres¬ 
ence of sources J(x) 


Z[j] = (0_oo|0+oo)j- 


The classical action is defined via the Lagrangian of the theory S [^] = 
J d^x£(^) and the generic field ^ may represent all fields in that theory. 

The effective action is the generating functional of IPI Green’s functions 
and follows from a Legendre transform 


r,ff[<P{x)]=w[J]- 


d^ xJ(x)^(x). 


(4.18) 


The classical field value (p minimizes Fgff with respect to J, 


- 5W[J] 

(b = -- such that 

^ 5J(x) 


5r,(f[<P] 


54>{x) 




-J{x). 


Spontaneous symmetry breaking occurs, if p 0 even for vanishing 
source J(x): 5r^ff/= 0, and Ggff is still minimized. The effective 
potential is a space-time independent quantity following from Fgff after 
division by the 4-volume 


ref{[4>] = 


d^x. 


J 


(4.19) 


The effective potential is the generating function (not functional) of 
connected IPI diagrams where the n-point vertex functions G^"^ (for zero 
external momenta) follow from the n-th derivative and > 0: 

OO 2 

= = (4.20) 

n=l 

In this way, V^ff is defined via the generating functional of IPI Green’s 
function and therefore, using this language, Vgff and I4pi are identical. To 
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resolve the convexity problem which was an issue in the early literature 
of effective potentials, we refer to Dannenberg [222] and his dictionary. 
The crucial point is, that is the double Legendre transform of 
which can be seen from Eq. (4.18) using W[J] as generating functional of 
connected Green’s functions, where 


w[j] = ripi[^(x)]- 


d^xJ {x)<P(x) 


J 


(4.21) 


itself is the Legendre transform of F jpj being the generating functional of 
IPI diagrams. Since Legendre transforms are convex functions, the double 
Legendre transform of T^pi is the convex envelope of this function (al). 

There is one further issue related to the imaginary part of the effec¬ 
tive potential, which is connected to a non-convex tree-level potential and 
shall be picked up later where appropriate. 


PATH INTEGRAL DERIVATION OF THE EFFECTIVE POTENTIAL lu the 
following, we briefly state the result of Jackiw’s path integral calculation 
of the effective potential [220] without going into details. In Ref. [220] 
not only the formal proof of the path integral derivation is given, also the 
first two-loop result for the effective potential of self-interacting scalars 
and a recalculation of the CW model (radiative breaking of massless, 
scalar QED) in arbitrary gauge was performed. The latter calculation al¬ 
ready raises the discussion about gauge dependence of the effective po¬ 
tential and its physical interpretation, see for a modern discussion and 
resolution of that problem [264-266]. 

The main argument uses the existence of the “classical” (x-independent) 
field value ^ which has the meaning of a vev. Erom the classical action 
S[(j)] = f d^xC[(j){x)], we obtain the shifted action S[(p + (p] which for¬ 
mally can be expanded as Taylor series around p = 0 

Sip + p] = s[<^] +S^[p]p + -S"[(^]</)^ + +..., 


where the prime denotes derivative (variation) with respect to p, S'[p] = 
. Then 


^S[P + P] 


4>=Q 


ai-Gffi'i’] - 


= Vpe 
'iPi 




(4.22) 


Vpe 




IPI 




where the integral is performed over all IPI configurations and £"{p) = 
— V"{p) with the tree-level potential V{p). 

The shift defines a new Lagrangian [220] via 


d^x^ (x) 


J 


5S[P] 

5p(x) 


<t>=4’ 


xF{p,p{x)). (4.23) 


J 
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According to 
Dannenberg, “the 
resolution of the 
convexity problem 
lies not in physics 
but in vocabulary”. 
Having this in mind, 
we use the terms 
“effective” and “IPI” 
synonymously 
throughout this 
thesis and never 
refer to the “true” 
effective potential 
(which is convex). 


The gauge 
dependence issue is 
not relevant for the 
purpose followed in 
the following, 
because we do not 
calculate the gauge 
part anyway. 

We thank Luminita 
Mihaila for pointing 
this and Eq. (4.22) 
out in her institute’s 
seminar talk in 
November 2014. 

This allows to nicely 
find Eq. (4.26) from 

= 

r exp(iS[<»]) ^ 

S[<^] + 

I lndet£"((i)). The 
effective potential 
has an additional 
minus sign. 

C!' in general is a 
matrix in the space 
of multiple scalar 
fields 4> = 4>a with 
a = 1,... ,n. 


S[p + p]-S[p]- 
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In the spirit of 
Jackiw [220], we 
explicitly write 
factors ofh to make 
obvious that the 
loop expansion is 
an expansion in h. 


In the h-expansion, 
isO{h) and the 
first term in the 
expansion of the 
exponential 
vanishes in between 
the angular 
brackets. 


Note that this 
causes the 
cosmological 
constant problem! 


The new Lagrangian £ comprises interaction terms with (^-dependent cou¬ 
plings. In this way, we decompose £ into a propagator and an interaction 
part according to [220] : 


d'^x£(0,^(x)) = 


d xd j-0^(x)iD^J(^;x,y)^f,(j) 


+ 


(4.24) 


d^x£^{(t),4>(x)), 


where we generalized the notation to an arbitrary number of scalars, 
a,b = l,...,n, and suppressed indices where possible. The propagator 
is to be calculated in an obvious manner 


^Dj{<p]x,y) = 


5^S[ct>] 


54>a{x)5(l)b{y) 


(4.25) 


4'=<t> 


Transformation into momentum space follows straightforwardly 




d'^xe'^^iD^^^(</);x,0), 


such that the effective potential can be obtained as 


Veff(0) =^o(^) - 


d^k 


{2n) 


Indet 


+ i(,(exp - 


(4.26) 


d'^x£i((^;(^(x)) 


where V'o(</)) denotes the tree contribution and the determinant operates 
on all indices a, b counting for internal and spin degrees of freedom. The 
last terms means calculation of the vev in between the angular brackets 
and averaging over space-time and is 0{h^) contributing to the two-loop 
potential. The one-loop contribution to the effective potential can then be 
evaluated (as an example for a single-scalar theory) 


/ -X 

Vi{4>) = -j 


d^k 


In [k^-m^{^)] 


with the field-dependent mass m{(p) = V^'{(j)). Introducing a cut-off A, 
the integration results in [220] 




6471" 


[yl^ln 1 + 






(4.27) 


-m"^(^)ln 1 + 


A^ 


m2(0)_ 


■A^m^((p)] 


up to an omitted overall constant A^. In dimensional regularization, the 
D-dimensional integral can be performed after Wick rotation to Euclidean 
space [264] (where D = 4 — 2£') 


--Q 
2 


2s 


^ d^ k r 2 1 

In [-k^ + m^] 


3 

In ^- 

Q2 2 


(27i)^ -^4 (471)^ 

with the renormalization scale Q and A^, = - — -f In 471. 


(4.28) 
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+ 


Figure 18: The one-loop effective potential for scalar 0^ theory. 


4.1.2 The “Coleman-Weinberg”potential 


What today is knotvn as CW potential is the result of a Feynman diagram¬ 
matic calculation of the one-loop effective potential. The result is identical 
to Jackiw’s path integral derivation and the integrated tadpole of Lee and 
Sciaccaluga. Moreover, CW shotved that the calculation is applicable not 
only to scalar 0^ theory (sample diagrams shotvn in Fig. 18), but also to 
scalar QED, non-abelian gauge theories and theories with fermions. And 
the basic formula is very simple and easy to use: 

Vi-ioop = +P;(0)] , (4.29) 


where the sum runs over all fields in the loop and the trace over the mass 
matrices Mi - The coefficients C; count the number of spin degrees of free¬ 
dom (and are negative for fermionic states). In the mass matrices, the 
vev is replaced by the classical field 0 and Q is the renormalization scale. 
In the definition of the effective potential, there is one arbitrariness left 
due to the choice of the renormalization scheme which expresses itself in 
certain polynomials Pi(0). 

Similar to the diagrammatic method described in the following section 
(which was used to explicitly derive the one-loop potential with third gen¬ 
eration squarks), CW sum over all IPI n-point Green’s functions and per¬ 
form the loop-integration after the summation. In Sec. 4.1.3, we do it the 
other way round and sum after integration. 

There are two issues to treat with care as CW point out in their pa¬ 
per [219]: one is the correct interplay of combinatorial factors (symmetry 
factors of the diagrams and statistical factors like a Bose | factor). Second, 
in the derivation there appear superficially infrared divergent diagrams 
that result in singularities at 0 = 0 for an initially massless theory. There 
is one confusing paragraph in [219] stating this infrared singularity be¬ 
comes obvious if one calculates radiative corrections to the propagator at 
= 0 which would behave like p^ Inp^ (which goes to zero as —> 0). 

Actually, the effective potential itself is calculated at vanishing external 
momenta to give the vacuum configuration of the theory. According to 
CW, this infrared singularity can be avoided staying away from 0 = 0. Is 
the effective potential invalid at the origin? 

The infrared behavior can be directly seen in the one-loop potential 


A 1 1 


d^k 


00 


E 

n=l 


1 

2n 


+ 


IS 


, (4.30a) 


The trace can be 
taken as the sum 
over the mass 
eigenvalues. 

In the MS scheme 
for a scalar field in 
the loop, 
e.g. P(0) = 
whereas for gauge 
fields P{f) = —|. 


The tree-level 
Lagrangian without 
counterterms for 
massless theory 
is given by 
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with renormalization constants B and C. After summation and Wick rota¬ 
tion to Euclidean space 


A 1 1 1 

4! 2 ^ 4! 2 




In 1 




. (4.30b) 


Now, performing the integral and introducing a momentum cut off A, CW 
obtain 


A 1 1 AA^ 

4r 2 ^ 4! ^ 64712 


2 j;4 




x^4) 


In 


A</)^ 


2567i2 I 2A2 


2 .t4.31) 


and after applying renormalization conditions to remove the constants 
B and C, the final result for the tree plus one-loop effective potential in 
scalar theory is given by 


^ -4 


A204 

256712 



(4.32) 


with the renormalization scale Q entering via the renormalization condi¬ 
tions. 


Details about the 
proper scale choice 
follow in the 
description of the 
RGI potential. 


Ref [268] is 
interesting in such a 
way that the 
authors derive the 
analogous 
expression to the 
scalar effective 
potential with 
Yang-Mills fields as 
external (and 
internal) particles. 


The small ip 
divergence 
corresponds to CW’s 
IR singularity and is 
absent with masses. 


d^y 
d ^2 


= 0 


(/)=0 


and 


d^y 


= A(Q). 

P=Q 


The renormalization scale Q is a priori an arbitrary number with mass 
dimension one—and has to be chosen in a proper way. This choice may 
be a simplification of the calculation (e. g. Q = ^ to make the logarithm 
vanish [267]) or a specific fixed scale choice in order not to make all the 
corrections vanish but intrinsically make them small in a neighborhood 
where the corrections shall vanish, short: Q = v such that for classical 
field values ^ = v at the minimum the logarithm vanishes. 

Similar considerations can be done and have been done for calculating 
effective scalar potentials in e. g. massless scalar QED and the electroweak 
SM [219] as well as the effective potential for a Yang-Mills field [268]. In 
any case, once gauge fields are added to the theory and the derivation 
of the effective potential, this object becomes explicitly gauge dependent, 
whereas physical observables that can be derived from it (e. g. tunneling 
rates, not the instability scale [264]) stay gauge invariant, especially the 
value of the potential at its minima [265, 266]. The issue of gauge in¬ 
variance is not relevant if one only considers scalar theories or, as done 
in Sec. 4.2, a scalar subset which gives the dominant contribution to the 
loop-corrected effective potential. 

The effective potential in the given description is a good approximation 
where ^ is close to the minimum determined by Eq. (4.2). However, for 
field values ^ » v, where the logarithm ln(0/v) becomes large, the ex¬ 
pansion becomes unreliable (the same holds for very small field values 
0 0). In order to avoid this, the effective potential has to be improved 

via the RG as already briefly stated in the introduction to this chapter 
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and discussed in greater detail in Sec. 4.1.4. What stays a bit unclear is 
the question what defines a large logarithm. This becomes clear, if largely 
separated scales are discussed as one compares the electroweak with the 
Planck scale [253, 256, 258, 259, 261, 262, 269-283] or the scale of 
heavy neutrinos as done in Chap. 5 and [274, 284-286]. To get a feel¬ 
ing: field values 3v as discussed in Sec. 4.2 produce ln3 1, which 
is not a large number (therefore we do not include the RGI potential in 
Sec. 4.2). On the other hand, ln((^/v) 47i for ^ 3 x lO^v indeed is 

a large number. This comparison is to justify the use of the RGI potential 
for the description of instabilities caused by heavy neutrino fields and not 
for SUSY scale instabilities. 


INTEGRATING THE TADPOLE The diagrammatic result of CW can be 
very elegantly obtained by the method of Lee and Sciaccaluga [221]. Es¬ 
pecially the evaluation of the combinatorial factors in the diagrammatic 
calculation can be very exhausting as shovra in the following section. The 
tadpole method exploits the utilization of the effective potential as gen¬ 
erating function of n-point functions. The loop-expansion gives the cor¬ 
responding loop counting for the n-point function. In this way, the first 
derivative of the one-loop potential gives the one-loop tadpole 

=G^^^ = Mcj)). (4.33) 

o<p 


Inverting Eq. (4.33) gives V^ff after integration of the tadpole with respect 
to 0. Eor a single scalar field with self-coupling potential P{4>), the tad¬ 
pole can be easily calculated in dimensional regularization yielding 




32712 




Q2 


-M2(0) 


(4.34) 


where the functions / and M are defined via derivatives of the potential 


M2(0) 


dcj)^ 


and/(</)) 

4>=f 


d^P{(j)) 

dcf)^ 


(t>^f 


d<p 


Therefore, the integration gives including a generic prefactor N 




N 


32712 




N 


327i2 

N 

647i2 




In 


Q2 

m2(0) 3' 

Q2 2 


(4.35) 


N already has to be 
present in the 
equation for T. 


The prefactor N collects all spin, color and charge degrees of freedom. A 
squark e. g. gets N = 3, where the corresponding quark has N = —3 ■ 2. 
Eq. (4.35) once again reproduces the generic GW formula. 
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+ 
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+ 





Figure 19: Summing up the one-loop potential for a triple scalar coupling up to infinitely many external 
legs. We refer to diagrams with an arbitrary number of external legs as sea urchin diagrams. 


All functions and 
series are 
well-behaved, so we 
can interchange 
sum and integral. 

Up to now, no 
sparticles have been 
found at the LHC, 
so squarks have to 
be heavy if they 
exist. Nevertheless, 
heavy stops are 
desired to get a 
reasonably heavy 
Higgs. Therefore, it 
is not unlikely to 
have multi-TeV 
squarks, one Higgs 
below 130 GeV and 
heavier Higgses in 
between. Such 
scenarios also allow 
to consider an 
effective 2HDM 
where SUSY 
corrections 
generate some of 
the generic Higgs 
couplings. Details 
on that were 
already given in 
Sec. 2.2.2 and 
follow even more in 
Sec. 4.2. 


4.1.3 Sea Urchins 


Coleman and Weinberg [219] performed their calculation of the effective 
potential by directly resumming all one-loop IPI diagrams, and calculated 
the loop integral after summation. The tadpole integration method by Lee 
and Sciaccaluga [221] includes the sum over all diagrams implicitly by 
the integral over the one-loop tadpole and Jackiw’s path integral deriva¬ 
tion [220] also implicitly includes all n-point diagrams in the determinant 
and the integration over the loop momentum follows after the implicit 
sum. Formally, summation and integration can only be exchanged with 
care. We shall, however, show that it makes no difference for the calcula¬ 
tion of the effective potential and explicitly resum the series of analytical 
n-point Green’s functions. The interesting point is, that we can trace back 
the loop-integral to derivatives of the tadpole functions. For the sake of 
generality (and because of later applicability to SUSY models), we con¬ 
sider triple and quadruple scalar couplings and restrict ourselves to the 
contribution from one class of heavy fields in the loop. This is the case 
when we calculate the Higgs effective potential by means of heavy squark 
fields in the loop and do not consider the effective quark potential, which 
is impossible to calculate (and analyze) analytically. 

The effective potential is the generating function of n-point diagrams. 
Sample diagrams are shovra in Fig. 19, where we only show triple scalar 
couplings (we also shall start our discussion with those). Corrections to 
renormalizable operators would stop immediately after the first diagram, 
which gives a correction to a certain self-coupling A; if the external lines 
are related to Higgs fields, say and/or H^. The loop expansion allows 
to attach more external lines, even an infinite number of such spikes. We 
refer to the spiky diagrams as sea urchin diagrams for obvious reasons. 

Each term in the expansion by the number of external legs defines an 
effective self-coupling, which can be obtained backwards from the deriva¬ 
tive of the effective potential. This allows to cross-check the methods. As 
we will see, the direct summation of sea urchins results in exactly the 
same structure as what follows immediately from the compact result given 
many times, cf. Eqs. (4.4), (4.29), (4.28), (4.35). Therefore, we define 


;t(4) 2 


2 

00 ^(2n) 


(4.36) 
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explicitly for a complex scalar field 4>. The bilinear term ~ is to be 
seen as a tree-level mass which gets renormalized by a UV counterterm. 
The self-couplings only receive finite corrections and in the following, we 
are only interested in the one-loop self-couplings. In principle, there is 
a which, however, is only an additive constant to the loop-induced 
quartic coupling. All others are exclusively loop-generated. We 

get the coupling of the 2n-point function by differentiating with respect 
to <P ' : 







4>t_4,=0 


(4.37) 


To be clear on factors n! and n and some more appearing in the following: 
the 1/(n!)^ in Eq. (4.37) eats up factors nl coming from the derivative and 
the n in the numerator is just the factor 1/n from the definition of self¬ 
couplings of the potential in Eq. (4.36). Likewise, we can calculate 
as one-loop diagrams and match the Green’s functions to the effective 
potential. Exploiting this fact, the loop calculation yields an expression 
like 


-iA 


(2n) 

loop 


. {n'-f r2n 


= 1 ■ 


1671^ 


iVL„(mi,m2), 


(4.38) 


where )' 2 n denotes the 2n-point vertex, some loop-2n-point-function 
and nil, m 2 ^^e the masses of two different fields running in the loop. N 
counts for example the numbers of color. 

Marrying (4.37) with (4.38): to perform the matching, we have to iden¬ 
tify the loop potential contribution with the result from the effective poten¬ 
tial. Given in Eq. (4.38) is the truncated Green’s function missing external 
fields. The effective potential follows by appending those fields $ and 
as 


00 2 
n=0 V / 

Now, Eqs. (4.39) and (4.36) describe the same object in the full or the 
effective theory, respectively. The mass term was omitted. 

The potential can be expressed in terms of the loop-derived quantity: 



N 

1671^ 


00 


^r^"L„(mi,m2), 

n=2 


(4.40) 


which will be breathed with life in the following. 

We perform a sample calculation in a scalar theory which can be mo¬ 
tivated from the MSSM. Gonsider a superpotential of two Higgs doublets 
and two chiral quark superfields 




(4.41) 


The triple coupling induces the top Yukawa coupling, which is important 
for the fermionic contribution to the effective potential, and a quadrilinear 
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The field takes 
the role of the 
“classical field” as it 
appears externally. 


$ has to be replaced 
by the Higgs field 
h° later on. 


The occurrence of 
the factor (n!)^/n 
in Eq. (4.38) can be 
seen from the 
symmetries of the 
diagram: the 
internal lines of the 
diagrams are 
different at each 
vertex. Contracting 
all vertices to a 
connected diagram 
then gives a 
symmetry factor 
n!(n- 1)!. By 
drawing only one 
diagram, we forget 
about the 
others—so this 
factor ends up in 
the numerator. 

There are no further 
symmetry factors, 
so we receive a 
prefactor of 


The superheld 
formalism was 
introduced in 
Sec. 2.2.2. Higgs 
fields and Q are 
SU(2) doublet 
fields. For the 
moment, we ignore 
the bottom Yukawa 
coupling. 
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Definitions of the 
loop functions Aq 
and Dq are given in 
App. A. There are 
also some (obvious) 
technicalities how 
to arrive at the 
derivative for equal 
masses. 


stop-Higgs coupling ~ RtL,R^u*^u> which will be discussed later. Leaving 
out soft SUSY breaking terms (the full case will also be discussed later 
and in Sec. 4.2), the only trilinear stop-Higgs coupling is due to the F- 
term contribution and couples to the “wrong” Higgs doublet, ~ /rY*t^tL^d' 
This vertex is quite unique since it is oriented (by the direction of complex 
fields) and couples three distinct fields. We calculate the one loop four 
point function (i.e. the one loop correction to the tree-level A^) for equal 
left and right stop masses. We work in the flavor basis, where the internal 
lines are no mass eigenstates and are propagators with the soft breaking 
left and right stop masses and /tir (iV^ = 3 is the number of colors) 


-iA^ =iW 


\j^ 

16n^ 


diL, dtp 




, liuTtp 1 f d Y ^ ^ 


(4.42) 


which can be easily generalized to the 2n point function resulting in A 


{2n). 
loop ■ 


loop c 1671^ n (2n —1)! ydm^y 


2n-l 


Ao(m), (4.43) 


where 


—1 

dm^y 


2n-l 


Mm) = (-1) 


, (2n-3)\ 

= -iM : ... 


M) 


2\2n-2 ■ 


Therefore, we can write the effective self-couplings in the following form: 


A(2") = - 


Nc 

16n^ 


fxYt 




2n 


{mY 


{nY 


(2n-3)l 

M-i)! ■ 


(2-6n-L4n2)-i 


(4.44) 


Now, we can combine Eq. (4.44) and the loop potential of Eq. (4.39) re¬ 
sulting in an infinite series that can be resummed to 


loop 


N^m 


4 00 




^2n 


In 


16712 ^2- 

n=2 


6n + 4n2 


(4.45) 


At the end, we are 
interested in the 
neutral 
component’s 
direction ~ 
Because of SU(2) 
invariance, we also 
get the full 
potential from the 
calculation of the 
part. 


where we defined the dimensionless quantities 



The sum of Eq. (4.45) can be evaluated yielding 

1 A n , , . Ox , 

V >4 =- {3x^ - 4x Artanh(x) - (1 +x2) ln(l -x^)) 

1671 ^ 2 ^ 2 

= (3^2 - (1 +x)2ln(l +x) - (1 -x)2ln(l -x)) , 

(4.46) 
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w^here the last replacement Artanh = |(ln(l+x)—ln(l—x)) can be 
done for lx| < 1 which is also the radius of convergence of the infinite 
sum (4.45). Likewise, for |x| > 1, log(l — x) develops an imaginary part 
which corresponds to the values of {/i°) for which one stop mass-squared 
eigenvalue drops negative—an obvious relation to a non-convex poten¬ 
tial at tree-level, in this case it is the stop potential. The imaginary part, 
however, is puzzling since it rises with x^ ~ On the other hand, the 
lightest stop is tachyonic as x > 1 and its mass^ drops with x. Given the 
superpotential of Eq. (4.41) and only including the /i*Yt-term and the soft 
breaking masses, the field dependent stop squared mass matrix is 


m?M) 


ml ^ 


and the two eigenvalues setting /til = /tir = m. 


2 = ± x). 


Such a situation would be excluded anyway (loop corrections to the stop 
masses may rescue the tachyonic mass in some regime) and is either a sign 
of a wrong theory or some missing ingredient. There is indeed something 
missing which is actually an important contribution to the stability of the 
potential, the |Yt/i°|^-term. Now, we could argue that this does not matter 
because it gives a contribution to a different direction in field space (~ /i° 
not /t°). To be prepared for full stop contribution to the Higgs effective po¬ 
tential want to be complete and also add the missing soft 

SUSY breaking trilinear coupling ~ The full (ignoring potential 

stop vevs, (u) = (Ir) — 0) stop mass matrix is then given by 




mr 


\Y,h' 


0|2 


and the eigenvalues by (/til = m^ = m) 
mi 2 = m^(l±x + y), 

with 


x2 = 


Whl-ii*Y,h\ 


0 * 


and 




-li*Y,hl*+A,hl 


y = 






(4.47) 


(4.48) 


(4.49) 


Using the generic formula (4.29), it is very easy to obtain the one-loop 
potential induced by stop fields in the loop: 




_N,m^ 
“ 32712 


[(l + x + y)2ln(l+x + y) + (l-x + y)2ln(l-x + y) 

-(x2+y2 + 2y)(3-21n(mVQ^)], 


The main question 
is, whether one can 
trust the expansion 
beyond the radius 
of convergence. For 
the discussion of 
vacuum stability, we 
take the analytic 
continuation where 
the potential 
becomes complex 
and ignore both 
issues related to the 
imaginary part and 
radius of 
convergence. 
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(4.50) 
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Figure 20: The series of n-point Green’s functions contributing to the Higgs effective potential with top 
squarks running in the loop. 
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where we kept the renormalization scale dependence (because we can). 
Eq. (4.50) follows directly from Eq. (4.46) with ±x —> ±x + y modulo 
the polynomial in y. The direct resummation of the series of IPI n-point 
Green’s functions, however, seems to be impossible for the generic setup, 
where the stop fields do not only couple to one Higgs field but also to 
a linear combination of and h^, h = /i°* —/t^Af/ (M*Tt) which appear 
mixed in the series shovra in Eig. 20. Note that the l/i^p coupling preserves 
“chirality”, where the trilinear coupling ~ h flips L —> R and vice versa. 

The summation includes mixed powers of |/t°|^ and \h\^ and has the 
form 

00 00 

Vi-ioop = -Y,J]^kn{h*h)HKXr. (4.51) 

fc=0n=0 


The three diagrams 
give the 
contribution to the 
same order in the 
expansion 
Eq. (4.51), fc = 1 
and n = 2. 


The coefficients can be calculated from the loop diagrams with 2k + 2n 
legs, where k counts the number of vertices (and the conjugated 

ones) and n counts the quadrilinear tRL-tR,L“^u*‘^u vertices. 

The calculation of the coefficient a^.Q has already been performed above, 
the symmetry factor is found to be effectively 1/fc. A similar result can 
be obtained for from the pure quadrilinear coupling diagrams. The 
dependence on the masses in the loop is in both cases a certain derivative 
of the tadpole function Ao(m^). Eor the mixed diagrams like 



\ / 

y 

4 4 

--►-H 



W 4 

, .K A- 

' 4 5 ' 

/ \ 


the combinatorics is more involved. We start with the trilinear coupling 
diagrams from above (Eig. 19) which have fc!(fc — 1)! internal symmetries. 
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Divided by (fc!)^ from the 2fc-th derivative of the effective potential, we 
get 1/fc. Now, we attach quadrilinear vertices to the fc-point trilinear cou¬ 
pling diagram. Remember that the quadrilinear coupling preserves chiral¬ 
ity, the trilinear one flips it. Therefore, the n-th quadrilinear coupling just 
“sits” on top of any propagator and we only have to count how many pos¬ 
sibilities are there to distribute one such coupling on an existing diagram. 
The same problem arises if one has a box of colored gummy bears having 
k colors (the existing k propagators) and one takes n bears out of the box. 
How many combinations are possible? The answer is 

(n + fc-1)! _ fn + k-l\ 
in-l)\k\ ^ k )' 


We can now give the explicit and generic coefficients aj.„ 




for fc > 1, 


aon = l^tP"- 
n 


2fc+2n^ O' + fc-l)! (n-J + k-1)! ^^ 
i!(fc-l)! (n-i)!(fc-l)! 

h+i,k+n-jiK,ml) for n,fc>l. 


(4.52a) 

(4.52b) 

(4.52c) 


The loop-functions 7p q(m^,m^) are result of the one-loop integrals with 
p propagators with mass and q propagators with 


Ip,q{K’K) 


N , 1 dp -^ 59-1 

1671^ (p - l)!(q — 1)! d{m^)P~^ 

(4.53a) 


Ao{ml) -Ao{ml) 


ml-ml 




for 

5'’-% (m2) 


n,o{~ ) (^n — 1)\ 5(m2)" ^ 


for n > 1, 


(4.53b) 

(4.53c) 


This series can be resummed and compared term by term with the deriva¬ 
tives following from the potential (4.50). The complete supersymmetric 
contribution to the one-loop Higgs potential contains also the fermionic 
part (which is easy) and, to be SU(2) invariant, the same diagrams from 
the (s)bottom sector. Once, the potential is fixed (SM minimum at the 
right positions, = v sin (3 and = v cos ^ and the lightest Higgs mass 
m^o = 125 GeV, via e. g. fixing A^), the influence of the bottom squark 
drives an instability below the SUSY scale (see [287] and Sec. 4.2). A 
similar result may hold for the top squark contribution, if Aj is not fixed. 


The results 
obtained in the 
section above are 
part of [287]. 
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4.1.4 Improving the potential 


In the n-loop expansion, the effective potential depends on the renormal¬ 
ization scale Q as tvell as the dimensionless and dimensionful couplings 
A;, see e. g. [219, 230]: 

Veft = V{^i,<P,Q). 


The full potential to 
all orders in 
perturbation theory 
shall not depend on 
the renormalization 
scale. To finite order 
in some expansion 
(like the loop 
expansion) one has 
to wangle for the 
desired behavior 
with the procedure 
described in this 
section. 
Eq. (4.55) is the 
Callan-Symanzik 
equation for an 
arbitrary n-point 
function generated 
byV. 


The renormalization scale Q, however, is arbitrary, though usually chosen 
in a way to keep the corrections (i. e. the logarithms) small. In general, 
the potential should not depend on the choice of this parameter, which 
can be expressed in the differential equation 

Q^V{Xi,4),Q) = 0. (4.54) 


The couplings and fields of the theory are renormalized and therefore also 
depending on the scale Q. Exploiting this fact, we get the renormalization 
group equation (RGE) for the effective potential by applying the chain 
rule to Eq. (4.54) 


Q 


5 




d 





(4.55) 


where we defined the /3 functions for the couplings as 


A(^i)=Q 


dA;(Q) 

dQ 


(4.56) 


and the anomalous dimension of the field (f) 

= (4.57) 

If we are now interested in the effective potential at some arbitrary scale Q 
knowing the potential at the value of the classical field , we easily obtain 
this by solving the RGE for the couplings (i. e. integrating the (3 functions) 
and replacing all couplings in Vgff by the RG improved ones evaluated at 
the scale Q. Eor convenience, let us define the logarithmic derivative 


^dQ dlnQ/Qo dt’ ^ ^ Qo’ 

for some arbitrary but fixed scale Qg. This scale conveniently is chosen to 
be the electroweak scale. Then the parameters of the effective potential 
can be written as 




with the field renormalization ^ expressed as 



(4.58) 


^(t) = exp 


(4.59) 
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Any coupling A;(t) can be determined from the coupling Aj(t = 0): 


r t 


A;(t) — Aj(0) + 




(4.60) 


Jo 

The t-dependence of the renormalization scale Q(t) follows directly from 
the definition of t: dln(Q/Qo)/ d t = 1, 


Q(t)=Qoe^ 


Note that not only the effective potential itself, but also its n-th derivatives 
in <p are scale-independent [267] : 


dV 


H 


eff 


d t 


= 0 , 


with 




d’^ 


(4.61) 


In that way, all the couplings derived from the potential according to 
Eq. (4.61) are scale-independent. 

The efficiency in the use of the RG improvement lies in the very simple 
fact (pointed out clearly by Kastening [288]) that solving the one-loop /3- 
functions and including them in the tree-level formula indeed gives a better 
approximation in the sense that the dependence on the renormalization 
scale vanishes. In a theory without m^-term, as elaborated by CW, the 
RGI effective potential can be written as 


where V{(j) =Q) = / 4! and 


(4.62) 


MQ) = ^ 


3A cj)'^ 
2(47i)2^^^ 


-I -1 


Eq. (4.62) is identical to the GW result 


2^4 






16(471) 


In 


Xcj)^ 


3 


after expanding in the logarithmic term (ln((^^/Q^) 0 around Qf^<p) 

and rescaling 


Q2 



The task to do for obtaining the RGI potential is the following: determine 
the effective potential and solve the RGE for the couplings and evaluate 
the potential at the scale Q(t) = (^(t). This procedure, however, gets 
complicated in case of more fields and more scales related to those fields 
as briefly discussed in the beginning of Ghapter 5. 


There is actually a 
typo in [288] in the 
line after Eq. (5). 
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The Higgs effective 
potential is needed 
to a high accuracy 
to give precise 
determinations of 
the minimization 
conditions, because 
we know that there 
are nontrivial 
minima. 


For the discussion 
about truncated 
series, see 
Sec. 4.2.2 


Also in [291] 
radiative breaking 
was suggested to 
cure the problem, 
but it persists. 


4.2 THE STABILITY OF THE ELECTROWEAK VACUUM IN THE MSSM 

The MSSM is a multi-scalar theory. Not only the two Higgs doublets (with 
also electrically charged directions in field space) but scalar superpartners 
of SM fermions influence the stability of the electroweak vacuum as well. 
The calculation of the full scalar potential at one (or even more) loop(s) 
is a formidable if not impossible task, where the pure Higgs effective po¬ 
tential is already knovra up to two-loop accuracy [289]. 

4.2.1 Distinguishing Between different Instabilities 

The scalar potential of the MSSM possesses any possible dangers, where a 
complete analysis should always include all directions in field space and if 
loop corrections are taken into account, a truncation of the potential after 
renormalizable terms may give a tAirong estimate. 

An exhaustive overview of tree-level instabilities in the MSSM scalar po¬ 
tential was given by the authors of [290] with a useful classification of con¬ 
straints. The existence of charge and color breaking (CCB) global minima 
in SUSY theories was already detected in the early 1980s in the context of 
electroweak breaking in supergravity models at tree-level [291] and via 
radiative breaking [292-294]. Estimates on metastable vacua comparing 
the tunneling time with the lifetime of the universe were also already pro¬ 
posed at that time [295]. However, it is genetically difficult to survey all 
possibly dangerous directions, so specific “rays” in field space are chosen 
(which limit the validity) like {/i°) = {ti) = (t^), which was generalized 
to more realistic configurations where the Higgs vevs do not coincide with 
the CCB minima in [296]. Constraints on soft SUSY breaking trilinear 
couplings A-terms were given by [297] and tan^ bounds in [298]. The 
most general MSSM is even more involved when the full flavor structure 
is taken into account and allows to set even stronger bounds on the flavor 
violating soft breaking terms than FCNC constraints [299]. 

There are in principle three constraints from instable electroweak vacua: 
limits from potentials being unbounded from below (UFB) (also known as 
triviality bounds), such from CCB minima and charge and color conserving 
deeper minima. The latter limit has not yet been discussed in the literature 
and is a genuine one-loop effect, see Sec.4.2.2. It might be that the global 
minimum in such configurations indeed is CCB, the detailed analysis is 
beyond the scope of this chapter. 

4.2.1.1 Unboundedness From Below 

Considering the scalar potential of the MSSM at tree-level 


Vo = Vp+VB+Vsoft, 


where the individual components were introduced in Sec. 2.2.2, one finds 
certain directions in field space (details omitted) where the potential drops 
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dow^n to minus infinity if certain parameters have some specific values. 
The field theoretical potential is related to the total energy density. If it 
is unbounded from below, conservation of energy is violated, thus such 
configurations are forbidden. Fixing those specific rays in field space, one 
obtains constraints on the parameters of the potential, e. g. [300, 301] 

+ — 2|m3|^ > 0, 

> 0 . 

UFB directions can only occur if the mass parameters are chosen inappro¬ 
priately, because quadrilinear and positive F -terms will take over and turn 
UFB to CCB directions [290]. In order to organize the various UFB direc¬ 
tions, the authors of [290] propose three classifications: including only 
(/t°) and (hj), taking one more field into account (where they choose to 
take the lepton doublet) and the (/i^) = 0 limit corresponding to tan^ —> 
oo. 


4.2.1.2 Charge and Color Breaking Minima 

The “traditional” CCB constraints follow from a minimization of the scalar 
potential in certain directions of field space, say (tp) = (u) = (/t°) ^ 0 
and can be expressed in inequalities constraining A-terms [297] 

Aj < 3 (/tIq + , (4.63) 

with similar relations for A^, and A^. Generalizations (which are more in¬ 
volved) for (tL) ^ (t^) were found in [302]. An optimized version of this 
class of bounds was given by [290], where the one-loop effective poten¬ 
tial was included to set the renormalization scale at which the tree-level 
potential has to be evaluated.^ The logic here follows [300], where the 
appropriate scale choice has to be such that the logarithms in Vi are small 
and the one-loop potential vanishes at that scale. The importance of all 
Vi contributions, even small ones, however was emphasized in a different 
context [303]. Eq. (4.63) gives very strong constraints that can be relaxed 
taking the vacuum tunneling rates into account compared to the lifetime 
of the universe [241]; this bound gets modified in view of the Higgs dis¬ 
covery [242]. CCB constraints are a powerful and often investigated tool 
to find limitations on MSSM parameters, see e. g. [304, 305]. Thus, it is 
also quite worthwhile working on an improvement of such conditions. 

Even more important are generalizations of Eq. (4.63) with general, 
flavor-violating trilinear soft terms [299]. The bounds read very similarly 



k = max(i, j). 


(4.64) 


Note the absence of the prefactor 3; the presence of the Yukawa coupling 
in front of the combination of soft masses appears or disappears with the 


3 We have some doubts on that method, because Vj as evaluated in Sec. 4.1.3 has a Q- 
independent part leading to a deeper minimum which was missed in earlier times. 


The parameters mi 
correspond to the 
mass parameters of 
the 2HDM. In SUSY 
theories, is 
given by the soft 
breaking Higgs By. 
term. 

Charged Higgs vevs 
play no role in the 
MSSM. It can be 
shown that the 
minimum in the 
MSSM always has 
{h+) = {h-)=0, 
this is not true for 
other scalars [290]. 


This inequality is 
generally ascribed 
to [291] in the 
literature, but never 
stated there. What 
is given instead, 

“A <3 ”, follows 
from equal soft 
masses (ttiq = 
fhi = m 2 = m) and 
a scaling of the 
trilinear soft 
breaking term with 
the same mass, 

Aj ~ m. 
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We found the same 
conclusions for our 
charge and color 
conserving minima. 

CCB saddle points 
with tachyonic 
squarks, however 
rather hint to 
another more global 
CCB minimum. 
Maybe this one can 
be found more 
efficiently using the 
convex hull. 


As we will see, the 
global minimum of 
the full scalar 
potential is again 
charge and color 
breaking. 


Vq{x) = -m^x^ + Ax^ 


convention for the A-terms. Casas and Dimopoulos [299] use the more 
convenient one for flavor physics tvhich is also used throughout this the¬ 
sis and given in Sec. 2.2.2, where the A-terms are independent of Yukawa 
couplings. Eq. (4.64) and its siblings for the down and charged lepton 
sector give very strong constraints on the flavor-violating soft breaking 
terms, more stringent than most FCNC observables and rule out large val¬ 
ues of Ajj by the demand for a stable electroweak vacuum. Metastability 
considerations weaken the bounds a bit, however, they are still one of the 
strongest constraint’s for flavor violation in the soft breaking sector [306]. 

To be complete, a full one-loop analysis for CCB minima has to be per¬ 
formed which makes it impossible to give analytic results [307, 308]. In 
general, one-loop minimization conditions tend to stabilize the potential 
and its vevs with respect to the renormalization scale as long as the vevs 
are “small” (< 1 TeV) and hint to a breakdown of perturbation theory 
for larger field values. Therefore CCB minima of such large values are 
not trustworthy [308]. Moreover, CCB extrema are found to be rather 
saddle points since one class of scalar mass squares is negative and the 
convex hull [309] shall be taken as approximation for the effective po¬ 
tential where it is non-convex [307]. Similar to our findings described in 
Sec. 4.2.2, no UFB directions occur using one-loop minimization. How¬ 
ever, it was stated that no “alternative MSSM minima” were found and no 
absolute CCB minima [310] using only nonzero stop vevs. We give argu¬ 
ments why this in general (including also sbottom vevs) is not the case. 

4.2.1.3 Charge and Color Conserving Minima 

Another possible class of instabilities in the MSSM scalar potential that 
have not yet been discussed in such a great detail as CCB minima are 
charge and color conserving deeper minima. Such minima preserve the 
gauge symmetries of the MSSM as they appear in the same direction of 
field space as the “standard” vevs, v^ or Vj. The principle occurrence of 
such minima was noted in [310]. We show a viable example of such a 
situation in Sec. 4.2.2 based on the one-loop potential given in Eq. (4.50). 
Adding to the tree-level potential with its standard minima, the 

loop-generated minimum may lie deeper. However, this point is related to 
a tachyonic squark mass eigenvalue and thus rather a saddle point in field 
space. The global minimum is supposed to lie in that direction. 

In Fig. 21, we show an example for a charge and color conserving 
deeper minimum which is loop-induced. The tree-level as well as the 
renormalizable one-loop potential up to ~ show only one minimum, of 
course. The position of that minimum is basically fixed by the quadratic 
term, v^ = m^/A, considering an ordinary potential. On the other 
hand, the pure one-loop potential of Eq. (4.50) itself has a minimum 
which always occurs at a position x > 1. If only the scalar trilinear ver¬ 
tices are taken into account, develops a quadratically rising imagi¬ 

nary part, shown by the dashed curve in Fig. 21. This imaginary part does 
not show up, when also the quadrilinear terms are included, which is a 
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Figure 21: Illustrations of the tree and one-loop potential according to Eq. (4.50). Clockwise from the 
upper left: (a) A generic tree-level potential V (x) = —w?x^ + Ax^. The position of the vev 
can be adjusted with if A is fixed, = m^/A. (b) Re for y = 0. Shown in dashed 
red is the imaginary part which rises for x > 1. The green curve is the real part (anal 5 U:ic 
continuation) of the complex effective potential, (c) with y = x^, for the definitions 

of X and y cf. Eq. (4.49). No imaginary part shows up. (d) Tree + one-loop potential with 
some appropriately chosen weighting factors to make both minima appear. The smaller one 
(X 0.7) corresponds to the minimum of the tree-level potential above, where the deeper 
minimum at x > 1 results from the one-loop potential. 


hint that the full scalar potential in use (Higgs and stop) is well-behaved. 

The complex nature would indicate a non-convex tree-level direction, in 
this case corresponding to a tachyonic stop direction. We show no imagi¬ 
nary part related to the non-convex Higgs potential at tree-level, because 
it is not included in the one-loop part. A tachyonic eigenvalue at position We only consider 

j^Q 

of the second minimum is undesirable and rather hints towards a CCB ^uoop’ ^i-ioop 
global minimum. with h° fields in the 

loop. 

4.2.2 Instable one-loop effective potential with squarks 

For the analysis of the tree plus one-loop effective potential in the MSSM, 
we only include the (s)top/(s)bottom contribution which is dominated by 
the large top (and for large tan (3 also bottom) Yukawa coupling. Contri¬ 
butions from gauge fields are neglected. The only place where we keep 
the gauge couplings is the self-coupling of Higgs fields stemming from 
D-terms. We also neglect D-term contributions in the squark masses. The 
Higgs quartics are needed to have a tree-level Higgs potential which is 
bounded from below and where electroweak symmetry breaking happens. 
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Charged Higgs directions are irrelevant for the discussion of stability [290] 
and by SU(2) invariance, it is sufficient to calculate the potential for 
the neutral components and only discuss instabilities in neutral direc¬ 
tions. Not to violate supersymmetry, we have to include the fermionic 
contributions, as well. Effective potentials for chiral superfields do not 
break SUSY radiatively [311]. The final result in the MS/DR scheme reads 


v,({ = Vo + vi + vl + vf + vl> 


(4.65) 
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(4.66b) 


The tree-level potential Vq is the neutral Higgs part of Eq. (2.29), which 
was written in a way that only 3 contribute to neutral Higgs phe¬ 
nomenology. We want to fix the tree-level mass parameters and 

m 22 in a way that the standard vacuum arises at v = 246 GeV and ex¬ 
pand the Higgs fields around that minimum: 


K = ^i^n + Vn + Un) , K = +^ Xd) ■ (4.67) 

The minimization conditions have to include the derivative of the one- 
loop potential as one-loop extension of Eq. (2.33) 
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Out of the potential, we obtain the mass matrices for the real (scalar) and 
imaginary (pseudoscalar) components of the Higgs field 
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One eigenvalue of pseudoscalar masses has to be zero tvhich corresponds 
to the Goldstone mode. The one-loop effective potential determination of 
the light Higgs mass m^o, however, is not sufficient to deal with the present 
precision. We therefore decided better to calculate the mass of the light¬ 
est Higgs boson with FeynHiggs 2.10.0 [66, 74-76] in order to include 
the dominant two- and three-loop contributions. The connection to the 
effective potential is then made via the pseudoscalar mass m^o which is 
not so sensitive to radiative corrections but rather governed by m ^2 that 
is adjusted to fit m^o from FeynHiggs. We work in the decoupling limit 
rrifyo, m^±, m^o » The trilinear soft breaking coupling in the stop sec¬ 
tor Aj is used to produce the right Higgs mass m^o 125 GeV, where A^, 
is of less importance and can be set to zero. It also has practically no in¬ 
fluence on a the formation of a deeper minimum, neither in /t° nor in 
direction. 

We evaluate the one-loop potential at the renormalization scale Q = 
and set again ihq = = m to simplify the discus¬ 

sion and reduce the amount of arbitrary parameters. A^ is fixed by the 
Higgs mass. A], = 0 and ju as well as tan (3 are kept as free parameters. 
In Sec. 4.2.3, we use the vacuum stability criterion to constrain the pa¬ 
rameter space in the fi-tan (3 plane. Because the scale of interest (Msusy) 
is close by to the electroweak scale, we also do not include the RG run¬ 
ning and use the “unimproved” potential in our discussion. Large values 
of tan 13 are knovra to give a sizable effect on the bottom mass, where we 
actually have to resum the tan/3 enhanced contributions to the bottom 
Yukawa coupling as 


ntb 


(4.70) 


with Aj, = ±ayQ = Ms^sy) in the limit of degenerate SUSY 

masses for the QGD contributions [312]. We are dominantly interested in 
the stop/sbottom contribution to the one-loop Higgs potential, so we take 
those sparticles to be rather light (~ 1 TeV) where especially the gluino 
is expected to be rather heavy—so Aj, is small and does not alter the 
bottom Yukawa coupling much. On the other hand, still gets signifi¬ 
cantly altered taking the SUSY threshold corrections into account. While 
the gluino contribution decouples with large gluino mass Mg, the same 
from the higgsino does not for increasing higgsino mass. Though the hig- 
gsino part is numerically smaller, it is not negligible. Both contributions 
sum up together, Af, = where [312-315] 


2.CC 

= -^MA^Gtan/ 3 Co(m 5 ^,m 5 ^,Mg), (4.71a) 

^higgsino ^ _^piAttan^Co(mE^,mf^,/i). (4.71b) 

There are also bino and wino contributions which we ignore due to our 
policy of ignoring gauge couplings in all loop contributions. As can be 
seen from Eq. (4.71), the bottom Yukawa coupling (4.70) gets enhanced 
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Likewise, the 
one-loop 
diagrammatic 
calculation gives 
too poor results. 


The electroweak 
contributions from 
charged 

Higgsino-stop and 
Wino are important 
for large pA^ and 
pM2. 


For illustration, we 
show the behavior 
of Af with the 
gluino mass: 


\ 



And the same for 

varying from 
-5TeVfo-|-5TeV 
(we take 

rrisoft = 1 TeV and 
At = -1.5 TeV) 














Veff/GeV' 
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Figure 22: The effective potential Vq + Vi in the Vu direction develops a second minimum for suitable 
configurations. Especially the /r parameter of the superpotential drives the steepness of this 
instability. From the requirement of both vacua being degenerate, we can derive a bound in 
the /r-tan/3 plane for which the standard electroweak vacuum is unstable. The parameters 
for the plot on the left side are fixed to tan [3 = 40, Mgusy = 1 TeV and m^o = 800 GeV to 
coincide with [287]. We have taken the negative AfSolution for the light Higgs mass (from 
— 1416GeV to —1468GeV in the displayed curves) whereat, = 0. 


The bottom mass 
and therewith the 
Yukawa coupling is 
a running (MS) 
parameter which is 
very sensitive to the 
scale choice. 

The assignment 
signAf = — sign p 
can be seen from 
Eqs. (4.70) and 
(4.71b) which 
enhances Yj,. 


for a negative sign of ju. This effect can be reverted with a negative Aj for 
the higgsino case. 

The bottom resummation including the higgsino contribution was mis¬ 
takenly ignored in [287], which is obviously wrong because it has an im¬ 
portant effect. Moreover, in [287] the tree-level bottom Yukawa coupling 
was not evaluated at the SUSY scale but at the scale of the bottom mass. 
Compared to the bottom Yukawa coupling at the SUSY scale, it was about 
a factor 1.6 too large. On the other hand, including signif¬ 

icantly enhance again and find the same observation of multiple vacua 
in roughly the same regime of ju-tan ^ —however, this only works if jj, and 
At have a different sign! In such a way, we can give quite complementary 
constraints to what is usually obtained. 

4.2.3 Constraining the Parameter Space by Vacuum Stability 

The one-loop effective potential shows a strong dependence on the p- 
parameter which is displayed in Fig. 22 on the left hand side. Small 
changes in this parameter obviously lead from a stable configuration to 
an unstable one. It is intriguing to figure out the “multiple point”, so the 
value of p for which the two vacua are degenerate. This corresponding 
curve in the /i-tan (3 plane gives an exclusion contour: everything to the 
upper right corner is excluded by the formation of a second, deeper charge 
and color conserving minimum. 


VACUUM DECAY Let US briefly address the question whether this at first 
sight unstable configuration is really instable or rather metastable in a 
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Figure 23; We approximate the potential barrier with a triangle according to [316]. We denote the 
position of the local maximum with and of the higher (instable) minimum with 0_|_. The 

global minimum is at </>_. Correspondingly, the values of the potential are denoted with 
and V. 


cosmological sense. How large (or small) is the decay time of our vacuum 
at 246/-/2GeV? Finding the bounce reduces to finding solutions of the 
bounce differential equation which can be complicated. Numerical tech¬ 
niques already suffer from the boundary condition = 0” together 

with (j)'{p)/p in the differential equation and may be cured using non¬ 
standard techniques. There is an easier-to-use possibility for an estimate 
on the bounce approximating the potential barrier with a triangle [316]. 

The decay probability per unit volume was given hy F/V = Ae“®^^ (in 
the following we work again with h = 1), where the coefficient A can be 
either estimated by (100GeV)^ or the barrier height, AV^ = “k'+- 

In [316] a pocket calculator formula is given for the bounce action 


2n^ - {A4>_fy 

~3 


(4.72) 


Note the scaling 
factor of yf2 
according to the 
definition in 
Eq. (4.67). 


We define 

A(p^ = (^jvL “ 0 + 

and 

A(p- = (p_- 


We determine B using the triangle method for the potential in Fig. 22 with 
p = 3910 GeV: 


X 

4>x 

ReVx(</)x) 

+ 

174GeV 

-5.1 X lO^GeV^ 

M 

373 GeV 

2.1 X 10® GeV^ 

- 

484 GeV 

-3.6x lO^GeV'^ 


and obtain Sg [(^bounce] ^ 18.8 which is a ridiculously small number com¬ 
pared to (meta)stable vacua with > 400. The decay time follows from 
r/V after multiplication with the volume of the past light-cone which is 
basically the age of our universe to the fourth power, T^. We have 

AV ] R 111 

^vac = ^ 2 X 10-111 Tu. 

If we were in such a configuration, our vacuum would immediately de¬ 
cay. Neighboring regions in parameter space behave the same. Especially 


Ridiculously, 
because we have to 
compare the 
exponents, 

gl8.8/g400 py 

3 X 10"!®^. 
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the difference — {A<p_)^ is roughly the same because the sec¬ 

ond minimum is driven by the one-loop vev. For deeper second minima, 
B gets even more reduced because AV^ increases. The transition from 
metastable to instable in cosmological terms happens very rapidly as a 
function of ju. 

CHARGE AND COLOR BREAKING GLOBAL MINIMUM Do We knOW^ "what 
the global minimum is? In general, it is arbitrarily difficult to find the 
global minimum of a multivariate polynomial—in our case the function 
contains also logarithms. Numerical algorithms may hang up in a local 
minimum or overshoot the global one. However, in our special case de¬ 
picted in Fig. 22, we can find out whether the deeper minimum is an 
impostor or not. 

We have argued that the tunneling to the deeper minimum happens im¬ 
mediately, so the theory has to be expanded around the new vacuum. A vev 
of about 600 GeV influences the masses. The fermion masses are directly 
~ (/i°), so the top quark is just heavier. But the squark masses have a non¬ 
linear dependency on the vevs: while the stop mass matrix scales in the 
diagonal as well as the off-diagonal elements with /t°, the sbottom mixing 
gets significantly enhanced via the juYb/t^-term in the off-diagonal—which 
drives one eigenvalue negative and therewith one eigenstate tachyonic. 

A negative sbottom mass squared means that the potential (its second 
derivative with respect to the fields, in this case the sbottom fields, gives 
the mass matrix) has a non-convex direction aligned with the correspond¬ 
ing mass eigenstate. Non-convex potentials, however, are a sign that we 
are expanding around a wrong point and the more global minimum is 
expected to lie in direction of the non-convexity. In this case, the global 
minimum is charge and color breaking and related to a sbottom vev (b-^). 

This observation is clear in regions excluded by the deeper second min¬ 
imum shown in Fig. 22 and gets unclear if the smaller vev is related to 
the deeper minimum. To figure out whether the CCB minimum in the 
tachyonic sbottom direction really is deeper, so if V((b^)) < V (vy), needs 
investigation of the full scalar potential and goes beyond the scope of this 
thesis. One point to check remains, which are the tree-level (and “em¬ 
pirical” or improved) CCB constraints of [241, 242, 297]. It is easy to 
check that for the parameter point in use (Mgysy = 1 TeV, ju = 3910 GeV, 
At = —1468 GeV) the CCB bounds of [242, 297] are fulfilled where the 
“empirical” one of [241] is violated. This conclusion about tree-level CCB 
minima stays unchanged compared with [287] despite of the updated 
numerics. We therefore assume to be safe from CCB minima of the tree- 
level potential without further investigation. Nevertheless, we want to 
stress that we get a loop-induced CCB minimum depending on the larger 
{/i°) > Vy which cannot be accessed using standard tools [317, 318]. 
A more careful and quantitative analysis of this property and potential 
stronger bounds on CCB minima has to be done as a follow-up of [287]. 


COMBINING NEUTRINOS AND THE VACUUM 


In Chapter 3 we have discussed neutrino flavor physics in the context of 
supersymmetric theories, where Chapter 4 was devoted to the question 
of the electroweak vacuum being stable under supersymmetric quantum 
corrections. We have restricted ourselves in Chapter 4 to the contributions 
from the third generation (s) quarks because their coupling to the Higgs 
fields is governed by large Yukawa couplings. The instability which we 
found appears close to the electroweak scale and can therefore be dis¬ 
cussed and analyzed without reference to any high scale physics. 

It is well-known that the SM effective potential reveals a metastable 
vacuum state in the light of present data [274] where the principle scale 
of instability is determined by the point where the Higgs quartic coupling 
turns negative [283, 319, 320]. Also it is well-known that the presence of 
seesaw neutrinos alter the statement about stability, instability or metasta¬ 
bility [274, 285]. 

In this Chapter, we shall examine the influence of heavy Majorana neu¬ 
trinos to the stability discussion we deduced for the SUSY corrections. 
Working with neutrinos, it is an important and necessary task to con¬ 
firm the stability of the low-energy theory. Eurthermore, it is an inter¬ 
esting proof to see whether some high-scale dynamics may render the 
low-energy vacuum unstable. Going this way, we get information from a 
very high scale which will never be directly accessible by experiment. 

There is, however, an issue to be treated with care. Eor the quark-squark 
one-loop contribution to the Higgs effective potential, we have identified 
a new vacuum arising just behind the corner at a scale where the loga¬ 
rithms of the one-loop potential are small anyway. This naive estimate fails 
as soon as particle masses in the loop differ very much from the (fixed) 
renormalization scale. Already for the stability discussion of the SM effec¬ 
tive potential evaluated for classical held values around the Planck scale, 
one has to resum the large logarithms by means of the renormalization 
group. In this approach, the renormalization scale is not taken at a fixed 
value (as was e. g. done in Chapter 4 where we chose Q = Msusy) but 
Q ~ <p to improve perturbation theory and reduce the size of the loga¬ 
rithm ln((^/Q) [267]. The several-scale approach was addressed in the 
literature and solved very elegantly and easy to implement with the de¬ 
coupling method [284, 322] which is briefly reviewed in the following 
and then applied to our problem of SUSY with type I seesaw. The decou¬ 
pling of heavier degrees of freedom allows us to consider at low scales 
only contributions from below the threshold. The scale dependence of the 
effective potential part stemming from those heavy fields is mild in a way. 



In general, in a 
one-field problem 
one chooses 
Q = M^{f) [321], 
which in massive 
-theory turns to 
Q~ <p at large 
values of the 
classical field ^. 
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However, already 
the light neutrino 
mass drastically 
destabilizes the 
potential 
introducing a 
UFB-direction as 
pointed out in 
Sec. 5.1. 


Decoupling is 
meant in a very 
simple sense: the 
corresponding part 
is just set to zero. 
Also, at the scale 
q2 = M^{4>), the 
logarithm itself is 
zero, so no 
discontinuity is 
introduced in the 
one-loop potential 
at the threshold. 


that it only enters through the parameters in the 1-loop part and therefore 
is of higher order. 

The issue of instability is related to large field values since the effec¬ 
tive potential at low field values is determined basically by the SM and 
new physics appearing at some higher scale alters the behavior of that po¬ 
tential. Heavy neutrinos (and sneutrinos) start to play an important role 
above their masses—below, the heavy states are integrated out and the 
light neutrino masses are suppressed by the scale in the spirit of the see¬ 
saw mechanism. If right-handed neutrinos are well separated in mass and 
show a strong hierarchy as well as when we want to incorporate standard 
SUSY contributions from the (s)top-(s)bottom sector, the choice of Mr (or 
say the heaviest Vr) as the overall scale seems unreliable because now we 
would introduce exactly large logarithms as log{MsusY/. The choice 
Q = Mr as universal scale choice for all renormalization scale dependent 
quantities is in such a case a bad choice. 

A proper and efficient way to cope with different scales and deal with 
the decoupling of heavy particles was suggested in Ref. [284]. The effec¬ 
tive potential contribution of fields with independent masses M^ and M 2 
can be given by 


Vi = 


6471 ^ 


^ M? „ 

iViM^ log 01 + log 02 


(5.1) 


where 0; = 0(Q —Mj(0)) is the step function taking care of the decou¬ 
pling below M; and the prefactors iV; account for degrees of freedom. The 
influence of those fields to the ^-functions comes in the same way: 

d A 

(5.2) 


where denotes the contribution from field i to the ^-function for some 
coupling A. 

In general, this can be generalized to an arbitrary number of heavy 
fields. For our purpose, however, two are enough, where one will be iden¬ 
tified with the scale of heavy neutrinos and the other one with the SUSY 
scale. 

The advantage of Eq. (5.1) compared to multi-scale approaches where 
for each heavy mass an individual renormalization scale is chosen [323] 
is the occurrence of only one global renormalization scale (which has to 
be properly adjusted) and a simple inclusion of /3-functions according to 
Eq. (5.2). The authors of [284] argue that two arguments lead to the 
preferred scale Q* = Q = min;[M;((^)}. Eirst, the complete potential 
V = Vq + to be scale independent, sody/dQ = 0 (which is the 

reasoning behind the RG improvement). Second, they want to choose Q in 
such a way that the loop expansion behaves best, especially AV = V — Vj 
is small (or Vi =0). At the scale Q*, tree-level and one-loop potential 
are identical and only the RG improved tree-level potential has to be eval¬ 
uated with running couplings at the scale Q*. Note that the ^-functions 
look differently in each regime. 
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The improvement of the one-loop potential V = Vq + Vi needs the 
use of two-loop RGE, which are available for the MSSM extended with 
right-handed neutrinos [164, 165], To get a basic feeling of the neutrino- 
sneutrino influence on the Higgs effective potential, we constrain our¬ 
selves in the following to one generation of leptons. 


5.1 NEUTRINOS DESTABILIZING THE EFFECTIVE POTENTIAL 


Neutrinos are fermions and therefore destabilize the effective potential 
anyway by the negative prefactor. It is very obvious and easy to see, that 
already light seesaw neutrinos significantly destabilize the Higgs potential. 
A clear observation is that light neutrinos enter dramatically, once the 
classical Higgs field value is not restricted to be at the SM vacuum. If 
we develop the potential for large (/i°), the neutrino contribution grows 
quadratically with the field as can be seen from the seesaw formula. Let us 
consider a simple type I seesaw with a Majorana mass for the right-handed 
neutrino, see Eq. (2.38): 


-£ 


V 

m 






+ h. c. , 


(5.3) 


using left-handed Weyl spinors, and the eigenvalues of My{h^) are 

(K) = l {mr ± • ( 5 . 4 ) 

We expand in the large right-handed mass Mr and obtain the well-known 
seesaw formula for the light neutrino, while the other one stays heavy 


0\2 




Mp 


^ Mr. 


(5.5) 


Below the scale Mr, only the light fields are active and only those con¬ 
tribute to the effective potential. The heavy contribution is kept away with 
a theta function in the spirit of Eq. (5.1). We want to have a light neutrino 
below 1 eV, so let us take for the moment Mr = lO^^GeV and Yy = 1 in 
order to have mY^(174GeV) 0.3 eV. The light neutrino Higgs effective 

potential is then given by 




327i2 


1 (Yyhir 

32n^ 


nr 


In- 









{Yyhy 

M2q2 



(5.6) 


We keep the 
notation of the 
2HDM, where one 
Higgs doublet 
couples to up-type 
fermions and 
therefore also to 
neutrinos, which 
we called Hy. Its 
neutral component 
is ?i°. For the 
discussion in this 
section where only 
one direction in 
scalar held space is 
of interest, it does 
not matter whether 
there are more 
fields. We also only 
take a look at the 
pure one-loop 
effective potential 
and its behavior at 
large field values. 
The tree potential 
can be estimated as 
in the large 
field regime. The 
question of 
(in)stability reduces 
then to the question 
whether |Vj^| > Vg 
at some point. 


There are some issues to be clarified: which renormalization scale we 
shall take and what to do with the RG scaling of the involved parame¬ 
ters. Moreover, the effective theory is not the correct description at scales 
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Figure 24: We plot the effective potential in a smart way inspired by [280], where V ^ = 1^0 + and 
Vq = On the left side, we only add the effective neutrino mass operator. The plot 

on the right side has the correct UV completion. In both plots, the red line shows the RGI 
tree-level potential, which on the left side is unaffected by the heavy scale because its /3- 
function does not feel the presence of the heavy scale. On the right side, the tree-improved 
potential drops negative beyond the scale of right-handed neutrinos. In both cases, the one- 
loop potential is unbounded from below as long as no new physics at the Planck scale (or 
below) rescues it. The other lines show specific choices of the renormalization scale: the 
fixed values Q = ITeV (green) and Q = = 10^"^GeV (blue). In the case where the RG 

evolution of A is not influenced by heavy fields, the scale choice in basically makes no 
difference and the instability scale is always at the same point with an uncertainty about one 
quarter order of magnitude. This changes once the running is altered: Obviously, the fixed 
scale choices either underestimate or overestimate the instability scale. If we take Q = 
(orange) as the value of the classical field, the RGI tree potential and the one-loop potential 
drop negative at roughly the same scale. 


The one-loop RGI 
potential actually 
has to be evaluated 
using two-loop RGE. 


larger than Mr, so in any case the UV completion has to be properly con¬ 
sidered in the stability discussion. For the latter point, we restrict our¬ 
selves to the simple type I seesaw of Eq. (5.3) where only one additional 
fermionic state appears at the high scale. To address the RG behavior of 
the one-loop effective potential (5.6), we have to take into account the 
self-coupling of the tree-level part. The neutrino mass operator can be 
written as = (/i°)^k:, where k has mass dimension —1 and a cer¬ 

tain running in the SM and MSSM [181]. Below the scale Mr we just have 
the /3-function for k in addition. Above, the other ^-functions get altered. 
The relevant RGE are given in App. B. 

Note that the fermionic contribution to the Higgs effective potential is 
always negative. So the light neutrinos as well as the heavy ones turn 
the potential negative at some point and the potential will be unbounded 
from below. We show the effective potential in some variations in Eig. 24. 

We compare in Eig 24 several choices for the renormalization scale at 
which running parameters are evaluated—and find agreement to what 
was genetically proposed in the literature (namely Q ~ ^). Eor that pur¬ 
pose, we discuss the tree and the one-loop potential (unfortunately only 
one-loop RGE are available for the Higgs self-coupling under presence 
of right-handed neutrinos) with running couplings evaluated at certain 
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scales. The high-energy (or high field value) behavior of the effective po¬ 
tential is governed by the quartic coupling, so we only have A(/t°)'^ as 
tree-level potential and therewith 


^eff(^S)=^0 + V 


0\4 


= ^(^u) 


32712 




!= 1,2 


Q2 


with given in Eq. (5.4). First of all, we cannot take Q = m{(p) 

to make the logarithm vanish as usually suggested because there are two 
different field dependent masses in the logarithm. The closest suggestion 
would be to set Q = Mr which, however, also produces large logarithms 
in the regime /t° ~ Mr (which grow for /i° > Mr) and is rather a fragile 
choice. Similar considerations hold for a fixed but low scale, Q = 1 TeV. 
The proper renormalization scale at which the effective potential has to 
be evaluated is indeed Q ~ /i°, for which the logarithms get smaller be¬ 
yond Mr. In Fig. 24, we show these three choices of the renormalization 
scale Q in the one-loop potential with running couplings evaluated at the 
same scale. Additionally, we display the pure tree-level potential with a 
running self-coupling. When we include only the effective light neutrino 
mass given in Eq. (5.5) (which is obviously wrong), the RGI tree potential 
is unaffected by any heavy scale. If we include the effect of right-handed 
neutrinos in the running of A, we already reproduce the same UFB be¬ 
havior as the one-loop potential shows—and the difference from the full 
inclusion of is rather mild in the shift of the instability. 


5.2 SNEUTRINOS STABILIZING THE EFFECTIVE POTENTIAL 

The one-loop effective Higgs potential in the presence of neutrinos is UFB 
below the Planck scale. For this main statement there is no difference if 
only an effective “light” neutrino Majorana mass is taken into account or 
a UV complete theory including heavy Majorana neutrinos. The statement 
may get altered if a different UV completion is considered like a type II 
seesaw inspired model with heavy scalars: scalars contribute to the effec¬ 
tive potential with a positive sign where fermions always give a negative 
contribution. 

How does the situation change when the theory is supersymmetrized? 
In a SUSY theory, bosonic and fermionic degrees of freedom are equal, so 
for exact SUSY one would expect the loop contribution of neutrino super- 
field components to be exactly canceled. As pointed out already before in 
Chapters 2 and 3, the seesaw mechanism is transferred to the sneutrino 
sector as well and we expect heavy sneutrinos at the scale Mr. We only 
consider one generation of neutrinos. A supersymmetrized version with 
heavy right-handed neutrinos at say Mr = 10 ^^ QeV is completely stable 
concerning the influence of neutrinos. However, at this stage we can say 
nothing about stability against sneutrino vevs. 


The alternative 
would be to take 
the multi-scale 
approach by [323] 
and choose 
different scales for 
both logarithms. 
Although the 
ren ormaliza tion 
scale is an 
unphysical scale it 
gets a more physical 
meaning via the 
relation to the held 
strength of the 
classical “external” 
Held, Q = (^. 


Our purpose, 
however, is not to 
survey all popular 
neutrino extensions 
of the SM but rather 
to check whether 
the simple SUSY 
seesaw type I 
discussed and 
exploited in 
Chapter 3 leads to a 
potentially unstable 
electroweak scale 
vacuum. 
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In the following, we show the neutrino-sneutrino effective potential, 
where the analytic calculation can be very easily done and deduce its sta¬ 
bility below the scale of right-handed neutrinos. The results are eventu¬ 
ally not surprising: what we shall see is that the large separation of scales 
(MgusY vs. Mr) makes any negative fermionic contribution to the effective 
potential vanish, because SUSY is more or less unbroken at the high scale 
and the splitting in the mass spectrum of ©(Msusy) is a small perturbation 
which actually plays no role compared to the large amplitude of the Higgs 
field, we are interested in. 

We set the stage with the following superpotential (we omit flavor in¬ 
dices and explicitly work with only one generation) 

W D pHd • Hu + Hu • LlHr - Hd • (5.7) 

from which the sneutrino mass terms can be calculated as the bilinears in 
the f-terms (see Sec. 2.2): 


VF{4>,r) 


SW'*' dw 


(5.8) 


where cj) is the scalar component of a chiral supermultiplet. The Higgs and 
lepton superfields were already introduced in Sec. 2.4. 

In the exact SUSY limit, the mass terms of the scalar neutrinos are de¬ 
termined by Eq. (5.8). However, since SUSY is broken, we introduce the 
common soft SUSY breaking Lagrangian for one generation of sneutrinos: 


^soft = VrV* + (a^ VlVr + + h. c. ) . (5.9) 

The mass matrix can then be written in a four-dimensional basis: 


-C= o (v^vL,VR,vy 


M 


LL 


{mIfT m 







where the 2x2 Higgs field dependent sub-matrices are given by 


= 


■^LR ~ 


■^RR ~ 


mr 


\yX\" 


0 

k*y:mr 


- nh[ 


ml + \Yyh' 


0 
d 

0|2 
v^ul 

2B^ 



h^Y M* 


(5.10a) 

(5.10b) 
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2{B^) 
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(5.10c) 


Note that we neglected D-term contributions ~ M| proportional to gauge 
couplings. In the following we will restrict ourselves to the one-family 
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case and also keep the parameter out of the discussion for convenience. 
Nevertheless, its presence does neither change the RGE for Mr nor does 
it directly contribute to the Higgs potential. It may be of interest once we 
are interested in sneutrino vevs at a high scale which is not our purpose at 
the moment. 

For a real mass matrix in the limit = 0, the four mass eigenvalues 
can be directly calculated 


(/t«) = - + 2M| + |/,o|2 ^ ^M^ + 4\himA^±M^Y^y 


ttt?, \'Ln, 

2/ 2 


(5.11) 


The = 0 scenario 
leaves the [i-term 
out of the 
discussion. 


assuming left and right soft masses being equal: M? = m? = m|. 

In the following discussion, we also ignore the trilinear coupling Ay 
which is relatively suppressed by 1/M^ as can be seen from Eq. (5.11). 
Assuming all SUSY-breaking parameters being 0{Mg), Ay would only be¬ 
come important if it was O (Mr) . 

The field dependent masses can be nicely written down in a way to 
estimate their importance compared to Mr using e = 4Y.^|/t2p/M^ as 

my^^ = ^[l±^/l + ~8y (5.12a) 

Ml f 2M| ^ . -A 

rti- — —— I 1 H-^—h ~ i Y 1 + e . (5.12b) 

-4 2 \ Ml 2 ^ j 


We are counting 
from light (1) to 
heavy (2 or 4). 


We interpret Mr and Yy as running parameters evaluated at some scale 
Q, where we take Q = /t°. The fermionic contribution to the effective 
potential gets a factor two due to the spin degrees of freedom. In our 
approximation, however, also and rriy^ = rriy^ and therefore 

the individual sneutrino contributions can be seen as one light and one 
heavy—garnished with a factor of two. So in the SUSY limit (Mg —> 0) 
both contributions indeed exactly cancel. With our knowledge of the multi¬ 
scale treatment of the effective potential, we state 

- K, { log(^vi / (^u)^) - 3/2) 

- (log(m2y (h0)2) - 3/2) 0 - Mr) ^ 

+ / (/t«)2) - 3/2) 0 (/t« - Mg) 

+ / (/t«)2) - 3/2) 0 (/t« - Mr) . 

Approximately, rtiy^ rtiy^ ph Mr and an extension to three degenerate 
(s)neutrinos can be made via multiplication with three. Eq. (5.13) does 
not contain a cosmological constant term because my^{h^) vanishes for 
/i° —> 0 and the other contributions are cut away from zero with the 0- 




3271 ^ 


For simplicity we 
use a fixed step in 
the 0-functions 
with either Mg or 
Mr. And ignore the 
discontinuities by 
thresholds in the 
following discussion 
anyway where we 
are only interested 
in the behavior of 
the full one-loop 
effective potential 
above the heavy 
neutrino threshold. 
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function. The analytic expression for the one-loop effective potential fol¬ 
lows from Eqs. (5.12a): 
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(5.14) 


Also the potentially 
large logarithm 
\n{M?/Ml) in the 
part is 
suppressed by the 
small prefactor 


For the discussion around the scale of right-handed neutrinos, we take 
= M^: Eq. (5.14) suggests this scale choice since the only large loga¬ 
rithms like ln(M|/M|) appear in the “cosmological constant” piece which 
is not present anyway because cut away by the 6-functions. We expand 
in £ = 4y^|/t°p/M^ to figure out the dominant behavior below the scale 
Mr— and find positive coefficients: 


3271^ 

Mi 


-v:-'' = 




(5.15) 


Beyond /t° = Mr or e = 1, the expansion in s breaks down and only the 
complete sum gives the appropriate result. Anyhow we do not need the 
expansion beyond: the potential is not driven into an instability before the 
heavy states enter the game and it shall not beyond. As stated before, we 
add and subtract basically the same and the one-loop potential is indeed 
given by the RGI tree-level potential. The running of the gauge couplings 
is not altered by the heavy Majorana neutrinos as was the running of the 
Higgs self-coupling in the non-SUSY theory of Sec. 5.1. We even cannot 
show any propaganda plot, because there is no propaganda to show. 
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The outcome of the preceding section is not overwhelming and even not 
surprising. However, we can state that a supersymmetrized version of a UV 
extension with SUSY broken at the TeV scale or any scale well below the 
scale of new physics is stable against further vacua in the SM vev direction. 
As we stated above, we cannot say anything about sneutrino vevs at the 
moment, especially when all the soft breaking A- and B-terms are taken 
into account. We dare to generalize our findings from type I seesaw with 
Majorana neutrinos to any not yet thought of theory. As long as SUSY is 
broken at a much smaller scale, the contributions to the effective potential 
at higher scales are canceled. 

Without knowledge of high-scale physics, especially without knowledge 
of any quantum-gravitational interaction around the Planck scale, there 
is no statement about absolute stability of the effective potential and the 
electroweak vacuum possible as was pointed out in [282]. However, an 
exact supersymmetric theory does not introduce further instabilities. More¬ 
over, the running of the Higgs quartic is determined by the running of the 
gauge couplings (squared) which never run negative. We therefore con¬ 
clude that any SUSY theory is expected to be stable beyond the Planck 
scale with respect to SM-like minima. Even softly broken SUSY is approx¬ 
imately exact up there and therefore no such second minimum at a high 
scale as in the SM [274] shows up. The dynamics of the electroweak 
vacuum at the electroweak scale is then only determined by SUSY scale 
physics. 


The inclusion of the 
full scalar potential 
in the (y)MSSM 
respecting all 
possible sfermion 
vevs is beyond the 
scope of the present 
discussion. 


The authors of 
[282] state that any 
higher-dimensional 
operator (of 
dimension six, 
eight) suppressed 
with the Planck 
mass changes the 
behavior below the 
Planck scale. 
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SUMMARY OF CHAPTERS 4 AND 5 We have explicitly recalculated the 
one-loop effective Higgs potential in the presence of third-generation squarks 
and have found the formation of a minimum at one-loop order in Chap¬ 
ter 4. In the combination of tree-level and one-loop potentials, the loop 
induced minimum may appear deeper than the standard minimum in the 
direction of one neutral Higgs component (we showed results in the /in¬ 
direction) . The new minimum, however, gives tachyonic sbottom masses 
which indicate a global charge and color breaking minimum with (b) ^0. 
From the requirement of both minima being degenerate, we can formu¬ 
late an exclusion limit on the parameters of the theory and we have 
explicitly shown such limits in the /i-tan^ plane. Heavier SUSY masses 
shift the limit to larger values of both ju and tan /3. The influence of tan (3 
resummation on the bottom Yukawa coupling leads to the requirement 
sign At = — sign /i on the relative signs of At and fi to produce this obser¬ 
vation. We are therefore quite complementary to existing bounds on these 
parameters from vacuum stability. 

The same calculation of the effective potential in presence of (heavy) 
neutrinos and sneutrinos has shown that no further instability is intro¬ 
duced in the SUSY theory. The non-SUSY description of neutrino masses 
(the effective theory as well as the UV completion with right-handed Ma- 
jorana neutrinos) instead results in an effective potential which is un¬ 
bounded from below up to the Planck scale. The SUSY version in contrast 
is well bounded from below. 


MIXING ANGLES FROM MASS RATIOS 



Be content with what you have; rejoice in the way things are. When you 
realize there is nothing lacking, the whole world belongs to you. 

—Laozi 

An amazing amount of work was put into a deeper understanding of 
our three matter families, their mixing and why their masses are so differ¬ 
ent [324-367, we for sure only refer to a small fraction of papers dealing 
with that topic]. In this concluding chapter, we shall address a minimalis- 
tic approach slightly orthogonal to what was discussed in Ghapter 3: how 
much do we need to know about mass matrices and how many assump¬ 
tions do we need to impose in order to get a viable description of flavor 
mixing. The results of this chapter have been published in [205], 

The gauge structure of the SM defines the largest global flavor symme¬ 
try which is allowed reflecting the remaining symmetry if Yukawa cou¬ 
plings are switched off: [U(3)]^. There is one U(3)-factor for each gauge 
representation; taking right-handed (singlet) neutrinos into account, we 
have 

[U(3)]'’ = U(3)q X U(3), X U(3)d x U{3), x U(3), x U(3)^. 

The Yukawa couplings, which break this maximal flavor symmetry group, 
are strongly hierarchical and differ over several orders of magnitude 

Yu : Yc: ^ 10 “^ : 10 “^ : 1 , Ya : Y^ : Yt : 10 “^ : 1 , 

Ye : Y^ : Y, : 1 . 

The smallness of the first and second generation Yukawa couplings allows 
to impose a smaller symmetry group, [U(2)]^, whereas only first genera¬ 
tion vanishing Yukawas lead to [U(l)]^. We propose a minimal breaking 
of maximal flavor symmetry by the following symmetry breaking chain: 

U(3) ^ U(2) ^ U(l) ^ nothing. (6.1) 

The scales A; at which the symmetry breaking occurs may be largely sep¬ 
arated: A 3 » A 2 » Ai similar to the hierarchy of Yukawa couplings 
Y 3 » Y 2 » Yi for third, second and first generation. 

We remark the similarity to radiative mass and flavor models briefly 
discussed in Sec. 2.2.3. On our approach, however, we try to be as model- 
independent as possible and only require Eq. (6.1) happening simultane¬ 
ously in the up and the down sector (i. e. = Af = A^). 

The principle of minimal flavor violation relying on U(3) and U(2) sym¬ 
metries and discrete subgroups was applied in many aspects of flavor phe¬ 
nomenology e. g. in [89, 104, 346, 366, 368-383] 


The application 
e. g. of RFV 
techniques may be 
even simpiified with 
the procedure 
presented here. 


In the gaugeiess 
iimit, we observe 
U(45) orU(48), see 
Sec. 2.1.1. 


Fermion masses are 
directiy 
proportionai to 
Yukawa coupiings, 
nix = vYx! 


We omit obvious 
trivial U (1) factors 
that are left: aiso in 
the SM with ail 
Yukawa couplings 
non-vanishing, 
there are still 
accidental 
symmetries known 
as baryon and 
lepton number. 
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The lower-rank 
approximation 
theorem is known 
as Eckart-Young- 
Mirsky or simply 
Eckart-Young 
theorem, though it 
is better to call it 
Schmid t-Mirsky 
theorem [388]. 


Singular values of 
mass matrices 
correspond to the 
physical fermion 
masses. Each vector 
of left- and 
right-handed 
fermions is rotated 
with the left- and 
right-singular 
matrix into the 
mass eigenbasis. 


Eq. (6.3) is a 
rank-one matrix, 
labeled with 

Inversion of 
Eq. (2.7), 
m = SimS^, gives 

m = S^itiSr. 


6.1 hierarchical mass matrices 


The hierarchy of fermion masses following from the minimal breaking 
of maximal flavor symmetry allows us to reversely approximate fermion 
mass matrices by matrices of lower rank [384-387] . The matrix rank gives 
the number of linearly independent columns or rows of a matrix. A generic 
three-generation mass matrix has rank = 3, a U(2)-symmetric mass ma¬ 
trix only rank = 1. We decompose the generic mass matrix by its 
singular value decomposition as of Eq. (2.7) with left- and right-singular 
matrices S-^ and S^, respectively, that are build up of their singular vec¬ 
tors 


c-f — 


-/ -/ 




pff pff pft 
^R,l’*R,2’*R,3 


and find the diagonal matrix ifi 
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■/ -/t 
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^L,l f-’R,! 


-/t 

^L,3*R,3 




m{, 


( 6 . 2 ) 


with the singular values that obey > 0. 

We now see that with a realistic fermion mass spectrum, » 

m^, the mass ratios in Eq. (6.2) are small, and we get a rank-one approxi- 

f f f f 

mation by neglecting both /nig and 


fh.^ = ' 

" r=l •^L,3'^R,3 


^0 0 0 ^ 
0 0 0 
0 0 1 


(6.3) 


where we divided by the largest mass nig. We keep the notation “ for nor¬ 
malization with respect to trig. Correspondingly, neglecting only /h-[, we 
find the rank-two approximation which has a potentially arbitrary 2x2 
sub-matrix left reverting the singular value decomposition: 


"V=2 *L,2'"2*R,3 ' *L,3*R,3 


f 0 
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0 nr 
0 ffr 
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22 
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32 
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23 


(6.4) 


A closer look at the 2-3 block of Eq. (6.4) shows a hierarchy in its elements 


m: 


/ |2 
331 


m: 


/ |2 U / |2 


23 I 


m: 


321 


m; 


/ |2 


22 I 


Confining ourselves only to an order-of-magnitude discussion, we neglect 
all contributions O ( 17122 ) “ ^ (|^ 23 1 ) work with the approximation 
^22 “ which does not alter the result for the mixing angle to leading 
order, see Appendix C of [205], Moreover, the off-diagonals can be con¬ 
strained as I ifigg I = I ifigg |, leading to the requirement for the mass matrix 












6.1 HIERARCHICAL MASS MATRICES 


101 


to be normal. Only the phases are unconstrained, so we impose for the 
2 x 2 submatrix 


mf = 


f 0 ^ 

V|m4|e^^32 


m4|e^^23A 

J 


(6.5) 


We then reparametrize the mixing matrix (i. e. the left-singular matrix— 
right-singular matrices are unobservable in weak charged current interac¬ 
tions) via the two invariants of the Hermitian product ', 

f f ^ f ^ l!?l^ I f 

Tr w = m 2 + ^3 =2 + ^33 , 

r f2 f 2 I f \4 

detm = m 3 = m 23 > 
and find 




and 



1 -m^. 


( 6 . 6 ) 


The mixing angle is found to be tan 9^^ = ^and the 
rotation matrix can be expressed in terms of this angle and one complex 
phase (other phases are unphysical and can be rotated away) 





(6.7) 


With Eq. (6.7), we have the left-singular matrix for an /-type fermion. 
The weak mixing matrix is composed out of up- and down-type mixing, 
^23 = ^ 23 ^ (^ 23 ^) ' • two-generation case we have 


V 23 = diag(l,e '^ 23 j 



Ce-" A 
j 


diag(l,e'^23j . 

( 6 . 8 ) 


The phase Sgg factored out can be absorbed in a global rephasing of third 
generation quarks. We keep the phases 5 and Sq to analyze their origin 
and to keep track of the phases of the individual rotations, so dgg and 
for later purpose. The parameters of Eq. ( 6 . 8 ) are found to be 




m^ + m^ - 2 ^m“m^ cos (d^g - 
(l + m-)(l + m^) 


/fig sin 1 

[5; 

^ - Sis) 

1 

/hg — /fig COS 1 

[^23 “ ^ 

i] 


sin ( 5 ^ 3 -5^) 

1 + cos ( 5^3 - 5 ^ 3 ) 


sin 023 = 


\ 


(6.9a) 


(6.9b) 


(6.9c) 


Expanding in the small mass ratios ffi^ we may obtain tanSg 0 and 
tan 5 tan ( 5 ^ 3 - 5 ^ 3 ). 


A normal matrix 
obeys 

mm' = m'm. 


See [333, 335, 345, 
357]. 


This can be seen 
from Eq. (2.9) with 
S“ = and 
s^ = (s«)Vt. 

It is important to 
note, that two 
unitary mixing 
matrices do not 
commute as real 
orthogonal 2x2 
matrices do, and 
the new mixing 
angle is not just 

®23 = ®23 ^ ®23 

well as the phase 

^ ^ ^23 ^ ^23- 
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6.2 FULL HIERARCHY AND THE NEED FOR CORRECTIONS 


Each individual 
rotation is given by 
S{9,5) = 
-15 


cq sge 
\5 


-Sge 


C9 


We have to keep in 
mind, that the 
object of interest is 


Comparing the 
masses, we 
conclude that 


O 




The description can 
also be applied to 
the 1-3 mixing with 


■J 


■ ifi{ and 

^ H)'- 


We construct the full hierarchy and full-rank picture following the con¬ 
secutive breakdown of U(3) symmetries. The two-generation description 
from Sec. 6.1 can be easily generalized to two-flavor mixings within three 
generations by filling up 3 x 3 matrices with zeros where appropriate. In 
view of minimal flavor symmetry breaking, the 2-3 mixing acts “first” on 
the mass matrix. Therefore, 

^L ~ ^12^ (®12’^12) ^13^ (^13’^13) ^23 (®23’^L)- 
The weak mixing matrices are then combined as 


( 6 . 10 ) 


V, 


CKM 


_ 5L, u^L, u^L, uZ^L, dATZ-L, dAlY_L, dAT , 

“ '^13 '^23 V'^23 J V'^13 J V'^12 J ’ 

L, V A1 


jj — cLj e Yck vA Y^L, vW ZpL, vA 

L'cKM - *12 *13 *23 1*23 J 1*13 J 1*12 J 


Building up the mixing matrices from lower-rank approximated mass ma¬ 
trices, we have to take care of all equally large contributions. The leading 
order 2-3 rotation diagonalizes the Hermitian product of a hierarchi¬ 
cal mass matrix of type (6.4) 

(o " " ^ 
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0 X" 
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0 X 1-X' 
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0 0 
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0 0 0 
0 0 1 


+ o (x^), 


where X is O fh^^. Now, the transition to the full-rank matrix comes 

along with new contributions O everywhere: in general, the lower- 

rank approximated matrix differs in all elements from the higher-rank 
one, 

f. . A 


with * = O 



The inclusion of the smallest mass leads also to nonzero matrix elements 
where previously neglected ~ X^ terms were present. For a reasonable 
description of the mixing matrices including all missing pieces already in 

the 2-3 rotation, O ~X^ as well as O ~X^. We include 

those “corrections” by correcting rotations like the composition 
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where the composed rotation angle is with e = m{, 
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For small S the mixing angle only changes slightly, as expected for a pertur¬ 
bation. This procedure appears to be equivalent to the adding of a much 
more complicated term in the leading order mass matrix: 




1 + 




(1-e). 


Collecting all corrections, we then find 


(2) ^^L,/ 


(m{ mQ [mQ . 


( 6 . 11 a) 


A similar discussion now holds for the succeeding 1-3 rotation, where 


o(m{) 


and the O 


eluded: 


[m{mQ 


rotation has also to be in- 


L.f (2) _ L,f 
■^13 '^13 


(m{ ihQ S\’/ S\’/ (m{) . 

Finally, the 1-2 rotation can be exactly solved and we have 


cC / _ f 

■^12 ~‘^12 




( 6 . 11 b) 


( 6 . 11 c) 


The only physical relevant phase can sit in the 1-2 rotation: following 
the rank evolution of mass matrices, where the rank-one mass matrix has 
a U(2) symmetry left in the 1-2 block. The other phases (which were 
omitted in Eqs. (6.11) are either zero or maximal (ti) corresponding to 
a rotation either in the same or the opposite direction than the leading 
order rotation so that there is no CP violation connected to them. 


We keep in mind 
that we actually 
describe rotations 
in three dimensions, 
so all matrices shall 
be 3 X 3 matrices 
where the 
nontrivial 
(i. e. nonzero and 
^1) entries are 
distributed to the 
right positions. 
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We motivated a description of individual rotations parametrized by mass 
ratios as physical parameters only. Fermion masses are the singular values 
of mass matrices and eigenvalues of the Hermitian product. We cannot 
set the overall scale (say the largest mass value) since a scale factor can 
always be multiplied out. The CP-phase in the weak mixing matrix has 
to be constructed in a similar manner out of the mass ratios. However, 
the origin of CP violation stays unclear. In principle, there occurs at each 
fundamental rotation one complex phase which we take either zero or n 
except for the 1-2 rotation. In any case, the phases always appear in pairs 
as in Eqs. (6.9) and therefore we only keep track of one of them. 

The existence of CP violation in quark flavor physics enforces at least 
one complex parameter. As we argued above, only 5^2 can carry any infor¬ 
mation about CP nonconservation, so we take it maximally CP violating: 
5 12 = f • This somewhat arbitrary choice gets obvious looking into the 
“data”, see Pig. 25. In any case, we do not want to take 612 or even as 
free parameter. We can show that it is not only possible to reproduce the 
mixing angles (or CKM elements) but also the CP-phase (or Jarlskog in¬ 
variant) with fermion masses only. A similar line of arguments is followed 
in [389, 390]. 


It is interesting to 
note that we want 
to describe four 
mixing parameters 
in the CKM and 
PMNS matrix (three 
angles and one 
phase) by four mass 

^ U.V u.v 

ratios, , m 2 , 

^ d.e _ j w d,e 

ana m 2 • 
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Figure 25: We show the upper right corner of the CKM matrix (V^f, over left and over right) 
as it results from Eqs. (6.12a) and (6.12b) with values for the masses plugged in. The small 
red dots scatter over the 1 <j regimes for the measured quark masses, where the blue crosses 
show the results from the central values. All phases in Eq. (6.12b) were allowed to take one 
of the values {0, ti}. 


We deconstruct the CKM and PMNS matrices as Vckm = (^l) ' 

t^PMNS = (^l) 


S“ = s\’^ ( — IS 



(6.12a) 


(6.12b) 
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Table 2: We show the input values of quark masses and their values at the Z scale 
Q = by virtue of the RunDec package [391]. The mass inputs corre¬ 
spond to the experimentally measured values while the outputs, evalu¬ 
ated at the Z pole, include the resummation of higher order corrections 
from QCD by the RG running. RunDec takes properly into account the 
decoupling of heavy quarks below their scale. Charm and bottom quark 
can also be simultaneously decoupled [392]. All masses are given in 
GeV. 


Input [23, PDG 2012] Output 


m,(2GeV) = 0.0023+°;“°J 

muiMz) 

= 0-0013+o.ooo3 

ma(2GeV) = 0.0048+°;°°ol 

ma{Mz) 

= 0.0028+°X' 

ms (2 GeV) = 0.095 ± 0.005 

ms{Mz) 

= 0.055 ±0.003 

me(mc) = 1.275 ±0.025 

m^(Mz) 

= 0.622 ±0.012 

mb(mb) =4.18 ±0.03 

m^{Mz) 

= 2.85 ±0.02 

mt(OS) = 173.07 ±1.24 

m,{Mz) 

= 172.16+[,\^ 


6.3.1 Quark mixing 


With the previous work, we can plug in the numbers for the quark masses 
as given in Tab. 2 and propagate the errors (the errors given in Eq. (6.13) 
are seen to be purely parametrical). The results have to be compared with 
the global fit values for CKM elements as shown in Eq. (2.11) and are 
found to be in an astonishingly good agreement within the errors: 



^ 0 Q74+0-004 

r\ oo[--L0.016 

O.^^u_0,oil 

0 0031^'^ °°^®^ 

U.UUOl_o 0015 


„prop 1 _ 

n O9ir“l“0-016 
'^*'^'^^-0.011 

0 974+0-004 

U.UO^^_0 004 

. (6.13) 


^U.UUo7_Q QQQg 

0.038+o°o°o4 

n 9997+0.0002 
^•^^^^-0.0001 J 



Our choice of 
phases in 
Eq. (6.12b).- 

e _ n 

Si2 - 2’ 

( 0 ) 


^13 
g(l) 

5 ( 2 ) 

g(0) 

"23 

g(l) 

‘^23 

g(2) 

‘^23 


= 0 , 
= rc, 
= ri, 
= 0 , 
= ri, 
= n. 


6.3.2 Lepton mixing 


The lepton case is a bit more involved because we only know the masses 
partially. However, we can revert the procedure and predict the neutrino 
mass spectrum out of the measured mixing matrix. Eor the 1-2 mixing we 
have 


+ m^i2 - 2^me^,m^i2Cos(5^2 “ ^ 12 ) 


(1 + me|j^)(l + myi2) 


\Ue2\^ 


(6.14) 
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mv/eV |Ue2] 




Figure 26: For low values of the neutrino mass spectrum is hierarchical in contrast to quasi¬ 
degenerate as shown on the left. The right plot shows the dependence of |t/e 2 l on the phase 
The experimentally allowed 3 ct ranges are indicated by the horizontal red lines. The 
best fit point lies surprisingly close to 5]'^ = ^ . Figure taken from [205]. 


with = m^/m^ and myi 2 = rtiy^/niy^. The neutrino spectrum can be 
obtained inverting Eq. (6.14) to get a solution myi 2 {U^ 2 ) ^nd therewith 



yAm2 ^/(l-m2i2), 


^vi = 1 

Jml^ - Aml^, 

(6.15) 

1 




We use Trig = 
mg/m^ = 0.00474 
and 

\Ug2\ = sin 012 = 

0 . 54 ... 0 . 56 . 


With the analogous choice (motivated by the observation of Fig. 26) of 
the 1-2 phase as in the CKM case, 5^2 ~ ^12 ~ f > 




|ffe2l^(l + hte) -hte 
l-|[/e2p(l+me) 


0.41... 0.45 


(6.16) 


Errors are added 
linearly to be 
conservative. 
However, this 
neutrino mass 
spectrum is clearly 
non-degenerate. In 
this way, the results 
of this chapter 
might be of 
importance if there 
is no neutrino mass 
measurement in the 
near future as 
desired for the 
corrections of 
Sec. 3.2. 


and estimate the neutrino mass spectrum to be 


ruy^ = (0.0041 ± 0.0015) eV, 
ruy^ = (0.0096 ± 0.0005) eV, 
TUy^ = (0.050 ± 0.001) eV. 


It is interesting to note that the sum of all light neutrino masses is ^ m,, = 
0.0637 ± 0.003 and thereby perfectly fine with cosmology, see [199]. 

The mass hierarchy in neutrinos, though existent, is not as strong as in 
the charged fermion case. However, the result of our description relying 
only on the minimal breaking of a maximal flavor symmetry shows also 
for the PMNS matrix an astonishingly good agreement with the global fit 
data shown in Eq. (2.44) 


\U 


prop I 
PMNS I 


0-83^o°.05 

0 54+0-06 

0.14±0.03 

0-38^o°o6 

0 57+0-03 
^•^^-0.04 

0.73 ±0.02 

0-41+0.06 

^•”-‘■-0.04 

0.67 ±0.02 


(6.17) 
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where we assigned the phases of Eq. (6.12d) to be 5^2 = §> ^ 13 ^ = 
^ 13 ^ = 71, 5 ^ 3 ^ = 71, 5^23 = ^ 23 ^ = ^ ^^id 5^^ = 0. We remark an 

exchange of «—> compared to the CKM case. 

Similar results for the neutrino masses were found by [200, 201, 360] 
for hierarchical charged lepton masses. However, the third mixing angle 
0^3 was predicted too low with 3°. 

6.3.3 CP violation 

A measure of CP violation in fermion mixing is given by the Jarlskog in¬ 
variant [393] 

J = lm{VijVklVX^?j. (6.18) 

A vanishing J indicates CP conservation. With the decomposition of mix¬ 
ing matrix elements in terms of mass ratios, we can give an approximate 
analytic formula for J for both quark (f = q and a = u, b = d) and lepton 
if = I and a = e, b = v) mixing [205] 

Jf cos Oj ’2 sin 9^2 sin @23 ®12 ^23 “ 1 “ ^13 “ ^ 12 ) > 

where 


sin 04 = 


y. 
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/12 
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Qa(b) _ 


13 




m 


a{b) 

12 
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i(fa) 


and 


12 


sin 0 


i(b) 


±\ m 


a{b) 

12 


^ a(b) „ a(b) 

^13 ^23 
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i{b) 


23 


13 


(i+-r)(i+ 

For the quark case, we have 


„ a(b) „ a(b) 
^13 ^23 


)(i+(<‘’)d 


J, = Im (v'„Acl.L*l,''c0 = x 1“"'. 

which has to be compared with Jg = (3.06 ^q'2o) x 10“^ [23] and due to 
our large error is found to be in agreement. In the lepton sector there has 
no CP violation yet been observed, so only an upper bound on the Jarlskog 
invariant can be set, = 0.033 ± 0.010 [132]. We find a result close 
to this maximal value, 

J, = Im [u, 2 U^ 2 U: 3 U; 2 ) = 0-031+°;“7^ 

This corresponds to a “prediction” of d^p 70° from the central values 
via Jf = J“^^sindcp. In any case, we are compatible with maximal CP 
violation (d^p = 90°) as well as d^p = 45° taking the lower error limit. 


The deeper reason 
behind this 
replacement has 
still to be found out. 


with fftT = m^ 
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The neutrino mass 
ratios, which 
govern the results 
for the mixing 
angles, are large! 

SUMMARY OF CHAPTER 6 We presented a formulation of fermion mix¬ 
ing based on the observation that in absence of Yukawa couplings, the SM 
fermion content intrinsically has an enhanced symmetry group, [U(3)]^. 
Breaking this maximal flavor symmetry in a minimal way by only intro¬ 
ducing the masses (diagonal Yukawa couplings), we could demonstrate 
that with hierarchical fermion masses as they are observed in nature, we 
are able to reproduce the mixing patterns for both quark and lepton mix¬ 
ing without further assumptions. Especially no specific patterns for the 
fermion mass matrices (“texture zeros”) are needed. The rather mild hi¬ 
erarchy for neutrinos, whose mass spectrum has been predicted by this 
minimal breaking of maximal flavor symmetry, gives a handle on large 
neutrino mixing. Moreover, the strong connection of fermion masses and 
mixing allows to predict the neutrino mass spectrum, taking neutrino mix¬ 
ing and the charged lepton masses as input. We have found a lightest 
neutrino mass m^y^ = 4.1 meV which is well below any direct detection 
possibility. 

The derived mixing matrices are in surprisingly good agreement with 
the experimentally measured flavor mixing. By imposing only one non¬ 
trivial CP-phase in the fundamental rotations, we are also in agreement 
with the effective CP-phase as determined via the Jarlskog invariant. It 
is a crucial point for further investigations, that 5-^2 = ^ and all other 
phases are either zero or n. This is also interesting in a model building 
perspective: we do not have to predict an arbitrary, continuous CP-phase 
but can choose phases of the individual rotations out of discrete values 

There can be a lot of work done with the results of this chapter in mind: 
what is the origin of CP violation? Is there maybe a discrete symmetry 
behind yielding discrete 5;^? How can we model the successive symmetry 
breakdown? Is there a SUSY description inspired by REV possible? Each 
symmetry breaking step occurs with the introduction of one Yukawa cou¬ 
pling for the corresponding fermion generation. Can we combine flavor 
breaking, electroweak breaking and SUSY breaking? 



CONCLUSIONS 


Ich kann's bis heute nicht verwinden, 
deshalb erzahl’ ich’s auch nicht gem: 
den Stein der Weisen wolit’ ich finden 
und fand nicht mal des Pudels Kern. 

—Heinz Erhardt 

We have addressed the influence of quantum corrections to neutrino 
masses and mixings and found enhanced corrections in case of a quasi¬ 
degenerate mass spectrum that have the power to completely generate 
the flavor mixing irrespective of any tree-level flavor model. Especially the 
case with trivial mixing at tree-level connected to exact degenerate neu¬ 
trino masses was shovra to give the observed deviations from the degener¬ 
ate spectrum and simultaneously generate large mixing angles. We have 
calculated threshold corrections in an extension of the Minimal Super- 
symmetric Standard Model including right-handed Majorana neutrinos in 
Chapter 3. The existence of right-handed neutrinos and Majorana masses 
at a very high scale can be motivated from grand unified theories. We have 
analyzed the impact of flavor non-universal SUSY breaking terms and pro¬ 
posed a description of neutrino mixing in terms of SUSY breaking param¬ 
eters in the context of the MSSM, extending previous ideas on REV in the 
literature. The driving force is given by contributions ^ij / ^susY that 
do not decouple with the SUSY scale. Even with very heavy SUSY spec¬ 
tra, this non-decoupling contribution to rather low-energy flavor physics 
persists. Additionally, we have shovra that the presence of hierarchical 
right-handed (s)neutrinos alone can substantially alter mixing patterns— 
a result which is also in line with what was found for the Standard Model 
with right-handed neutrinos before. We have performed a Dyson resum¬ 
mation of the enhanced contributions to neutrino mixing and found a 
stabilization of the description with respect to the neutrino mass. The 
same combination AT / Mgusy was found to give an equally good descrip¬ 
tion of neutrino mixing for the full neutrino mass range: This is a kind 
of “non-decoupling” contribution with respect to the neutrino mass. With¬ 
out resummation, smaller corrections are needed to give the same mixing. 
Once resummation is switched on, the same parameters give an equally 
well generation of neutrino mixing irrespective of the neutrino mass. 

SUSY threshold corrections also affect the effective Higgs potential and 
alter conclusions about electroweak symmetry breaking. We have found 
and described a new class of vacua in Chapter 4, where the term vacuum 
describes the ground state of the theory or a minimum of the effective 
potential. Einding the global minimum is a challenging task. However, we 
discovered regions in the SUSY parameter space that treat the electroweak 


I have not yet 
overcome, that’s 
why I don’t like to 
tell it gladly: I 
wanted to find the 
philosopher’s stone, 
and I haven’t even 
found the core of 
the poodle. 
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minimum with v = 246 GeV as false vacuum. These findings are essential 
to constrain parameter regions. Transitions to the deeper minimum were 
found to happen nearly instantaneously, so those constraints can be seen 
as strict without referring to metastability and the life-time of our universe. 
Moreover, the results hint towards a charge and color breaking global min¬ 
imum and therefore extend existing bounds at tree-level. The full analysis 
including colored directions in field space has not yet been performed and 
is left for future work. 

We have combined the discussion of vacuum stability in the MSSM with 
the right-handed neutrino extension in chapter 5. Without Supersymme¬ 
try, the existence of the effective neutrino mass operator as well as the 
extension with singlet neutrinos at a high scale render the Higgs poten¬ 
tial unbounded from below—a behavior that better shall be avoided. The 
supersymmetric version, however, rescues the potential and also does not 
induce further instabilities neither at the SUSY scale nor at the scale of 
heavy neutrinos. At the high scale, SUSY appears to be effectively exact 
and SUSY masses are only a small perturbation; so the contributions to 
the effective potential exactly cancel summing up fermionic (neutrino) 
and bosonic (sneutrino) contributions. 

Finally, we have elaborated on a different view of fermion mixing and 
how to parametrize the mixing angles in terms of mass ratios only. We 
found a minimal breaking of the maximal flavor symmetry group [U(3)]^ 
together with hierarchical masses sufficient to reproduce the observed 
mixing patterns for quarks as well as leptons. Moreover, by inverting the 
procedure, we predict the lightest neutrino mass to be = (4.1 ± 
1.5) meV. The large error is result of the still large uncertainties in neu¬ 
trino physics (Am^ as well as the PMNS matrix itself). By the same de¬ 
scription, we also predict the CP violating phase in lepton physics to be 
rather large (close to maximal, ^). 

The main results of this thesis are a description of neutrino mixing at 
quantum level with the restriction to a quasi-degenerate physical mass 
spectrum—and genuine one-loop bounds from vacuum stability on SUSY 
parameters. Also the description of fermion mixing angles in terms of mass 
ratios may be lined with a radiative origin of the minimal flavor breaking 
chain, which is kept for future work. 


TECHNICALITIES 


A.l SPINOR NOTATION AND CHARGE CONJUGATION 


We collect the necessary notations and rules to deal with left-handed Weyl 
spinors and charge conjugation. The Weyl spinor notation turns out to be 
more convenient. A four-component Dirac spinor lE is written in the chiral 
basis with left- and right-handed Weyl spinors ip^ and t/ir as 


xp = 



(A.1) 


We use an index free notation (there is no need to introduce dotted indices 
for our purpose and write a bit sloppy as a column). The left- and 
right-handed fields can be projected out with Pl = | (1 — Ts) and Pr = 

I (1 + Ts)^ 


>El = ^ 




Charge conjugation is defined as 


C: ip ^ ip" = C>P^, (A.2) 

where C = iy^y® with C = = C~^ = —C and ^ = 'P'y°. The charge 

conjugation of a left-handed field gives a right-handed one and vice versa: 


C : *Pl ^ (<Pl)^ = (>P‘^)r = <P^, vpR ^ ('Pr)^ = (*P^)l = 'P^ (A.3) 

It follows clearly that 

= c<P^ = I 1, 

Uu 

and therewith (j/ir)'^ — which is the left-handed component of the 
charge conjugated Dirac spinor. 

Dirac mass terms are = 'i/'r'i/ir + 'ipi.'ipRy Majorana mass terms 'P^ClP = 
i’hi’L + <^with t/ir = xpl). 


a.2 NEUTRALINO and charging mass AND MIXING MATRICES 

The chargino and neutralino mass matrices follow basically from the soft 
breaking Lagrangian Eq. (2.28), the superpotential Eq. (2.26) which con¬ 
tributes the pi-term • Hy, and the gauge interaction (after spontaneous 
symmetry breaking, so with ^ d). 



A good compilation 
of two- and 
four-spinor notation 
and its applications 
can be found 
in [394], 


We have 

• n 1 9 

Ts = irr r r 

with the Dirac 
y-matrices that 
satisfy 

{y^,y''} = 2gf^r 


By abuse of 
notation we write 
i'^)x = with 

a: = l,r. 


We decided to write 
(V'r)" = i/'r in the 
definition of the 
MSSM superfields 
in Chapter 2 and 
label with r the 
gauge 

representation 
there and wherever 
applied. 


Ill 
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The Lagrangian part defining the mass basis for charginos is given by 

= + h. c. ) 

- (M 2 A+A“ + /ih+h” + h. c. ) CA.4) 

= - (^ip~YMc'4’^ + h. c. , 

where xp^ = (A+,/t+)^ and xp~ = (X~,h~y. The mass matrix 

f M 2 •\/2My^f sin 13 

\^v2MiyCos/3 /r 

can be diagonalized with a bi-unitary transformation 

= Mq = diagonal. (A. 6 ) 

The charged mass eigenstates (“charginos”) are defined as 

X+ = Zl[X+,h+y and I-= Z! 

Analogously, we have for the neutral electroweakinos and higgsinos 

^mass = -y ^3 (vd^d “ + y ^0 (vd^d “ ^u^u) 

+ — -M 2 A 3 A 3 — -M^XqXq + h. c. (A. 7) 

=-l{ryM^xp^+h.c., 

where xp^ = (Ag, A 3 ,/i°,/i°). The mass matrix 

Ml 0 —CpSy^Mz CpSyyMz 

0 M 2 CpCyyMz —SpCyy 

— CpSyyMz CpCyyMz 0 /I 

SpSwM^ —SpCy/Mz /I 0 

is Takagi-diagonalized with a unitary matrix 

Mlf ® = zlMnZn, (A.9) 

such that 

Xo = Z*{Xo,X3,h%~K) 

are the neutral mass eigenstates (“neutralinos”). 
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A.3 SFERMION MASS AND MIXING MATRICES 


UP AND DOWN SFERMION MASS MATRICES IN THE MSSM We restrict 
ourselves to the derivation of selectron and sneutrino mass matrices in the 
MSSM, where the squark sector follows analogously. The more involved 
situation with Majorana neutrinos is given below. We include right-handed 
neutrinos already without Majorana masses. 

The scalar mass Lagrangian is contained in the soft breaking Lagrangian 
and the F- and D-term potential, 


y'=yL 


+ ''/ + v/. 


We collect the individual contributions from Eqs. (2.28), (2.21) and (2.23), 


‘-soit ~ 


Ki ij + Sr,! ij 2rj + ij 

+ [^d • + h,i ■ + h. c. ] , 


v/ = 


yi = 


1 


dW 

J/7 


dcjii ^ 


j f*f 


1 

-D“D“ 

2 a a 


J/7 




(A.10) 
(A. 11) 


J /*/ 

(A.12) 


The exhausting part are the F -terms as derivatives of the superpotential. 
Sfermion mass terms are obtained after electroweak breaking, so each 
occurrence of ^ has to be replaced by its vev. 


Vi 


* ~ * TT-V* ~ 

= M Vd-VrJ^. Vrj 


+ y^n-~etiYtr^R,j 


i i 

+ S I + J] I VuY,]*VK,j 

i i 

mass ~ T7-e~* ~ i i_ 

= - FVdYij v^^iVRj - FVuYije^ ie^j + h. c. 


Mcot/Sv^m^VR 


/xtan^SejiTieep^ 


+ + ^XrVi^jT^j 






+ X<,rYtXXj +• 


(A.13) 




v^m^'m^VR 


The D-terms can be easily read off and are contributions ~ 2 




^gi (141^ - 141^) Xi - 2|eR^iP) 

i 

+ \sl {hd^ha + hlYK) £[ 


Mass terms are the 
bilinear terms in the 
Lagrangian, so only 
~ f*f are needed. 
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The objects are 
the third 
components of 
weak isospin, so the 
SU(2)l charge of 
sfermions f. 


where Tj/2, i = 1,2,3 are the generators of SU(2 )l. 

Combining to mass matrices, we get 

j^2 m2 + 2/31 + m^m®'!' -m^/icot^ + 

" V -m>*cot^ + vX m2+mfm^ j’ 

(A. 14a) 


X? = 


fhj + M|(r|L — QeSin^ 6^) cos2^1 + —me/itan/3 + v^A* 

^ — m^pi*tan/3 + V(jAg + M|Qe cos 2^ sin^ 0^1 + m^m 


e"‘ey 
r A 1 AU^ 


in a basis / = such that pMjf. 


SNEUTRINO SQUARED MASS MATRIX AND MIXING MATRIX Extending 
the MSSM by right-handed neutrinos and giving them a Majorana mass 
leads to a seesaw-like mechanism in the sneutrino sector. Similar to the 
seesaw-extended Standard Model, where the neutrino spectrum gets dou¬ 
bled, the sneutrino spectrum gets quadrupled. Why that? The MSSM con¬ 
tains only three sneutrino states. Including right-handed fields, the num¬ 
ber of states get doubled, although half of them are singlets under the SM 
gauge group. Moreover, due to Dirac and Majorana masses, the physical 
spectrum gets even more enlarged. Effectively, we are left with six light, 
more or less active states, and six heavy singlet-like states. A priori, the 
sneutrino squared mass matrix is therefore a 12 x 12 matrix, which can 
be perturbatively block-diagonalized similar to the neutrino mass matrix. 
The complete procedure is described in great detail by [140]. 

We choose the following basis: N = (vl, v*, v^, (such that —= 
N ' and classify chirality conserving iLL,RR') and chirality chang¬ 

ing blocks: 




with 




f 


V 




■^RR ~ 


( Ml, 


Ik«)' 



rfi^ + ^M|cos2/31 


■m>D' 

2 ^T 

I 


(m?)'^ + COS 2/31 


m>R 


-[JL cot pm 


-VnK 


-fjL cot/3m!;j’+ VuA^ 
(m2)T + m°‘m 
2B 


D* 

’’y 

dT d* 

}2 


m^Ml 


M*Mr 


A 


D* DT 


J 

(A. 15a) 
(A.15b) 


mt 


2{Bp 


(A. 15c) 
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A.4 FEYNMAN RULES FOR THE TYPE I SEESAW-EXTENDED MSSM 


The relevant vertices for the lepton flavor changing self energies are triple 
vertices for the lepton-slepton-gaugino and -higgsino interactions: 




• pi _ 

^ -t Vf — 




~ gl^lfc)^i^*(t^PMNs)i/] Pl (A.16a) 

+ Pk|, 




‘A" =j^\ (te24 + Si4t)»''-y54»'A-3,(AT6b) 


- ('^8i4i<+3,. +yn^?iK) ^4- 


= -1 


• r-> • 

if^,. 


'f 

iXk) 


e* 




's2Z^z^;p^-ytjZ^zlP^ 


(t^PMNs)i/^L> (A. 16c) 
(A.16d) 


where summation over double indices is understood. 

The vertices of Eqs. (A. 16) are given for an incoming standard model 
fermion, outgoing chargino or neutralino as well as sfermion. They gener- 
ically follow from an interaction Lagrangian like 


Ant=/i^|^'' Xkff + h. C. 


Each vertex comes along with the corresponding chirality projector: 


j-'fsXk _ j-'fsXk n _i_ r^fsXk n 

V; -hjfP^ + ^RjfP^- 

The mixing matrices diagonalize the mass matrices in the following 
manner: 

• Chargino mixing: z4A1c2_|_ = 

• Neutralino mixing: zJjA4nZn = 

• Slepton mixing: W^A4?We = 

• Sneutrino mixing: 

= Wlv^'MlVW^ = 

such that Zy = VWy diagonalizes the original mass matrix Ad? and 
therefore: 

= + and 

appear in the vertices of Eqs. (A. 16). 
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• Neutrino mixing: The PMNS mixing matrix can be determined from 
the neutrino mass matrix ihy = IfpMNs'^T^PMNS’ where rtiy is the 
effective light neutrino mass matrix and the charged lepton masses 
can be taken diagonal (otherwise there would be a contribution to 
the PMNS mixing similar to the CKM mixing from both up and down 
sector: l/pMNs = ^eL^T,L> where Vgp rotates the left-handed elec¬ 
tron fields). 


a.5 loop functions 


The dimensionality 
of space-time is 
“reduced” to 
D = 4-2£. 
Physical results are 
then obtained in the 
limit e —>■ 0. 


The evaluation of the one-loop potential as well as the supersymmet¬ 
ric threshold corrections to neutrino masses and mixing needs standard 
scalar loop integrals. In all cases, only solutions with either = 0 (ap¬ 
proximation for neutrino self-energies with heavy SUSY loops) or no exter¬ 
nal momenta at all (effective potential) are needed. The n-point integrals 
in the sea-urchin derivation of the effective potential can be reduced to 
derivatives of tadpole integrals. We follow the notation of [216]. 

In dimensional regularization [395] the tadpole integral is given by 


r 


Q 


d^k 


2s r 




1 




iO 

1 


(27r) + iO ’ 


(A. 17) 


with the arbitrary renormalization scale Q. The i 0 prescription evades the 
singularities on the real axis and is omitted in the following propagators. 
Evaluation of the D-dimensional integral then yields 


Ao(m) =m2 ( Ag-ln^ + 1 I +0(e), 


(A.18) 


Ye is the Euler 
constant, 
Ye^ 0.577. 


with Ag = ^ — Te + log 471. Where necessary, we work in the modified 
minimal subtraction (MS) scheme and subtract Ag (not only the pole -) 
with the counterterms. 

Similarly, the two-point integral is given by 


Bo(p^;mi,m2) = 
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1 71^ 
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The momentum independent three- and four-point functions are conse¬ 
quently 
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with Dj = k^ — mf. 
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We finally list the MS-subtracted expressions for the loop functions: 
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Results partially 
transferred from 
[89] and [85]. 


In Sec. 4.1.3, we calculate explicitly the 4-point function (and later n-point 
functions) for equal masses, llo(m, m, m, m). Obviously, there is a relation 
withAo(m) = sj 
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This easily can be generalized to an n-point scalar integral: 
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RELEVANT RENORMALIZATION GROUP EQUATIONS IN 
THE SM AND BEYOND 


B.l SM WITH HEAVY SINGLET NEUTRINOS 

Unfortunately, all RG including heavy Majorana singlets without SUSY are 
only available to one-loop order. In the MSSM, two-loop results exist. We 
nevertheless exploit the two-loop RGE for the SM (without heavy singlets) 
and add the neutrino part above the threshold to one-loop order. The 
RGE used for the study of the RGI potential are given by the following 
equations with t = ln(Q/ Qq) with some arbitrary but fixed scale Qq. 
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We are using the p 
function for the 
Higgs self-coupling 
as provided by 
[396] with 
h = l/16n:2, 
dj. =3 ,Ng= 3, 

C/ =4/3. 


54A2g2 + l8A2g2 - 156A3 - dj^Y^Vf - ^4Y2g4 + 

- ^d^Yfg^, - 7d«Y2AY2 + ^d^Y2^g2 + ^d^YfXg^ 

- 24d^YfX^ - d^YfY^ - ^-dj,Yfgl - \dj,YfX + Sd^Yf 


- 4Cfd,glYf + lOC^d^^g^AY^ + lOC^d^g^Y^A - 4C^d^g32Y4 


+ 


167r2 


A1;"-2Y; 


119 






120 


RELEVANT RENORMALIZATION GROUP EQUATIONS IN THE SM AND BEYOND 


d K 

d t 




+ A+ 


dMp 


eiQ-M^)3Y^ Jk, 
1 


d t 

1671^ 

dYy 

1 

d t 

1671^ 


- ^ 




2Y^"Mr0(Q-Mr), 




^2 ^Y 


dYe 1 /5 „ „ „ , 1 „ 

- - - + 3^t - e (Q - Mr) 


d t 1671 


4 S 1 - 4 S 2 


d t 


9 3 


16^2 I 2^b + 2 V + + e (Q - M^)Y^ 


1 2 9 2 „ 2 ) 

-4gi-4g2-8g3 - 


d t 


+ 0 (Q - MR)y, 


dgi 
d t 
dg2 
d t 
dg3 
d t 


17 2 9 2 2 ) 

-^gi-4g2-8g3j. 

1 41 3 

■ 167i2 10^^’ 

1 19 3 

^2- 


1671^ 16 
1 


1671^ 


7g3- 


(B.lb) 


(B.lc) 

(B.ld) 


(B.le) 


(B.lf) 


(B.lg) 


(B.lh) 

(B.li) 

(B.lj) 



B.2 MSSM WITH HEAVY SINGLET NEUTRINOS 


121 


B.2 MSSM WITH HEAVY SINGLET NEUTRINOS 


Singlet neutrinos do not influence the gauge couplings at one loop. How¬ 
ever, there is a significant impact for the two-loop RGE [397]: 
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Two-loop RGE for neutrino parameters in various seesaw models can be 
found in [165], especially the equations for the MSSM extended by singlet 
neutrinos: 
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